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I.  INTRODUCTION 

A.  BACKGROUND 

The  MEO  UHF  Satellite  Constellation  for  Robust  Assured  Telecommunications 
(MUSCRAT)  communications  system  design  described  in  this  chapter  was  created  as  a 
project  for  AA4871,  Spacecraft  Design  II  at  the  Naval  Postgraduate  School.  The  project 
provided  experience  in  the  design  of  a  complex  satellite  at  both  the  system  and  subsystem 
levels.  Each  student  in  the  project  team  was  assigned  responsibility  for  a  major  subsystem 
or  design  support  function.  Design  and  analysis  of  each  subsystem  was  performed  over  a 
period  of  eleven  weeks  using  facilities  available  at  the  Naval  Postgraduate  School. 
Additional  support  was  obtained  from  government  and  commercial  sources  as 
documented  in  this  report.  [Ref.  40] 

The  MUSCRAT  project  was  sponsored  by  the  Naval  Space  Command.  The  goal 
of  this  project  was  to  design  a  spacecraft  bus  and  payload  to  meet  the  essential 
requirements  for  a  Medium  Earth  Orbit  (MEO)  satellite  communications  system  which 
will  replace  the  aging  UHF  Follow-on  (UFO)  Constellation.  Supporting  trade  studies, 
analyses,  and  options  for  increasing  spacecraft  capability  are  presented  along  with 
associated  impacts  on  spacecraft  mass,  power,  cost  and  schedule.  [Ref.  40] 

1.  Mission  Overview  and  Requirements 

Advances  in  information  technology  place  us  at  the  point  of  a  revolution  in  military 
operations  -  one  in  which  the  ability  to  collect,  process,  and  disseminate  information 
causes  new  thinking  in  doctrine,  strategy,  tactics,  and  procedures.  Information  and  the 
dominant  battlespace  knowledge  that  it  imparts  allows  a  smaller,  mobile  joint  force 
structure  to  be  increasingly  lethal  and  survivable  in  the  battlespace.  Success  in  modern 
warfare  depends  heavily  upon  achieving  information  dominance  and  that  dominance 
depends  upon  the  connectivity  provided  by  communications  systems.  [Ref.  40] 

There  is  a  compelling  need  for  space-based  communications  assets  to  help  meet 
the  growing  information  needs  of  the  advanced  warfighting  and  supporting  systems  the 


Department  of  Defense  (DoD)  is  investing  in  for  today  and  the  future.  Although  land,  sea, 
and  air-based  communications  systems  are  vital  to  the  warfighter,  they  cannot  meet  all  of 
the  warfighters'  information  needs.  Space  based  communication  assets,  due  to  their  instant 
accessibility,  survivability,  coverage,  flexibility,  and  ability  to  support  mobile  subscribers, 
are  in  some  cases  the  only  practical  means  to  support  mobile  warfighter  platforms,  such 
as  ships,  submarines,  aircraft,  and  land  combat  vehicles.  As  a  result,  the  demands  on 
space-based  communications  are  growing.  Recent  DoD  studies  and  industry  assessments 
project  at  least  a  three  to  five-fold  growth  in  warfighter  information  requirements  that 
must  be  satisfied  by  SATCOM  systems  in  the  next  decade.  [Ref.  40] 

Under  the  direction  of  the  Deputy  Under  Secretary  of  Defense  for  Space 
(DUSD(Space)),  a  Mobile  Users  Study  (MUS)  was  commissioned  on  19  February  1997 
to  define  the  requirements,  architecture,  and  cost  for  a  replacement  to  UFO.  The  Naval 
Space  Command  was  tasked  as  the  lead  for  the  Navy-directed  MUS  Requirements 
Working  Integrated  Product  Team  (WIPT).  This  product  team,  consisting  of  52  members 
from  all  military  staffs,  the  Joint  Chiefs  of  Staff  (JCS-J65),  the  National  Reconnaissance 
Office  (NRO),  and  the  Defense  Information  Systems  Agency  (DISA),  reviewed  over  100 
general  requirements  for  low  data  rate  netted  communications  and  narrowed  them  into  a 
basic  set  of  8  high  level  needs: 

a.  Assured  Access 

Assured  access  to  SATCOM  services  is  the  most  fundamental  SATCOM 
need  of  the  warfighter.  The  SATCOM  study  process  confirmed  that  deployed  and  mobile 
warfighters  are  largely  dependent  on  SATCOM  to  satisfy  their  most  critical  beyond-line- 
of-sight  information  transfer  needs.  Unified  Combatant  Commanders  have  repeatedly 
requested  the  ability  to  access  SATCOM  on  demand  and  to  control  the  resources 
apportioned  to  them  by  the  Joint  Chiefs  of  Staff  (JCS).  A  warfighter's  access  to 
SATCOM  support  must  therefore,  be  available  on-demand  when  and  where  needed  for 
the  duration  of  the  mission.  [Ref.  40] 


b.  Netted  Communications 

A  communications  topology  should  be  provided  that  allows  multiple  users 
to  transmit  (one  at  a  time)  and  receive  (all  users  less  the  one  transmitting)  with  others  on  a 
common  data  stream.  Netted  communications  on  demand  exists  when  flexible  and 
reconfigurable  nets  enable  the  net  control  to  bring  on  new  members  to  the  net  without 
difficulty.  [Ref.  40] 

c.  Communications  on  the  Move 

This  capability  shall  provide  the  ability  of  the  warfighter  to  move,  and  talk, 
at  the  same  time.  This  capability  enables  the  user  to  have  voice  and  data  communications 
in  any  natural  environment  to  include  a  double  canopy  /  jungle,  in  rain  at  rain  rate  mode 
"H",  and  in  a  sea  environment.  [Ref.  40] 

d.  Jointly  Interoperable  Communications 

Current  and  emerging  warfighting  doctrine  require  joint  operations  to 
include  inter-service  compatibility  and  operating  with  allies.  Joint  service  and  allied 
operations  require  jointly  interoperable  communications.  [Ref.  40] 

e.  Worldwide  Coverage 

Worldwide  coverage  is  defined  as24  hours  a  day  communications  service 
from  65  degrees  N  to  65  degrees  S  latitude  without  gaps  in  geographical  coverage.  [Ref. 
40] 

f.  Point  to  Point  Communications 

The  system  shall  provide  a  communications  topology  that  allows  a  single 
terminal  to  connect  to  another  single  terminal.  Examples  are  communications  link  between 
a  commander  and  a  single  subordinate  element,  a  control  platform  or  a  remote  device,  or 
even  a  normal  telephone  connection.  [Ref.  40] 


g.  Broadcast  Communications 

The  system  shall  provide  a  communications  topology  that  allows  a  single 
user  to  transmit  while  multiple  users  receive.  A  "one  to  many"  broadcast  capability  is 
particularly  valuable  for  the  transmission  of  intelligence  products  or  any  information  that 
needs  to  be  distributed  from  a  single  source  to  multiple  users.  [Ref.  40] 

h.  Polar  Coverage 

Polar  coverage  is  defined  as  24  hours  a  day  communications  service  for  the 
area  of  the  Earth's  surface  that  is  above  65°  N  and  below  65°  S  latitude.  Communications 
service  coverage  primarily  directed  at  the  northern  polar  regions  is  important  for  reasons 
of  national  security.  [Ref.  40] 

2.  Management  Responsibilities 

Management  functions  and  responsibilities  for  MUSCRAT  will  mirror  those 
assigned  to  current  Military  Satellite  Communications  (MILSATCOM)  systems  as  set 
forth  in  CJCS  MOP-37.  Figure  1-1  illustrates  the  organizational  structure  for 
MUSCRAT.  Table  1-1  outlines  the  responsibilities  of  each  of  the  organizations  involved. 
[Ref.  40] 


Users 
(Components,  JTFs) 


COMNAVSPACECOM 
System  Operational  Manager 


Satellite 

Pay  load 

Health 

Control 

AF< 

NCTAMS/ 
JFTOC 

Figure  1-1 :  Organization  Structure  for  MUSCRAT 


Organization 

Responsibilities 

Joint  Staff 

1 .  Approve  requirements  as  part  of  the  Integrated  Communications  Data 
Base  (ICDB)  process 

2.  Apportion  satellite  resources  to  CINCs 

3.  Adjudicate  resource  apportionment  conflicts  not  resolved  at  lower  levels 

4.  Approve  the  Operational  Requirements  Document  (ORD) 

Defense  Information 
Systems  Agency  (DISA) 

1.     Administer  the  ICDB 

Unified  Commanders-in- 
Chief  (dNCs) 

1 .  Review,  validate,  and  forward  requirements  to  the  Joint  Staff 

2.  Collect,  validate,  and  prioritize  satellite  access  requests  from  their 
components 

3.  Manage  satellite  resources  apportioned  to  them  by  the  Joint  Staff 

Commander-in-Chief, 
United  States  Space 
Command 
(USCINCSPACE) 

4.     Provide  management  oversight  of  communications  payload  control  at 
the  direction  of  the  Joint  Staff  and  in  coordination  with  the  System 
Manager 

Chief  of  Naval  Operations 
(CNO  N63) 

1 .     Carry  out  the  Department  of  the  Navy's  responsibilities  as  the  System 
Manager  for  MUSCRAT 

Commander,  Naval  Space 

Command 

(COMNAVSPACECOM) 

1 .  COMNAVSPACECOM  is  MUSCRAT's  system  operational  manager 

2.  Develop  and  coordinate  procedures  for  payload  control 

3.  Manage  communications  payload  configuration  to  support  CINC 
requirements 

4.  Review  requirements  and  perform  analyses  to  determine  feasibility  of 
satisfaction  via  the  MUSCRAT  system 

5.  Recommend  apportionments  to  the  Joint  Staff 

6.  Provide  system  use  and  status  information  to  the  Joint  Staff,  DISA,  and 
USCINCSPACE 

Naval  Satellite  Operations 
Center  (NAVSOC) 

1.  Coordinate  MUSCRAT  payload  health  and  welfare  with  the  AFSCN 
ground  stations  and  perform  payload  control  functions 

2.  Perform  anomaly  detection  and  resolution  in  conjunction  with  the 
AFSCN  ground  stations 

3.  Act  as  controlling  authority  for  the  transmission  security  (TRANSEC) 
crypto  key  system,  and  direct  uploading  of  the  new  key  segment 
variables  to  the  onboard  TRANSEC  device  as  required 

Air  Force  Satellite 
Control  Network 
(AFSCN) 

1 .     Coordinate  with  NAVSOC  concerning  satellite  health/welfare  and 
anomaly  detection/resolution 

Naval  Computer  and 
Telecommunications  Area 
Master  Stations 
(NCTAMS)  and  Joint 
Fleet  Telecommunication 
Operations  Centers 
(JFTOCs) 

1.  Coordinate  with  NAVSOC  concerning  satellite  payload  control 

2.  Provide  connectivity  to  the  Defense  Switching  Network  (DSN)  and 
Public  Switching  Telephone  Network  (PSTN) 

Users 

1 .  Obtain  approval  of  their  requirements  through  the  ICDB  process 

2.  Adhere  to  the  network/terminal  configurations  and  precedences  assigned 
by  the  CINC  communications  manager 

3.  Ensure  that  payload  resources  are  used  as  assigned 

Table  1-1:  Responsibilities  of  Organizations  Involved  with  MUSCRAT 


3.  Concept  of  Operations 

MUSCRAT  will  be  launched  out  of  Cape  Canaveral,  Florida  onboard  either  a 
Delta  III  or  one  of  the  Atlas  family  of  launch  vehicles  (HAS,  IIAR,  IIARS).  Operational 
checkout  is  expected  to  last  for  approximately  one  month,  after  which  time,  the  satellite 
will  be  placed  in  a  non-operational  mode  until  the  entire  constellation  has  been  placed  in 
orbit.  The  MUSCRAT  system  will  then  be  activated. 


Figure  1-2:  MUSCRAT  Concept  of  Operations 


Figure  1-2  illustrates  MUSCRAT's  concept  of  operations.  The  system  uses 
current  Air  Force  Satellite  Control  Network  (AFSCN)  stations  for  controlling  satellite  bus 
operations,  as  well  as  for  monitoring  standard  satellite  telemetry  parameters.  Control  of 
payload  operations  will  be  accomplished  through  the  current  Naval  Computer  and 
Telecommunications  Area  Master  Stations  (NCTAMS)  and  Joint  Fleet 
Telecommunications  Operations  Centers  (JFTOCs)  via  an  S/EHF  link.    The  MUSCRAT 


design  allows  for  connectivity  of  payload  control  commands  when  a  satellite  is  not  within 
line  of  sight  of  a  NCTAMS/JFTOC  via  EHF  satellite  crosslinks  (60/61  GHz). 

MUSCRAT  provides  direct  UHF  connectivity  between  users  (uplink:  290- 
320  MHz,  downlink:  243-270  MHz).  If  one  of  the  users  are  non-mobile  and  is  using  the 
Defense  Switching  Network  (DSN)  or  the  Public  Switching  Telephone  Network  (PSTN), 
gateways  located  at  the  NCTAMS/JFTOCs  will  be  used.  This  greatly  increases 
MUSCRAT's  capacity  and  flexibility.  The  gateways  connect  to  the  satellite  via  an  S/EHF 
link  (downlink:  20  GHz,  uplink:  44  GHz).  The  EHF  satellite  crosslinks  will  also  be  used 
for  payload  communications  connectivity.  Examples  of  these  different  modes  are  shown 
in  the  figure. 


4.  Spacecraft  Description 

MUSCRAT  will  provide  an  assured,  netted,  global  (including  the  North  and  South 
Poles)  UHF  communications  link  between  both  mobile  and  fixed-site  users.  The 
MUSCRAT  constellation  consists  of  ten  operational  satellites  divided  into  two  orbital 
planes  separated  by  180°  in  longitude  of  ascending  node.  Their  circular  orbits  have  an 
altitude  of  1 1,000  km  and  an  inclination  of  48  degrees.  MUSCRAT  can  be  launched  on 
Atlas  IIAR,  Atlas  HAS,  and  Delta  III  launch  vehicles.  The  constellation  guarantees  global 
coverage  with  at  least  one  satellite  within  view  of  a  user  with  a  ten  degree  elevation  angle. 
[Ref.  40] 

Figures  1-3  through  1-5  depict  MUSCRAT's  deployed  and  stowed  configurations. 
MUSCRAT  is  basically  rectangular  in  shape,  with  a  bus  height  of  six  meters  and  a  base 
measuring  two  meters  square.  The  panels  for  the  main  "box"  are  constructed  of  aluminum 
honeycomb  and  are  supported  with  aluminum  shear  webs.  MUSCRAT  also  uses  an 
aluminum  central  cylinder  for  load  support.  [Ref.  40] 


Figure  1-3:  MUSCRAT  Deployed  Configuration 
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Figure  1-4:  MUSCRAT  Deployed  Configuration 
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Figure  1-5:  MUSCRAT  Stowed  Configuration  Inside  Delta  Fairing 
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The  bus  has  a  modular  design  (Figure  1-6).  The  upper  module  contains  the  UHF 
components.  The  middle  module  houses  the  EHF  downlink  components,  solar  arrays, 
attitude  control  equipment,  and  control  electronics.  The  lower  module  includes  the  EHF 
satellite  crosslink  components  and  batteries.  The  aluminum  central  cylinder  runs  the 
length  of  the  bus,  and  encloses  one  pressurant  and  two  propellant  tanks  (Figure  1-7). 
[Ref.  40] 


UPPER  MODULE: 

UHF  Payload  Components 
Heat  Pipes 


MIDDLE  MODULE: 
Central  Computer 
EHF  Gateway  Components 
ADCS  Components 
EPS  Solar  Arrays  and 
Control  Electronics 


LOWER  MODULE: 

EHF  Satellite  Crosslink  Components 
EPS  Batteries 


Figure  1-6:  MUSCRAT's  Modular  Bus  Design 

The  most  prominent  feature  of  MUSCRAT  is  its  8.5  meter  diameter  reflector.  It 
deploys  on  a  graphite  truss  structure  5.3  meters  above  the  bus.  The  ground  plane,  with  a 
diameter  of  four  meters,  contains  19  crossed  dipole  antennas.  Twenty-six  high  heat 
dissipating  solid  state  power  amplifiers  (SSPA's)  are  mounted  on  the  north  and  south 
panels  of  the  spacecraft.  Seven  of  these  are  for  redundancy.  The  remaining  UHF 
components,  except  the  data  processing  units,  are  located  throughout  the  upper  module. 
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Figure  1-7:  MUSCRAT's  Central  Cylinder 

MUSCRAT  uses  an  EHF  link  to  connect  to  gateway  ground  stations.  These 
stations  are  used  to  link  communications  to  the  Defense  Switching  Network  (DSN)  and 
the  Public  Switching  Telephone  Network  (PSTN).  Two  EHF  antennas  (plus  one  for 
redundancy)  accomplish  this  connection.  MUSCRAT  also  uses  EHF  satellite  crosslinks  to 
ensure  timely  netted  communications.  Four  EHF  0.45  meter  dish  antennas  are  located  on 
the  lower  corners  of  the  bus. 

Control  of  spacecraft  operations,  with  the  exception  of  attitude  determination  and 
control,  is  accomplished  using  a  central  data  processor.  It  contains  eight  Motorola  Power 
PC  603e  microprocessors  mounted  on  modified  Radstone  PCP2  conduction  cooled 
integrated  circuit  cards. 
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MUSCRAT  is  three-axis  stabilized  and  yaw  steered.  Pointing  requirements  of  0.1° 
are  met  using  three  reaction  wheels,  with  a  fourth  reaction  wheel  for  redundancy. 
Positional  and  three  axis  attitude  data  are  acquired  using  GPS  and  star  trackers.  The 
reaction  wheels  are  desaturated  using  magnetic  torque  rods  and  0.5  lb  thrust ers  as  primary 
and  secondary  systems  respectively.  The  attitude  determination  and  control  subsystem 
(ADCS)  is  controlled  by  its  own  redundant  processor.  Other  supporting  ADCS 
components  include  a  magnetometer,  two  inertial  reference  units  (IRU's),  sun  trackers 
located  on  the  solar  arrays,  and  two  earth/sun  sensor  assemblies  used  to  initially  orient  the 
spacecraft  once  it  gets  to  orbit. 

The  telemetry,  tracking,  and  control  (TT&C)  subsystem  provides  for  both 
autonomous  operations  and  direct  ground  control  of  MUSCRAT.  The  satellite  has  two 
patch  TT&C  omni-directional  antennas  for  assured  telemetry  download  and  command 
reception.  Onboard  TT&C  operations  will  be  controlled  by  redundant  spacecraft  control 
electronics.  TT&C  ground  functions  will  take  place  at  Air  Force  Satellite  Control 
Network  (AFSCN)  Space-Ground  Link  System  (SGLS)  stations. 

The  electric  power  subsystem  (EPS)  is  a  dual  bus,  with  a  fully  regulated  bus 
voltage  of  5 1  volts.  MUSCRAT  requires  8083  W  of  power  in  non-eclipse  periods.  Two 
GaAs  solar  arrays  with  a  combined  area  of  30  m2  provide  this  power,  and  it  is  stored  in 
four  NiH2  batteries,  each  having  27  cells.  Shunt  regulators  are  located  on  the  solar  array 
yokes. 

The  propulsion  subsystem  is  comprised  of  two  monopropellant  (hydrazine)  tanks 
and  a  nitrogen  pressurant  tank  housed  inside  MUSCRAT' s  central  cylinder.  Propellant 
will  be  used  for  satellite  re-phasing,  orbital  insertion  corrections,  as  a  backup  for  reaction 
wheel  desaturation,  and  for  satellite  deorbit  at  the  end  of  mission  life.  Twelve  0.5  lb 
thrusters  are  used  for  reaction  wheel  desaturation,  and  two  5.0  lb  thrusters  are  used  for 
gross  orbit  maneuvers.  All  thruster  valves  have  dual  series  seat  valves  for  redundancy 

MUSCRAT' s  thermal  subsystem  provides  the  required  temperature  control 
throughout  the  satellite's  lifetime.  The  primarily  passive  system  design  uses  heat  pipes  in 
the  north  and  south  panels  to  cool  concentrated  heat  sources  (i.e.  SSPA's)  and  minimize 
panel  temperature  gradients.    These  panels  are  radiators,  and  are  covered  with  OSR. 
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Multi-layer  insulation  (MLI)  blankets  are  also  used  to  maintain  proper  temperatures. 
Heaters  are  used  to  maintain  the  batteries  and  hydrazine  system  within  required 
temperature  limits. 

Table  1-2  is  a  summary  of  MUSCRAT's  mass  budget.  The  total  launch  mass  of 
2380.5  kg  is  broken  down  by  subsystem  and  includes  a  15%  mass  margin,  propellant,  and 
the  spacecraft  adapter. 


Subsystem 

Mass  (kg) 

%  Dry  Mass 

Structure 

214.6 

12.1  % 

Payload/TT&C 

871.9 

49.4  % 

EPS 

336.6 

19.0% 

Thermal 

50 

2.8  % 

ADCS 

90.7 

5.1  % 

Propulsion 

45.2 

2.6  % 

Integration 

110.1 

6.2  % 

Mechanisms 

50 

2.8  % 

Dry  Mass 

1769.1 

— 

1 5%  Margin 

265.3 

— 

Propellant 

246.1 

— 

Spacecraft  Adapter 

100 

___ 

Launch  Mass 

2380.5 

— 

Table  1-2:  MUSCRAT  Mass  Summary 

Table  1-3  summarizes  MUSCRAT's  power  budget  during  eclipse  and  non-eclipse 
periods.  It  includes  a  5%  growth  margin  in  sun,  a  10%  growth  margin  in  eclipse,  and  an 
additional  10%  solar  array  design  margin. 
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Subsystem 

Non-Eclipse  Requirements  (W) 

Eclipse  Requirements  (W) 

Payload 

5661 

5661 

EPS 

1102 

105 

ADCS 

205.7 

205.7 

Thermal 

30 

30 

Propulsion  * 

108 

108 

Margin 

1084.7 

600.2 

Total 

8083.4 

6601.9 

*  Not  included  in  total  power  requirements  due  to  extremely  low  average 
power  requirements 

Table  1-3:  MUSCRAT's  Power  Budget 

Table  1-4  lists  MUSCRAT's  propellant  budget.     The  total  propellant  mass  of 
246.1  kg  includes  a  10%  margin. 


Propellant  Use 

Mass  (kg) 

Orbit  Injection  Correction 

12.69 

Orbit  Rephasing 

72.54 

Reaction  Wheel  Desaturation 

17.28 

Deorbit 

118.98 

10%  Margin 

24.61 

Total  Propellant  Mass 

246.1 

Table  1-4:  MUSCRAT's  Propellant  Budget 
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B.  ATTITUDE  DYNAMICS  AND  CONTROL  SUBSYSTEM  - 
CONCEPTS  AND  DESIGN  PROCESS 

The  attitude  determination  and  control  subsystem  (ADCS)  stabilizes  the  vehicle 
and  orients  it  in  desired  directions  during  the  mission  despite  the  external  disturbance 
torques  acting  on  it.  This  requires  that  the  vehicle  determine  its  attitude,  using  sensors, 
and  control  it,  using  actuators.  The  ADCS  is  often  tightly  coupled  to  other  subsystems  on 
board,  especially  the  propulsion  and  navigation  functions.  [Ref.  6:  p.  340] 

This  paper  begins  by  discussing  several  useful  concepts  and  definitions,  including 
mass  properties,  disturbance  torques,  angular  momentum,  and  reference  vectors.  The 
mass  properties  of  a  spacecraft  are  key  in  determining  the  size  of  control  and  disturbance 
torques.  And,  a  knowledge  of  the  location  of  the  center  of  mass  or  gravity  (eg)  as  well  as 
the  elements  of  the  inertia  matrix  is  a  requirement.  As  well  as  the  moments  and  products 
of  inertia  about  chosen  reference  axes.  The  direction  of  the  principal  axes,  those  axes  for 
which  the  inertia  matrix  is  diagonal  and  the  products  of  inertia  are  zero,  are  also  of 
interest.  Finally,  we  need  to  know  how  these  properties  change  with  time,  as  fuel  or  other 
consumables  are  used,  or  as  appendages  are  moved  or  deployed.  [Ref.  6:  p.  340] 

A  body  in  space  is  subject  to  small  but  persistent  disturbance  torques  (i.e., 
10"4  Nm)  from  a  variety  of  sources.  These  torques  are  categorized  as  cyclic,  varying  in  a 
sinusoidal  manner  during  an  orbit,  or  secular,  accumulating  with  time,  and  not  averaging 
out  over  an  orbit.  These  torques  would  quickly  reorient  the  vehicle  unless  resisted  in 
some  way.  An  ADCS  system  resists  these  torques  either  passively,  by  exploiting  inherent 
inertia  or  magnetic  properties  to  make  the  "disturbances"  stabilizing  and  their  effects 
tolerable,  or  actively,  by  sensing  the  resulting  motion  and  applying  corrective  torques. 
[Ref.  6:  p.  340] 

Angular  momentum  plays  an  important  role  in  space,  where  torques  typically  are 
small  and  spacecraft  are  unconstrained.  For  a  body  initially  at  rest,  an  external  torque  will 
cause  the  body  to  angularly  accelerate  proportionally  to  the  torque  -  resulting  in  an 
increasing  angular  velocity.  Conversely,  if  the  body  is  initially  spinning  about  an  axis 
perpendicular  to  the  applied  torque,  then  the  body  spin  axis  will  precess,  moving  with  a 
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constant  angular  velocity  proportional  to  the  torque.  Thus,  spinning  bodies  act  like 
gyroscopes,  inherently  resisting  disturbance  torques  in  two  axes  by  responding  with 
constant,  rather  than  increasing,  angular  velocity.  This  property  of  spinning  bodies,  called 
gyroscopic  stiflhess,  can  be  used  to  reduce  the  effect  of  small,  cyclic  disturbance  torques. 
This  is  true  whether  the  entire  body  spins  or  just  a  portion  of  it,  such  as  a  momentum 
wheel  or  spinning  rotor.  [Ref.  6:  p.  340] 

Conservation  of  vehicle  angular  momentum  requires  that  external  torques  change 
the  system  angular  momentum.  Thus,  external  disturbances  must  be  corrected  by  external 
control  torques  (i.e.,  thrusters  or  magnetic  torquers)  or  the  resulting  momentum  buildup 
must  be  stored  internally  (i.e.,  by  reaction  wheels)  without  reorienting  the  vehicle  beyond 
its  allowable  limits.  The  momentum  build  up  due  to  secular  disturbances  ultimately  must 
be  reduced  by  applying  compensating  external  control  torques.  [Ref.  6:  p.  340] 

Often,  in  addition  to  rejecting  disturbances,  the  ADCS  must  reorient  the  vehicle  (in 
slew  maneuvers)  to  re-point  the  payload,  solar  arrays,  or  antennas.  These  periodic  re- 
pointing  requirements  may  drive  the  design  to  larger  actuators  than  would  be  required  for 
disturbance  rejection  alone.  [Ref.  6:  p.  341] 

To  orient  the  vehicle  correctly,  external  references  must  be  used  to  determine  the 
vehicle's  absolute  attitude.  These  references  include  the  Sun,  the  Earth's  IR  horizon,  the 
local  magnetic  field  direction,  and  the  stars.  In  addition,  inertial  sensors  (gyroscopes)  also 
can  be  carried  to  provide  a  short-term  attitude  reference  between  external  updates. 
External  references  (i.e.,  Sun  angles)  are  usually  measured  as  body-centered  angular 
distances  to  a  vector.  Such  measurements  provide  only  two  of  the  three  independent 
parameters  needed  to  specify  the  orientation  of  the  spacecraft.  This  results  in  the  need  for 
multiple  sensor  types  on  board  most  spacecraft.  [Ref.  6:  p.  341] 

The  Attitude  Dynamics  and  Control  Subsystem  design  process  is  summarized  very 
nicely  with  Table  1-5.  Much  of  the  input  information  for  this  table  comes  from  other 
subsystems  and  the  overall  concept  of  operations  for  the  satellite.  This  can  be  seen  in 
Figure  1-8.  [Ref.  6:  p.  344] 
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Step 

# 

Description 

Inputs 

Outputs 

la) 

-  Define  control  modes. 

-  Mission 

-  List  of  different  control 

lb) 

-  Define  or  derive 

requirements,  mission 

modes  during  mission  (table 

system  level 

profile,  type  of 

7-2). 

requirements. 

insertion  for  launch 

-  Requirements  and 

vehicle 

constraints  (table  7-3). 

2 

-  Select  type  of 

-  Payload,  thermal  & 

-  Method  for  stabilizing  and 

spacecraft  control  by 

power  needs 

control:  three  axis,  spinning, 

attitude  control  mode 

-  Orbit,  pointing 

or  gravity  gradient 

and  subsystem 

direction 

functional  description 

-  Disturbance 
Environment 

3 

-  Quantify  disturbance 

-  Spacecraft 

-  Values  for  forces  from 

environment 

geometry,  orbit, 

gravity  gradient,  magnetic 

solar/magnetic 

aerodynamics,  solar 

models,  mission 

pressure,  internal 

profile 

disturbances,  and  powered 
flight  effects  on  control  (eg 
off-sets,  slosh) 

4 

-  Select  &  size  ADCS 

-  Spacecraft 

-  Sensor  suite:  Earth,  Sun, 

hardware 

geometry,  pointing 

inertial,  or  other  sensing 

accuracy,  orbit 

devices 

conditions,  mission 

-  Control  actuators,  i.e., 

requirements,  lifetime, 

reaction  wheels,  thrusters, 

orbit,  pointing 

or  magnetic  torquers 

direction,  slew  rates 

-  Data  processing 
electronics,  if  any,  or 
processing  requirements  for 
other  subsystems  or  ground 
computer 

5 

-  Define  determination 

-  All  of  above 

-  Algorithms,  parameters, 

and  control  algorithms 

and  logic  for  each 
determination  and  control 
mode 

6 

-  All  of  above 

-  Refined 
requirements  and 
design 

-  Subsystem 
specifications 

Table  1-5:  Attitude  Dynamics  Control  Sub-System  Design  Process 
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Figure  1-8:  The  Impact  of  Mission  Requirements  and  Other  Subsystems  on  the  ADCS 
Subsystem.  Direction  of  arrows  show  requirements  flow  from  one  subsystem  to  another. 
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The  first  step  in  designing  the  Attitude  Dynamics  and  Control  Subsystem  is  to 
define  the  control  modes.  Each  control  mode  is  related  to  the  different  phases  of  the 
spacecraft  life.  The  ADCS  for  this  satellite  will  have  eight  control  modes  which  are 
described  in  chapter  2.  The  sensor  and  actuator  requirements  for  each  control  mode  may 
vary  drastically.  As  a  result,  the  designer  of  the  ADCS  must  ensure  the  requirements  of 
each  control  mode  are  satisfied  while  at  the  same  time  producing  as  simple  and  reliable  a 
subsystem  as  possible. 

Next,  a  thorough  understanding  of  the  system  level  requirements  is  required. 
These  requirements  usually  originate  from  the  other  spacecraft  subsystems  and  flow  down 
to  the  ADCS  designer.  For  example,  the  payload  designer  in  this  design  specified  the 
UHF  antenna  must  maintain  a  pointing  accuracy  of  0.1°.  Then,  the  systems  engineer  and 
structures  designer  determined  the  UHF  antenna  would  be  fixed  to  the  spacecraft.  This 
can  then  be  translated  into  a  spacecraft  attitude  control  accuracy  of  0.1°  along  the  roll  and 
pitch  axes. 

Then,  based  on  mission  requirements  and  inputs  from  other  subsystems,  a  decision 
has  to  be  made  on  the  type  of  spacecraft  control  to  be  used.  For  example,  a  spin  stabilized 
spacecraft  may  be  most  suitable  for  some  missions  while  a  three-axis  control  momentum 
wheel  system  may  be  best  for  others.  When  the  type  of  spacecraft  control  is  determined 
an  understanding  of  the  type  of  sensors  and  actuators  should  be  developed.  This  is  not 
actually  selecting  the  components,  but  instead  performing  trade  studies  for  determining 
which  types  of  sensor  and  actuator  are  best  suited  for  the  spacecraft  mission.  An  example 
might  be  determining  whether  to  use  a  star  sensor  or  earth  sensor.  Chapter  three 
discusses  these  issues  and  presents  a  functional  description  of  the  attitude  dynamics  and 
control  subsystem. 

The  next  step  in  the  design  process  is  to  quantify  the  disturbance  environment. 
There  are  generally  four  external  disturbances  to  consider:  Gravity  Gradient,  Solar 
Pressure,  Aerodynamic  Pressure  and  Geomagnetic  field  disturbances.  These  disturbances 
should  be  estimated  so  that  appropriate  actuators  and  sensors  can  be  selected  to 
counteract  their  effects.  Chapter  four  presents  a  thorough  explanation  on  quantifying  the 
external  disturbance  torques. 
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With  information  from  the  analysis  discussed  above  a  selection  of  ADCS 
components  was  performed.  The  components  of  the  ADCS  include  sensors,  actuators, 
and  control  processors.  Chapter  5  presents  this  selection  process  and  provides  a  detailed 
description  of  each  component  selected. 

The  final  portion  of  the  design  process  is  the  definition  of  determination  and 
control  algorithms.  This  also  includes  simulation  and  modeling  of  the  ADCS.  This  begins 
in  chapter  six  which  focuses  exclusively  on  yaw  controlling  and  control  algorithms. 
Chapter  7  then  focuses  on  simulating  the  motion  of  the  spacecraft  and  presents  a  thorough 
explanation  of  the  differences  between  a  non  yaw  steering  satellite  and  a  yaw  steering  one. 
The  thesis  concludes  with  several  spacecraft  attitude  control  simulations. 


21 


This  Page  Intentionally  Left  Blank 


22 


II.  DEFINE  CONTROL  MODES  AND  SYSTEM  LEVEL 

REQUIREMENTS 

Attitude  control  modes  are  related  to  the  different  phases  of  the  spacecraft  life. 
Basically,  the  ADCS  must  be  designed  to  meet  the  requirements  of  each  mode.  These 
ADCS  requirements  are  tightly  coupled  to  other  spacecraft  subsystems  and  the  overall 
mission  needs.  Since  the  requirements  change  drastically  with  each  mission  phase,  the 
ADCS  designer  may  be  challenged  to  built  one  system  that  meets  all  requirements. 

A.  DEFINE  CONTROL  MODES 

1.  Orbit  Insertion  Mode 

The  orbit  insertion  is  the  period  during  and  after  boost  while  the  spacecraft  is 
brought  to  final  orbit.  Options  in  this  control  mode  include  no  spacecraft  control,  simple 
spin  stabilization  of  solid  rocket  motor,  and  full  spacecraft  control  using  liquid  propulsion 
system. 

2.  Acquisition  Mode 

This  is  normally  where  the  initial  determination  of  attitude  and  stabilization  of  the 
vehicle  occurs.  But,  attitude  determination  and  stabilization  may  also  need  to  occur  due 
to  recovery  from  power  upsets  and  /  or  emergencies. 

3.  Sun  Acquisition  Mode  /  Sun  Hold  Mode 

After  attitude  determination  and  stabilization  is  complete  Sun  acquisition  takes 
place  in  order  to  ensure  an  adequate  power  supply. 

4.  Earth  Acquisition  Mode 

A  critical  aspect  of  this  mode  is  acquiring  a  nadir  pointing  payload  orientation 
while  simultaneously  keeping  the  solar  arrays  tracking  the  sun.  Information  from  the  sun 
sensors  and  additional  navigational  aids  allow  the  spacecraft  to  re-orient  such  that 
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simultaneous  Nadir  payload  pointing  and  solar  array  sun  tracking  is  attained.   Much  more 
will  be  discussed  on  this  topic  later. 

5.  Normal,  On-Station  Mode 

The  Normal,  On-Station  mode  is  used  for  the  vast  majority  of  the  mission.  As  a 
result,  the  requirements  for  this  mode  will  drive  the  system  design.  And,  the  majority  of 
discussion  will  focus  on  this  mode  of  operation. 

6.  Slew  Mode 

This  mode  occurs  when  a  reorientation  of  the  vehicle  is  required  (i.e.  the 
spacecraft  executes  a  slew  maneuver). 

7.  Contingency  /  Safe  Mode 

This  mode  of  operation  is  used  in  emergencies  for  when  the  regular  mode  fails  or 
is  disabled.  Also,  this  mode  may  be  selected  to  meet  unexpected  requirements  like  power 
or  thermal  constraints;  at  the  expense  of  normal  operations. 

8.  Orbit  Maintenance  /  Orbit  Transition  Mode 

This  mode  occurs  when  the  spacecraft  needs  to  be  repositioned  in  the  orbit.  This 
either  could  be  due  to  an  accumulation  of  error  over  several  years;  or  due  to  a  satellite 
failure  which  requires  selected  satellites  to  be  repositioned  in  the  orbit. 

B.  SYSTEM  LEVEL  REQUIREMENTS 

Following  are  the  system  level  requirements  for  the  Attitude  Dynamics  and 
Control  Subsystem  of  the  MUS  communications  satellite. 

1 .  The  ADCS  must  be  fully  operational  without  the  aid  of  ground  base  control. 

2.  The  operational  lifetime  of  the  ADCS  must  meet  or  exceed  ten  years. 

3.  The  ADCS  must  maintain  a  UHF  antenna  beam  pointing  accuracy  of  less  than 
or  equal  to  0.1°. 
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4.  The  ADCS  must  have  the  capability  to  provide  continuous  yaw  steering  and 
solar  array  tracking  such  that  the  solar  array  pointing  accuracy  is  better  than  or  equal  to 
5°. 

5.  The  spacecraft  must  be  fully  operational  during  eclipse. 

6.  The  spacecraft  must  be  stabilized  and  under  control  during  all  phases  of  the 
mission. 

7.  The  ADCS  system  should  have  system  redundancy  such  that  there  are  no  single 
point  failures. 

8.  The  ADCS  should  have  a  ground  override  capability. 

9.  The  spacecraft  must  have  a  cross-link  antenna  pointing  accuracy  of  0.07°  which 
requires  a  spacecraft  attitude  knowledge  requirement  of  0.02°  in  roll,  pitch,  and  yaw. 

10.  The  Attitude  Dynamics  and  Control  Subsystem  must  have  a  positional 
knowledge  requirement  of  350  meters. 
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III.  FUNCTIONAL  DESCRIPTION  AND  TYPE  OF  SPACECRAFT 
CONTROL  BY  ATTITUDE  CONTROL  MODE 

A.       ADCS    FUNCTIONAL    DESCRIPTION    AND    COMPONENT 
OVERVIEW 

A  Zero  Momentum  three  axis  system  was  selected  which  can  provide  any 
combination  of  slew  maneuvers  in  the  pitch,  roll  and  yaw  directions.  The  Attitude 
Dynamics  and  Control  Subsystem  will  maintain  Nadir  pointing  accuracies  to  within  0.1°. 
At  the  same  time  it  will  continuously  rotate  in  the  yaw  direction  while  the  solar  arrays  are 
simultaneously  rotating  in  order  to  keep  a  perpendicular  orientation  towards  the  sun.  This 
will  greatly  decrease  the  complexity  of  the  solar  array  control  mechanisms.  In  a  zero- 
momentum  system,  reaction  wheels  respond  to  disturbances  on  the  vehicle.  For  example, 
a  vehicle-pointing  error  creates  a  signal  which  speeds  up  the  wheel,  initially  at  zero.  This 
torque  corrects  the  vehicle  and  leaves  the  wheel  spinning  at  low  speed,  until  another 
pointing  error  speeds  the  wheel  further  or  slows  it  down  again.  If  the  disturbance  is  cyclic 
during  each  orbit,  the  wheel  may  not  approach  saturation  speed  for  several  orbits.  Secular 
disturbances,  however,  cause  the  wheel  to  drift  toward  saturation.  We  then  must  apply  an 
external  torque  with  a  thruster,  a  magnetic  torquer  or  both,  to  force  the  wheel  speed  back 
to  zero.  This  process  is  called  desaturation,  momentum  unloading,  or  momentum 
dumping  and  can  be  done  automatically  or  by  command  from  the  ground.  [Ref.  6:  p.  348] 

The  Zero  Momentum  Three  Axis  Control  system  chosen  was  one  of  three  options 
seriously  considered.  The  other  two  were  Dual  Spin  Stabilization  and  a  Zero  Momentum 
Control  Momentum  Gyro.  The  Dual  Spin  Stabilization  was  dismissed  because  the  antenna 
configuration  did  not  fit  well  on  the  despun  platform  and  the  power  requirements  were 
greater  than  a  Dual  Spin  Stabilization  could  provide.  The  Control  Moment  Gyro  (CMG) 
would  work  for  this  system,  but  it  was  dismissed  because  it  is  heavier  and  would  not  be 
conducive  to  yaw  steering.  Yaw  steering  enables  a  decrease  in  mass  and  complexity  of  the 
solar  array  drive  assembly.  In  addition  the  Zero  Momentum  Three  Axis  Stabilization 
system  would  be  lighter.    Also,  the  Control  Moment  Gyro  is  best  suited  for  high  torque 
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demands  and  fast  slewing  requirements.  Neither  of  which  we  have.  It  should  be  pointed 
out  however,  that  the  three  axis  stabilization  has  more  moving  parts  and  could  therefore 
be  considered  more  complex.  So,  a  legitimate  argument  could  be  made  for  either  system. 
Figure  3-1  is  a  general  block  diagram  of  the  Zero  Momentum  Attitude  Dynamics  and 
Control  Subsystem.  [Ref.  6:  p.  348,  Ref.  21:  p.  46] 
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Figure  3-1 :  Block  Diagram  of  Attitude  Determination  and  Control  Subsystem 
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1.  TT&C 

Telemetry,  Tracking  and  Control  will  take  place  at  the  seven  Air  Force  Satellite 
Control  Network  ground  stations  with  one  of  the  stations  exercising  overall  control. 
Signals  will  be  sent  to  the  spacecraft  as  required  where  they  will  be  processed  in  the 
Spacecraft's  main  computer.  Those  commands  relating  to  Attitude  Dynamics  and  Control 
will  be  sent  to  the  ADCS  central  processor  and  then  executed.  In  addition,  data  will  be 
sent  from  the  ADCS  system  to  the  ground  stations  in  order  to  monitor  the  health,  welfare 
and  position  of  both  the  satellite  and  the  ADCS.  The  TT&C  subsystem  will  consist  of 
dedicated  antennas,  a  receiver,  a  transmitter,  and  a  command  and  data  handling  system. 
[Ref.  15:  p.  13] 

2.  Reaction  Wheels 

The  Zero  Momentum  reaction  wheel  system  will  counteract  disturbance  torque 
and  provide  the  momentum  to  perform  slew  maneuvers  as  required.  Reaction  wheels  are 
characterized  by  being  light  weight  and  capable  of  producing  a  continuous  range  of 
torque.  This  will  work  well  for  our  satellite  since  we  will  continuously  need  to  track  the 
sun  in  order  to  attain  maximum  capacity  with  our  solar  arrays.  [Ref.  21 :  p.  47] 

3.  Magnetic  Torque  Rods 

Magnetic  Torque  Rods  will  be  used  to  efficiently  dump  angular  momentum  in  the 
reaction  wheels  produced  by  secular  disturbance  torques.  By  providing  electric  current 
through  the  torque  rod  coil  a  magnetic  dipole  moment  results.  This  moment  causes  a 
torque  that  tries  to  align  the  torque  rod  with  the  earth's  magnetic  field.  When  this  torque 
is  produced  an  equal  and  opposite  torque  will  be  produced  by  the  Reaction  Wheels;  thus 
allowing  a  means  to  dump  momentum.  [Ref.  21 :  p. 5 5] 

4.  Star  Sensor 

Based  on  the  crosslink  antenna  pointing  requirements  of  0.02°  some  type  of  star 
sensing  system  will  have  to  be  used.    In  order  to  provide  a  continuous  star  sensing  three 
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star  sensors  will  be  required.  In  the  worst  case  scenario  two  of  the  three  sensors  may  be 
obstructed  by  the  sun  and  moon.  If  this  occurs  the  Inertial  Reference  Unit  will  be  able  to 
provide  the  required  data  for  the  short  period  of  time  required. 

5.  Global  Positioning  System 

A  Global  Positioning  System  (GPS)  is  required  on  the  spacecraft  in  order  to  allow 
for  an  ability  to  accurately  determine  orbit  position  without  the  aid  of  a  ground  station. 
And,  the  GPS  will  be  a  back-up  for  attitude  determination.  [Ref.  21 :  p.  53] 

6.  Inertial  Reference  Unit 

The  Inertial  Reference  Unit  will  provide  continuous  angular  velocity  information. 
Then,  with  the  knowledge  of  angular  velocity,  angular  acceleration  can  be  obtained  using 
a  pseudorate  modulator.  In  the  event  of  star  tracker  failure  or  temporary  un-availability, 
the  Inertial  Reference  Unit  will  also  provide  continuous  attitude  information  for  a  limited 
period  of  time.  [Ref.  21:  p.  52] 

7.  Three  Axis  Magnetometer 

Magnetometers  are  simple,  reliable,  lightweight  sensors  that  measure  both  the 
direction  and  size  of  the  Earth's  magnetic  field.  When  compared  to  the  Earth's  known 
field,  their  output  helps  establish  the  spacecraft's  attitude.  But  their  accuracy  is  not  as 
good  as  that  of  star,  sun  or  horizon  references.  The  Earth's  field  can  shift  with  time  and 
cannot  be  calculated  precisely  in  the  first  place.  To  improve  accuracy,  data  is  often 
combined  with  the  data  from  Sun  or  horizon  sensors.  When  the  vehicle  is  using  the 
magnetic  torquers  and  passes  through  magnetic  field  reversals  during  each  orbit,  the 
magnetometer  will  be  used  to  control  the  polarity  of  the  torquer  output.  And,  the 
torquers  will  have  to  be  turned  off  while  the  magnetometer  is  sampled  to  avoid  corrupting 
the  measurement.  [Ref.  6:  p.  361] 
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8.  Central  Processor 

The  Attitude  Dynamics  and  Control  Subsystem  Central  Processor  will  be  a 
commercial  off  the  shelf  product  produced  by  Honeywell.  The  Central  Processor  is  the 
brain  of  the  ADCS.  As  can  be  seen  from  Figure  3-1,  input  data  is  fed  into  the  central 
processor  from  the  TT&C,  Star  Sensors,  GPS,  IRU,  Three  Axis  Magnetometer  and 
Reaction  Wheels.  Based  on  this  input  data  the  central  processor  will  calculate  the 
spacecraft's  position,  velocity,  orientation,  angular  velocity  and  angular  acceleration. 
Then,  the  central  processor  will  calculate  the  required  spacecraft  orientation,  angular 
velocity  and  angular  acceleration  profile  to  maintain  Nadir  earth  pointing,  an  adequate  sun 
pointing  profile,  momentum  limits  and  orbit  repositioning/transfers.  Finally,  appropriate 
commands  will  be  sent  to  the  reaction  wheels,  torque  rods,  thrusters  and  solar  array  drive 
mechanisms. 

9.  Sun  Sensor 

A  sun  sensor  is  a  visible  light  detector  which  can  measure  one  or  two  angles 
between  their  mounting  base  and  incident  sunlight.  One  sun  sensor  assembly  will  be 
mounted  on  or  near  each  solar  array  assembly  and  will  be  used  to  track  the  suns  rays. 

10.  Earth  Sensor 

An  Earth  Sensor  is  an  infrared  device  that  detects  the  contrast  between  the  cold  of 
deep  space  and  the  heat  of  the  Earth's  atmosphere.  These  sensors  provide  Earth  relative 
information  directly  for  earth  pointing  spacecraft,  which  simplifies  on  board  processing. 

B.  TYPE  OF  SPACECRAFT  CONTROL  BY  ATTITUDE  CONTROL 
MODE 

1.  Orbit  Insertion  Mode 

The  launch  vehicles  we  have  selected  will  allow  for  direct  orbit  insertion  with  no 
need  for  an  apogee  kick  motor  or  a  perigee  kick  motor.     As  a  result,  under  normal 
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condition  there  will  be  no  spacecraft  control  until  after  the  spacecraft  is  in  or  very  near  its 
final  orbit  and  has  separated  from  the  launch  vehicle. 

2.  Acquisition  Mode 

When  the  spacecraft  is  separated  from  the  launch  vehicle  the  acquisition  mode  of 
operation  will  begin.  At  this  time  the  spacecraft  will  normally  have  a  spin  about  the  z  axis. 
Each  of  the  two  sun  sensor  assemblies  have  two  perpendicular  fan  shaped  fields  of  view, 
ax  and  a2 .  Each  sensor  generates  an  output  pulse  as  the  sun  passes  through  the  field  of 
view  of  the  light  sensitive  element.  Then,  the  time  of  arrival  of  the  pulse  is  detected  in  the 
ADCS  central  processor.  In  addition,  the  earth  sensor  provides  a  pulse  which  is  also 
detected  by  the  central  processor.  The  central  processor  is  then  able  to  take  these 
measurements  and  compute  the  attitude  of  the  spacecraft  with  respect  to  the  sun  and 
earth. 

3.  Sun  Acquisition  Mode  /  Sun  Hold  Mode 

Using  angular  position,  angular  rate  and  data  from  the  pitch  sun  sensor,  the 
spacecraft  is  able  to  lock  into  and  hold  the  sun  in  pitch.  At  this  time  the  yaw  and  roll  axes 
are  under  gyro  reference  control.  Then,  a  slow  yaw  rotation  occurs  in  order  to  acquire  the 
sun  with  the  yaw  sun  sensor.  When  both  sensors  have  locked  onto  the  sun,  the  sun  hold 
mode  automatically  begins.  Then,  with  the  sun's  position  firmly  established,  the  solar 
arrays  are  deployed  and  the  remaining  ADCS  components  are  brought  on  line  and 
thoroughly  tested. 

4.  Earth  Acquisition  Mode 

The  earth  acquisition  mode  occurs  at  either  0600  or  1 800  local  time.  At  this  time 
the  earth  line  and  sun  line  relative  to  the  spacecraft  are  perpendicular.  With  the  yaw  and 
pitch  axis  locked  onto  the  sun  a  gentle  roll  maneuver  takes  place  until  the  earth  is  acquired 
by  the  earth  sensors.  When  an  adequate  Nadir  pointing  accuracy  is  attained  by  the  earth 
sensors,  both  roll  and  pitch  control  transfer  from  the  Gyros  over  to  the  earth  sensors. 
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5.  Normal,  On-Station  Mode 

The  Normal/On-Station  operating  mode  is  by  far  the  most  significant  of  the 
Attitude  Control  Modes.  The  vast  majority  of  the  spacecraft's  time  will  be  spent  under 
normal  operation.  As  a  result,  the  design  of  the  spacecraft's  Attitude  Dynamics  and 
Control  Subsystem  was  driven  by  this  mode.  And,  that  is  why  a  vast  majority  of  the 
discussion  will  be  on  the  Normal,  On-Station  operating  mode.  The  other  mode  operations 
will  largely  feed  from  how  the  Normal/On-Station  operations  are  conducted.  [Ref.  6:  p. 
340] 

6.  Slew  Mode 

We  will  utilize  the  Zero  Momentum  Three  Axis  Control  System  to  simplify  the 
solar  array  design.  The  spacecraft  will  continually  slew  about  the  unconstrained  yaw  axis 
in  order  to  keep  the  solar  array  axis  perpendicular  to  the  solar  radiation.  Thus,  the 
reduced  array  size  possible  with  two  degrees  of  freedom  can  be  achieved  with  one  array 
axis  drive  and  continual  spacecraft  rotation  about  the  yaw  axis.  It  is  important  to  ensure 
the  reaction  wheel  selected  can  deliver  the  required  torque.  Slewing  about  the  roll  and 
pitch  axis  will  not  normally  occur.  If  however,  commanded  to  do  so,  the  process  and 
dynamics  of  slewing  about  the  roll  and  pitch  axis  will  be  the  same  as  that  about  the  yaw 
axis.    [Ref.  6:  p.  351] 

7.  Contingency  /  Safe  Mode 

This  mode  is  enabled  in  emergencies  or  when  a  certain  level  of  uncertainty  exists 
concerning  the  operational  status  of  the  spacecraft.  In  the  Contingency/Safe  Mode  all 
non-essential  equipment  and  actuators  are  turned  off.  Then,  when  given  approval,  the 
spacecraft  will  enter  the  Sun  Acquisition  Mode  followed  by  extensive  testing  as  the 
various  components  and  actuators  are  brought  back  on  line. 

8.  Orbit  Maintenance  /  Orbit  Transition  Mode 

If  the  spacecraft  needs  to  be  repositioned  due  to  accumulated  orbital  error  or  a 
newly  assigned  orbit  position  it  will  enter  the  Orbit  Maintenance/Orbit  Transition  Mode. 
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The  first  step  in  this  mode  is  to  shut  down  non-essential  operations.  Then,  the  spacecraft 
will  orient  itself  such  that  the  yaw  axis  of  rotation  points  in  the  desired  direction  of 
movement.  Then  a  spin  up  will  take  place  followed  by  the  firing  of  appropriate  control 
thrusters.  When  the  spacecraft  has  reached  its  new  position  satisfactorily,  it  will  enter  into 
the  acquisition  mode. 

Table  3-1  summarizes  the  type  of  spacecraft  control  and  the  sensors  utilized  by 
each  attitude  control  mode. 
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MODE 

SENSORS 

ACTUATORS 

FUNCTION 

Orbit  Insertion 

-  None 

-  None 

Insert  spacecraft 
into  final  orbit 

Acquisition 

-  Sun  Sensors  (Pri) 

-  Torque  Rods  (Pri) 

Determine 

Mode 

-  Earth  Sensor  (Pri) 

-  Thrusters  (Sec) 

Spacecraft 

-  Sun  Sensors  (Sec) 

attitude  and 

-  Earth  Sensor  (Sec) 

stabilize  space- 
craft 

Sun 

-  Sun  Sensor  #1  (Pri) 

-  Reaction  Wheels  (Pri) 

Acquire  and  or 

Acquisition  / 

-  Sun  Sensor  #2  (Sec) 

-  Torque  Rods  (Sec) 

hold  a  sun 

Sun  Hold 

-  Thrusters  (Sec) 

pointing 

Mode 

orientation 

Earth 

-  Earth  Sun  Sensor 

-  Reaction  Wheels  (Pri) 

Acquire  and 

Acquisition 

Assembly  #1  (Pri) 

-  Torque  Rods  (Sec) 

maintain  Nadir 

Mode 

-Sun  Sensor  #1  (Pri) 

-  Earth  Sun  Sensor 
Assembly  #2  (Sec) 

-  Sun  Sensor  #2  (Sec) 

-  Thrusters  (Sec) 

Earth  pointing 

Normal,  On- 

-Sun  Sensor  #1  (Pri) 

-  Reaction  Wheel  (Pri) 

Maintain  Nadir 

Station 

-  Earth  Sun  Sensor  #1 

-  Torque  Rods  (Sec) 

Earth  pointing, 

(Pri) 

-  Thrusters  (Sec) 

Solar  Array 

-  Star  Tracker  (Pri) 

tracking  and 

-  Earth  Sun  Sensor  #2 

attitude  /  position 

(Sec) 

accuracy 

-  GPS  (Sec) 

requirements 

-  Sun  Sensor  #2  (Sec) 

Slew 

-  Star  Trackers  (Pri) 

-  Reaction  Wheels  (Pri) 

-  Reorient 

-  IRU  (Sec) 

-  Torque  Rods  (Sec) 

Spacecraft  as 

-  Thrusters  (Sec) 

required 

Contingency  / 

-Sun  Sensor  #1  (Pri) 

-  Reaction  Wheels  (Pri) 

Ensures  spacecraft 

Safe 

-  Earth/Sun  Sensor  #1 

-  Torque  Rods  (Sec) 

remains  within 

(Pri) 

-  Thrusters  (Sec) 

limited  operating 

-  Sun  Sensor  #2  (Sec) 

constraints 

-  Earth/Sun  Sensor  #2 

(Sec) 

Orbit 

-  GPS  (Pri) 

-  Thrusters  (Pri) 

spacecraft 

Maintenance  / 

-  IRU  (Sec) 

-  Thrusters  (Sec) 

reposition  / 

Orbit 

change  in  orbital 

Transition 

parameters 

Mode 

Table  3-1 :  Type  of  Spacecraft  Control  by  Attitude  Control  Mode 
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C.       THE    MOTION    OF    THE    SATELLITE    AND    THE    SOLAR 
ARRAYS 

As  a  spacecraft  decreases  in  altitude  the  time  spent  in  eclipse  becomes  a 
considerable  part  of  the  total  lifetime  of  the  satellite.  For  instance,  at  1 1 ,000  km  the 
maximum  eclipse  will  occur  when  the  Sun  is  in  the  orbital  plane  and  the  eclipse  time  will 
be  12%  of  the  orbital  period.  [Ref.  37:  pp.  75-78] 

To  take  maximum  advantage  of  the  solar  radiation,  the  solar  arrays  must  be  sun 
pointed.  At  the  same  time,  the  payload  requires  a  constant  nadir  pointing  configuration. 
At  every  position  in  the  orbit  these  two  demands  must  be  satisfied.  This  will  be 
accomplished  with  a  combination  of  yaw  steering  and  solar  array  pointing.  Assuming  the 
satellite  has  an  arbitrary  earth-pointing  attitude,  it  must  first  rotate  about  the  Z  s  -axis  (axis 
parallel  to  the  antenna;  yaw-axis)  with  an  angle  v|/,  so  that  the  axis  about  which  the  solar 
arrays  rotate  (BAPTA-axis)  is  in  the  plane  perpendicular  to  the  direction  of  the  solar 
radiation.  The  second  part  of  the  motion  consists  of  a  rotation  of  the  solar  arrays  about 
the  BAPTA-axis  with  an  angle  T|,  so  that  the  arrays  point  to  the  Sun.  These  two  motions 
will  now  be  looked  at  in  much  greater  detail.  [Ref.  37:  pp.  75-78] 

The  angles  r|  and  vj/  are  functions  of  time.  They  also  depend  on  the  position  of  the 
satellite  in  the  orbital  plane  and  the  position  of  the  Sun  with  respect  to  the  orbital  plane. 
To  aid  in  our  analysis  the  following  definitions  are  made:  First,  (3  is  defined  as  the  angle  in 
the  orbital  plane  between  the  point  with  minimum  distance  (sub-solar  point)  to  the  Sun 
and  the  satellite,  and  8  as  the  angle  between  the  orbital  plane  and  the  direction  of  the  Sun. 
This  is  illustrated  in  Figure  3-2.  Because  of  the  direct  relation  between  P  and  time,  \\f  and 
r\  will  be  expressed  as  functions  of  (3  and  5.  [Ref.  37:  pp.  75-78] 
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Figure  3-2:  Definition  of  the  angles  \\t,  r\,  P  and  5. 

The  derivation  of  these  functions  i|/(P,8)  and  r|(P,5)  will  be  done  below  by  means 
of  vector  algebra.  [Ref.  37:  pp.  75-78] 

The  direction  of  satellite  axes  X  s ,  Y  s ,  Zs  and  of  the  sunlight  s  in  the  inertial 
XYZ  coordinate  system  are: 
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Xs 


-sin(/?)*sin(y) 
cos(/?)*sin(y/) 
cos(y/) 


Ys  = 


sin(/?)*cos(y) 

~-cos(/?)~ 

cos(/?)*cos(y) 

;        Zs  = 

-sin(/?) 

-s\n(y/) 

0 

s  = 


-cos(S) 

0 
sin(^) 


While  the  vector  of  the  sunlight  s  stands  perpendicular  to  the  BAPTA-axis,  and 
because  the  BAPTA-axis  and  the  Ys-axis  are  parallel,  Ys   is  also  perpendicular  to  s. 

Now  \|/(P,5)  can  be  determined: 

Ys  •  s  =  -  sin(/?)  cos(^)  cos(5)  -  sin(y/)  sin(S)  =  0 


tan(^)  = 


-sin(y?) 


and 


y/  =  tan 


-/-sin(^) 


tan(J)  J  '  V  tan(£) 

The  direction  of  the  solar  panels  in  the  X  s ,  Y  s ,  Z  s  satellite  coordinate  system 


(see  figure  III. 2)  is:   ps  = 


-cos(rj) 
0 
Vsin(77)   J 


The  direction  of  the  solar  panels  in  the  XYZ  system  is  obtained  by  performing  the 
matrix  multiplication    p  =  \XS  Ys  Zsl*[/?5]: 


=>     P  = 


sin(/?)  *  sin(  y/  )*  cos(  if)  -  cos(/?)*  sin(  rf) 
- cos(/?)*  sin( y/ )* cos( rf)  -  sin(/?)*  sin( rj) 
-cos(^/)*cos(^) 


Since  p  and  s  have  to  point  exactly  in  opposite  direction,  also   rj(/3,5)  can  be 
determined: 
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p  x  s  =  det 


ex 

ey       ez 

P\ 

Pi       Pi 

S\ 

S2        h 

(-cos(/?)*sin(^)*cos(77)  -  sin(/?)*  sin(  rj))*  sm(S) 
(-sin(jfl)*sin(»*cos(77)  -  cos(/?)*sin(^))*sin(J)  +  cos(i//)*cos(tj)*cos(S) 
-(cos(/?)*sin(»*cos(/7)  +  sm(fi)*sm(?]))*cos(S) 


=  0 


tan(  tj)  = 


sin(^/) 
tan(/?) 


77  =  tan 


tan(/?)J 


^(/?)  ««^  ^(>0)  are  plotted  in  Figures  3-3  and  3-4.  From  these  figures  it  is  clear 
that  when  5  ->  0°  (when  the  Sun  is  almost  in  the  orbital  plane)  the  satellite  has  to  make 
radical  short-time  movements.  Especially  the  motion  of  the  satellite  about  the  yaw  axis 
nearby  p  =  0°  and  p  =  1 80°  is  critical  because  of  the  relatively  large  moment  of  inertia 
about  this  axis  (1^  )  and  the  high  rotation  speed  (see  Appendix  C).  [Ref.  37:  pp.  75-78] 

For  this  reason  a  limit  value  5  L  was  chosen.  If  the  satellite  is  in  an  orbit  with  5  > 
I  8  J  then  an  alternative  motion  will  be  executed.  The  yaw-axis  will  keep  a  constant 
attitude:  \\i  =  90°,  so  that  the  BAPTA-axis  is  perpendicular  to  the  orbital  plane.  In  this 
case  the  solar  arrays  would  rotate  360°/orbit  in  one  direction.  The  power  lines  would 
therefore  turn  around  the  BAPTA-axis  and  after  a  while  make  the  rotation  of  the  arrays 
impossible.  The  only  way  to  anticipate  this  is  rotating  them  in  the  opposite  direction 
during  the  eclipse.  This  alternative  motion  is  plotted  in  Figure  3-5.  [Ref.  37:  pp.  75-78] 
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Figure  3-3:  Motion  of  Yaw-axis  (vj/)  and  solar  array  (rj)  during  one  orbit,  starting  at  the 
sub-solar  point  (i.e.,  Noon)  in  terms  of  the  sun  rotation  angle. 
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Figure  3-4:  Motion  of  Yaw-axis  (v|/)  and  solar  arrays  (r\)  during  one  orbit,  starting  at  the 
sub-solar  point  in  terms  of  the  time  in  minutes. 
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Figure  3-5:  Alternative  motion  of  yaw-axis  (\\i)  and  solar  arrays  (r|)  during  one  orbit, 
starting  at  the  sub-solar  point,  for  |5|  >  5°. 
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IV.  QUANTIFY  DISTURBANCE  ENVIRONMENT 

In  this  part  of  the  design  we  must  understand  what  disturbances  are  acting  on  the 
spacecraft  and  then  try  to  quantify  them.  Disturbance  torques  can  be  either  external  or 
internal.  External  disturbance  torques  consist  of  gravity  gradient  effects,  magnetic  field 
torques  acting  on  the  vehicle,  the  effects  of  solar  radiation,  and  aerodynamic  effects. 
These  external  disturbance  torques  are  most  likely  the  largest  source  of  disturbance  to  the 
spacecraft  and  may  be  easier  to  quantify.  Internal  disturbance  torques  come  from  within 
the  spacecraft.  As  a  result,  we  have  a  great  amount  of  control  over  these.  [Ref.  6:  p.  352] 

In  this  section  an  investigation  of  both  external  and  internal  disturbance  torques 
will  be  performed.  And,  a  knowledge  of  both  magnitude  and  direction  of  disturbance 
torques  and  angular  momentum  will  be  attained.  This  information  will  be  very  useful  in 
selecting  and  sizing  the  ADCS  hardware.  [Ref.  6:  p.  352] 

A.  EXTERNAL  DISTURBANCE  TORQUES 

From  the  standpoint  of  stability  and  control,  the  principal  external  torques  come 
from  four  environmental  effects:  solar  radiation  pressure,  gravity  gradients,  aerodynamics 
and  magnetics.  Solar  radiation  pressure  is  generally  a  significant  source  of  attitude  and 
trajectory  errors  for  high  altitude  (>  1000  km).  Gravity  gradients,  which  result  from  the 
extended  dimensions  of  the  spacecraft,  may  either  cause  disturbing  or  perturbing  torques 
or  provide  restoring  torques  when  the  effect  is  used  for  attitude  control.  Gravity 
gradients,  as  well  as  the  magnetic  torques  caused  by  the  interaction  of  the  spacecraft's 
magnetic  materials  with  the  planetary  magnetic  field  are  most  significant  at  low  altitudes 
(<  1000  km).  Similarly,  aerodynamic  effects  are  significant  only  below  500  km  altitude 
and  are  generally  negligible  above  1000  km.  Finally,  aerodynamic  torques  acting  on  the 
spacecraft  are  functions  of  spacecraft  geometry,  attitude  and  altitude.  Figure  4-1 
illustrates  the  change  in  external  disturbance  as  a  function  of  altitude.  [Ref.  22:  pp.  77-78] 
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Figure  4-1 :  Typical  torques  on  a  satellite  as  a  function  of  altitude 

1.  Solar  Pressure  Torques 

The  solar  pressure  torque  is  the  major  long-term  disturbance  torque  for  a 
spacecraft  at  higher  altitudes.  The  solar  radiation  forces  are  due  to  photons  impinging  on 
the  spacecraft  surfaces.  In  general  a  fraction,  ps,  of  the  impinging  photons  will  be 
specularly  reflected,  a  fraction,  pd,  will  be  diffusely  reflected,  and  a  fraction,  pa,  will  be 
absorbed  by  the  surface  (see  Figure  4-2).  The  solar  pressure  moment  Ms  is  given  by 
equation  4-1.  [Ref.  25:  pp.  133-134] 
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Surface  Area  =  A 
Solar  Pressure  =  P 


Specularly       **~ 
reflected 
photons  <■ 


\       Incoming    <■ 
photons 


Figure  4-2:  Solar  Radiation  Force  on  a  Spacecraft  Surface 


Ms  =  PA(nS)rx 


(l-ps)S+2^ps  +  ^pdjn 


(4-1) 


where  r  is  the  vector  from  the  center  of  mass  of  the  spacecraft  to  the  center  of  pressure  of 
a  given  area,  A  .  The  solar  radiation  pressure,  P ,  is  generally  assumed  to  be  constant  and 
to  have  the  value  4  -644  x  l"6  N  /  m2 .  [Ref.  25:  pp.  133-134] 

2.  Gravity  Gradient  Torques 

A  spacecraft  body  experiences  a  gravity-gradient  torque  due  to  the  variation  of 
distances  between  the  spacecraft  mass  points  and  the  center  of  the  mass  of  the  earth.  The 
gravity  gradient  has  been  used  on  early  low  earth-orbit  satellites  to  maintain  the  earth 
pointing  of  antennas  or  other  instruments.  The  gravitational  torques  about  the  principal 
body  axes  d,  e2,  and  e3  are  given  in  equation  4-2.  [Ref.  25:  p.  131] 
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Tis)  =  -K(I33  -  722)cos(02)sin(03)sin(02)         =  0    for     G2  =  03  =  0 

T2(g)  =  K(IX ,  -  /33)cos(6>3) cos(02 ) sin(02 )  (4-2) 

r3(g)  =  -K(I22  -  /u)cos(03)[cos(02)]2  sin(03). 

3.  Aerodynamic  Torques 

The  satellite  will,  in  general,  pass  through  an  atmosphere  of  density  p ,  with  a 
velocity  v  .  The  magnitude  of  the  aerodynamic  force  F<a)is  then  given  as 

F^=±pv.vACd  (4-3) 

where  A  is  the  reference  area  of  the  satellite  (such  as  the  cross  section  along  v  )  and  Cd 
is  the  total  drag  coefficient.  Torque  is  given  by 


Pa)  =  ±-pv2lSCd  (4-4) 

2 


or 

T^=^pCdAv2{uxJcp)  (4-5) 


where  /  is  the  length  of  the  perpendicular  from  the  mass  center  to  the  force  line  of  action 
and  u  is  the  unit  vector  in  the  velocity  direction.  [Ref.  22:  p.  90] 

Difficulties  arise  in  determining  p,l,and  Cd.  At  satellite  altitudes,  p  is  highly 
dependent  on  the  time  of  day  and  the  level  of  solar  activity.  For  example,  at  600  km  the 
solar  daytime  maximum  density  and  the  nighttime  minimum  density  may  differ  by  a  factor 
of  100.  Table  4-1  gives  typical  values  of  the  daytime  maximum  air  density  as  a  function  of 
orbital  height.  It  is  important  to  note  that  as  altitude  increase  the  aerodynamic  pressure 
effects  decrease.  As  a  result,  our  satellite  orbit  of  11,000  km  will  have  negligible 
aerodynamic  pressure  effects.  [Ref.  22:  p.  90] 
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^able  4-1 :  Density  as  a  function 

of  height  above  the  earth's  surfa 

For  a  satellite  having  a  spherical  shape,  an  average  value  of  Cd  -  2.2  can  be  taken 
which  is  computed  assuming  "free  molecular  flow"  (i.e.,  the  molecular  mean  free  path  is 
assumed  to  be  large  compared  to  the  size  of  the  satellite  and,  therefore,  inter-particle 
collisions  are  ignored).  For  a  cylinder  Cd  may  be  taken  as  3.  [Ref.  22:  p.  90] 

4.  Disturbance  Torques  Due  to  the  Geomagnetic  Field 

a.  Geomagnetic  Field 

The  geomagnetic  field  can  be  represented  to  a  first  order  approximation  as 
a  magnetic  dipole  with  the  axis  inclined  to  the  Earth's  spin  axis  by  approximately  11.5°. 
Observations  have  shown  that  the  sources  of  the  geomagnetic  field  are  in  the  core  and  at 
the  surface  of  the  Earth  as  well  as  in  the  upper  atmosphere.  The  primary  source  is 
presumed  to  be  a  system  of  electric  currents  in  the  Earth's  molten  interior.  The  changes  in 
these  currents  account  for  the  migration  of  the  geomagnetic  poles  on  the  Earth's  surface. 
[Ref.  22:  p.  84] 

The  geomagnetic  field,  B ,  is  a  function  of  altitude,  longitude,  and  latitude. 
And,  can  be  expressed  in  units  of  gauss,  gamma  (y) ,  or  tesla  (webers  per  square  meter). 
A  one  gauss  magnetic  field  intensity  acting  on  a  unit  magnetic  pole  produces  a  force  on 
one  dyne.  The  magnitude  of  the  field  is  given  by  equation  4-6.  [Ref.  22:  p.  84] 
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|5j  =  ^f(l  +  3sin6g 


1/2 


(4-6) 


where  9m  is  the  magnetic  latitude  measured  from  the  geomagnetic  equator.  At  the  earth's 
surface  along  the  geomagnetic  equator,  6m  =  0 ,  and  B0  =  0.3 1 1  gauss.  At  the  magnetic 
pole,  B0  =  0.622  gauss  or  twice  the  equatorial  value.  [Ref.  22:  p.  84] 

A  plot  of  the  Earth's  equatorial  magnetic  field  intensity  B^  in  tesla  units  as 

a  function  of  altitude  is  shown  in  Figure  4-3.  The  effect  of  latitude  is  given  in  Figure  4-4, 
where  for  any  latitude  B0  =  aBeq .  [Ref.  22:  pp.  87-88] 
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Figure  4-3:  Earth's  Magnetic  field  intensity  vs  altitude  above  the  magnetic  equator. 
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Figure  4-4:  Multiplication  factor  a  for  magnetic  field  intensity  as  a  function  of  magnetic 

latitude.  For  any  latitude  B0  =  aBeqyaU)rial . 

b.  Magnetic  Torques 

Magnetic  torques  acting  on  a  spacecraft  can  result  from  the  interaction  of 
the  spacecraft's  residual  magnetic  field  and  the  geomagnetic  field.  Thus,  if  M  is  the  sum 
of  all  magnetic  moments  in  the  spacecraft  the  torque  acting  on  the  spacecraft  is 


f(m)  =  MxB 


(4-7) 


where  B  =  the  geomagnetic  field  vector.  In  general,  M  can  be  caused  by  permanent  and 
induced  magnetism  or  by  spacecraft-generated  current  loops.  This  is  how  Magnetic 
torque  rods  operate  and  can  be  used  for  momentum  dumping.  The  units  of  M  may  be 
gauss-cm3 ,  ampere-m2 ,  or  pole-cm.  For  example,  if  Mis  in  ampere-m2 ,  then  Pm)  is  in 
Newton-meters.  [Ref.  22:  p.  89] 


49 


The  concept  of  Magnetic  Torquers  is  simple,  but  there  is  a  difficulty  in 
quantifying  external  disturbance  torques  due  to  the  geomagnetic  field.    M  is  the  magnetic 
field  produced  by  the  spacecraft  and  is  generally  unknown  until  the  spacecraft  is  built  and 
tested.  Once  tested,  an  estimate  for  the  magnetic  disturbance  torque  can  be  made.  [Ref. 
22:  p.  89] 

B.       MATLAB    SIMULATION    OF    EXTERNAL    DISTURBANCE 
TORQUES 

To  attain  an  accurate  quantifiable  value  for  external  disturbance  torques  a  matlab 
simulation  was  developed  for  each  of  the  external  disturbances.  A  brief  explanation  of 
each  scenario  follows: 

1.  Solar  Pressure  Torques 

1.  Developed  the  Simulink  model  appropriate  for  the  Solar  Pressure  Simulation 
(see  Appendix  E). 

2.  Established  the  initial  conditions  for  the  velocity  and  position  vector  of  the 
spacecraft  relative  to  the  earth  centered  inertial  system  (see  Appendix  E). 

3.  Defined  various  constants  and  solved  for  v  and  8  as  a  function  of  time,  v  is  the 
earth  centered  true  anomaly  and  5  is  the  sun  orbital  element  analogous  to  v  (see  Appendix 
E). 

4.  Established  numerous  direction  cosine  matrices.  This  part  of  the  process  was 
very  tedious.  The  details  are  excluded  because  it  is  assumed  the  reader  is  familiar  with  the 
various  orbit  frames  of  reference  (see  Appendix  E). 

5.  Computed  the  unit  vector  from  the  spacecraft  to  the  sun  (see  Appendix  E). 

6.  Defined  unit  vectors  that  are  perpendicular  to  each  exterior  piece  of  the 
spacecraft  (see  Appendix  E). 

7.  Defined  dr  vectors  from  the  center  of  mass  to  each  external  piece's  center  of 
pressure  (see  Appendix  E). 

8.  Performed  various  vector  operations  (see  Appendix  E). 
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9.  Calculated  the  moment  due  to  solar  pressure  as  given  by  equation  4-1  [Ref.  25: 
p.  133]: 


Ms  =  PA{n*S)r  x 


(l  +  ps)S  +  2(ps  +  [^pdn 


where, 

Ms  =   The  solar  radiation  pressure  moment. 

f  =  The  vector  from  the  center  of  mass  of  the  spacecraft  to  the 
center  of  pressure  of  a  given  area,  A. 

P  =   The  solar  radiation  pressure. 

A  =   Area 

ps  =  Fraction  of  photons  impinging  on  the  spacecraft's  surface 
that  are  specularly  reflected. 

pd  =  Fraction  of  photons  impinging  on  the  spacecraft's  surface 
that  are  diffusely  reflected. 

1 0.  Computed  the  M, ,  M2 ,  and  M3  for  solar  pressure  torques  for  each  instant  in 

time  (see  Appendix  A- 10  to  A- 17). 

1 1 .  Solved  for  angular  velocities  iteratively  (see  Appendix  E). 

12.  Utilized  Poisson's  equations  to  solve  for  the  dcm  components  from  the  body 
frame  to  the  orbit  frame  iteratively  (see  Appendix  E). 

13.  Solved  for  r  and  v  iteratively  utilizing  state  space  variables  (see  Appendix  E). 
Figures  4-5  through  4-7  are  the  result  of  the  calculations  described  above.  A  brief 

explanation  of  each  figure  follows: 
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a.  Figure  4-5 

Figure  4-5  is  a  graph  of  the  three  components  of  torque,  (i.e., 
TX,T ,  and  Tz )  acting  on  the  spacecraft  due  only  to  solar  pressure  disturbances,  when  the 

spacecraft  is  not  being  controlled  by  the  Attitude  Dynamics  and  Control  Subsystem.  Also, 
the  torque  components  are  with  respect  to  the  spacecraft's  body  fixed  frame. 

b.  Figure  4-6 

Figure  4-6  is  a  graph  of  the  three  components  of  the  spacecraft's  angular 
velocity  due  only  to  solar  pressure  disturbances  when  the  spacecraft  is  not  being 
controlled  by  the  Attitude  Dynamics  and  Control  Subsystem.  Also,  the  angular  velocity 
components  are  expressed  with  respect  to  the  spacecraft's  body  fixed  frame. 

c.  Figure  4-7 

Figure  4-7  is  a  graph  of  the  three  components  of  torque,  (i.e., 
Tx,  Ty,  and  Tz)  acting  on  the  spacecraft  due  only  to  solar  pressure  disturbances,  when  the 

spacecraft  is  being  controlled  by  the  Attitude  Dynamics  and  Control  Subsystem  and  is 
maintaining  a  nadir  pointing  configuration.  Also,  the  torque  components  are  with  respect 
to  the  spacecraft  body  fixed  frame. 

It  should  be  noted  that  the  solar  radiation  torque  component  along  the 
pitch  axis  is  periodic  and  secular,  with  a  period  equal  to  the  orbit  period.  Hence  the  net 
effect  of  this  torque  is  a  build  up  of  Angular  Momentum.  On  the  other  hand,  the 
components  along  the  other  axes  have  cyclic  components,  resulting  in  no  net  change  in  the 
spacecraft  angular  momentum  over  a  full  orbit. 
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Figure  4-5:  The  three  components  of  torque  due  to  solar  pressure  as  a  function  of  time, 

with  no  control. 
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Figure  4-6:  The  three  body  components  of  angular  velocity  due  to  solar  pressure  as  a 

function  of  time,  with  no  control. 
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Figure  4-7:  The  three  components  of  torque  due  to  solar  pressure  as  a  function  of 

time,  with  active  control. 

2.  Gravity  Gradient  Torques 

1 .  Developed  the  simulink  diagram  (see  Appendix  E). 

2.  Defined  constants  and  user  friendly  conversions  (see  Appendix  E). 

3.  Determined  initial  conditions  for  r,  v ,  and  the  dcm  between  orbit  frame  and 
principal  inertia  frame  (see  Appendix  E). 

4.  Solved  for  angular  velocities  iteratively  (see  Appendix  E). 

5.  Utilized  Poisson's  equations  to  solve  for  the  dcm  components  from  the  body 
frame  to  the  orbit  frame  iteratively  (see  Appendix  E). 

6.  Solved  for  r  and  v  iteratively  utilizing  state  space  variables  (see  Appendix  E). 

7.  Calculated  the  moment  due  to  Gravity  Gradient  as  given  by  equation  4-11. 

^=K(I33-I22)a2^ 

T^=K(IU-I^)aua„ 
T^=K(122-Iu)aua„ 
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Figures  4-8  and  4-9  are  the  result  of  the  calculations  described  above.    A  brief 
explanation  of  each  figure  follows: 

a.  Figure  4-8 

Figure    4-8    is    a    graph    of   the    three    components    of  torque,    (i.e. 

Tx,Ty,  and  Tz)  acting  on  the  spacecraft  due  only  to  gravity  gradient  disturbances,  when 

the  spacecraft  is  not  being  controlled  by  the  Attitude  Dynamics  and  Control  Subsystem. 

Also,  the  torque  components  are  with  respect  to  the  spacecraft's  body  fixed  frame  and 

increase  over  time  such  that  they  cause  the  spacecraft  to  become  unstable. 

b.  Figure  4-9 

Figure  4-9  is  a  graph  of  the  three  components  of  the  spacecraft's  angular 
velocity  due  only  to  gravity  gradient  disturbances  when  the  spacecraft  is  not  being 
controlled  by  the  Attitude  Dynamics  and  Control  Subsystem.  Also,  the  angular  velocity 
components  are  expressed  with  respect  to  the  spacecraft's  body  fixed  frame 
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Figure  4-8:   The  three  components  of  torque  due  to  gravity  gradients  as  a  function  of 

time,  with  no  control 
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Figure  4-9:  The  three  body  components  of  angular  velocity  due  to  gravity  gradients  as  a 

function  of  time,  with  no  control. 

3.  Aerodynamic  Pressure  Torques 

1 .  Calculating  the  Aerodynamic  Pressure  Torques  was  very  similar  to  calculating 
solar  pressure  torques.  Nearly  the  same  technique  was  used  except  for  the  following  (see 
Appendix  E). 

2.  The  fundamental  equation  used  to  compute  is  given  by  equation  4-5  [Ref.  6:  p. 
309]: 


M  =  [^jpCdAv2(uvxIcp) 


where. 


M  =   Torque  due  to  Aerodynamic  Pressure 

p  =   Atmospheric  density 

Cd  =   Coefficient  of  drag 

A  =   Area  perpendicular  to  uv 
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uv  =   Unit  vector  in  velocity  direction 
Icp  =     Vector  distance  from  center  of  mass  to  center-of- 
pressure 

3.  Computed  the  density  (p)  as  a  function  of  height  above  the  earth  (see 
Appendix  E). 

4.  Computed  the  cross  sectional  areas  of  the  spacecraft  in  the  three  body  frame 
coordinate  axis  directions  (see  Appendix  E). 

5.  Used  the  same  Euler  rates  and  Poisson's  equations  to  solve  iteratively  (see 
Appendix  E). 

6.  Utilized  state  space  variables  to  solve  for  r  and  v  iteratively  (see  Appendix  E). 
Figures  4-10  through  4-12  are  the  result  of  the  calculations  described  above.    A 

brief  explanation  of  each  figure  follows: 

a.  Figure  4-10 

Figure  4-10  is  a  graph  of  the  three  components  of  torque,  (i.e., 
Tx,  Ty,  and  Tz)  acting  on  the  spacecraft  due  only  to  aerodynamic  pressure  disturbances, 

when  the  spacecraft  is  not  being  controlled  by  the  Attitude  Dynamics  and  Control 
Subsystem.  The  torque  components  are  with  respect  to  the  spacecraft's  body  fixed  frame. 
Also,  these  torque  magnitudes  are  orders  of  magnitude  less  than  gravity  gradient  and  solar 
pressure  disturbances.  And,  catastrophic  cancellation  (i.e.,  computer  noise)  most  likely 
occurred  when  calculating  Ty .   This  was  probably  due  to  the  fact  that  the  torque  values 

were  so  small.  As  a  result,  the  torques  due  to  aerodynamic  pressure  are  inconsequential. 

b.  Figure  4-11 

Figure  4-11  is  a  graph  of  the  three  components  of  the  spacecraft's  angular 
velocity  due  only  to  aerodynamic  pressure  disturbances  when  the  spacecraft  is  not  being 
controlled  by  the  Attitude  Dynamics  and  Control  Subsystem.  The  angular  velocity 
components  are  expressed  with  respect  to  the  spacecraft's  body  fixed  frame.    Also,  the 
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aerodynamic  torques  are  orders  of  magnitude  less  than  gravity  gradient  and  solar  pressure 
disturbances.  Therefore,  the  effects  are  inconsequential. 

c.  Figure  4-12 

Figures   4-12   is   a   graph   of  the   three    components   of  torque,    (i.e., 
Tx,  Ty,and  Tz)  acting  on  the  spacecraft  due  only  to  aerodynamic  pressure  disturbances, 

when  the  spacecraft  is  being  controlled  by  the  Attitude  Dynamics  and  Control  Subsystem 
and  maintaining  a  nadir  pointing  configuration.  The  torque  components  are  with  respect 
to  the  spacecraft  body  fixed  frame.  Also,  the  effects  of  aerodynamic  pressure  disturbances 
are  orders  of  magnitude  less  than  gravity  gradient  and  solar  pressure  disturbances.  Also, 
catastrophic  cancellation  (i.e.,  computer  noise)  may  have  occurred.  As  a  result,  the 
torques  due  to  aerodynamic  pressure  are  inconsequential. 
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Figure  4-10:  The  three  components  of  torque  due  to  Aerodynamic  Pressure  as  a  function 

of  time,  with  no  control. 
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x  10'2B  Angular  Velocities  from  Aerodynamic  Pressure  Torques  only 
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Figure  4-11:  The  three  body  components  of  angular  velocity  due  to  Aerodynamic 
Pressure  as  a  function  of  time,  with  no  control 
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Figure  4-12:   The  three  components  of  torque  due  to  aerodynamic  pressure 
disturbances  as  a  function  of  time,  with  active  control. 
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4.  Combined  External  Disturbance  Torques 

Next,  calculations  were  made  for  when  all  three  torques  acted  on  the  spacecraft 
simultaneously.  This  required  combining  the  three  separate  Matlab  programs  into  one 
program.  The  motion  of  the  spacecraft  was  then  simulated  and  the  sum  of  the 
environmental  disturbance  torques  computed  (see  Appendix  E). 

Finally,  an  estimation  of  the  total  torque  imparted  to  the  spacecraft  over  one  year 
for  each  of  the  four  situations  was  performed.  To  do  this  for  gravity  gradient  and 
aerodynamic  required  taking  data  from  several  orbits,  then  computing  the  average  torque 
over  that  period  and  finally  multiplying  by  the  time  expired.  This  method  is  good  for 
gravity  gradient  and  aerodynamic  as  long  as  data  from  several  orbits  is  computed. 
Computing  the  total  torque  from  the  solar  pressure  is  more  complicated.  Since  the  Beta 
angle  changes  with  time  throughout  the  year,  several  orbits  will  not  be  adequate.  So,  the 
simulation  was  ran  for  a  delta  t  of  one  year,  then  the  total  torque  was  computed.  The 
results  can  be  seen  in  Table  4-2. 


SOLAR 
PRESSURE 

GRAVITY 
GRADIENT 

AERO- 
DYNAMIC 

COMBINED 

Average  Tx 

9.33e-10Nm 

*0 

*0 

-2.5e-6  Nm 

Average  T2 

1.33e-6Nm 

1.2e-6Nm 

9.55e-29  Nm 

1 .4e-4  Nm 

Average  T3 

5.03e-8  Nm 

«0 

6.6e-ll  Nm 

-2.78e-6  Nm 

Maximum  Torque 
(T     ) 

V     max  / 

1.8e-6Nm 

le-4Nm 

7e-llNm 

2.2e-4  Nm 

Estimated  //,  (1  year) 

0.0296  Nms 

*0 

*0 

0.1718  Nms 

Estimated  H2  (1  year) 

42.23  Nms 

3155  Nms 

3.03e-21  Nms 

3197.23  Nms 

Estimated  //3  (1  year) 

1.594  Nms 

*0 

0.0021  Nms 

0.19033  Nms 

Table  4-2:  Summary  of  External  Disturbance  Torques  and  Angular  Momentum  Imparted 

to  the  Spacecraft 


Figures  4-13  through  4-15  are  the  result  of  the  calculations  described  above.    A 
brief  explanation  of  each  figure  follows: 
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a.  Figure  4-13 

Figure  4-13  is  a  graph  of  the  three  components  of  torque,  (i.e. 
Tx,Tv,and  Z)  acting  on  the  spacecraft  due  to  solar  pressure,  gravity  gradient  and 

aerodynamic  disturbances,  when  the  spacecraft  is  not  being  controlled  by  the  Attitude 
Dynamics  and  Control  Subsystem.  Also,  the  torque  components  are  with  respect  to  the 
spacecraft's  body  fixed  frame. 

b.  Figure  4-14 

Figure  4-14  is  a  graph  of  the  three  components  of  the  spacecraft's  angular 
velocity  due  to  solar  pressure,  gravity  gradient  and  aerodynamic  pressure  disturbances 
when  the  spacecraft  is  not  being  controlled  by  the  Attitude  Dynamics  and  Control 
Subsystem.  Also,  the  angular  velocity  components  are  expressed  with  respect  to  the 
spacecraft's  body  fixed  frame. 

c.  Figure  4-15 

Figures  4-15  is  a  graph  of  the  three  components  of  torque,  (i.e. 
Tx,Ty,and  Tz)  acting  on  the  spacecraft  due  to  solar  pressure,  gravity  gradient  and 

aerodynamic  pressure  disturbances,  when  the  spacecraft  is  being  controlled  by  the  Attitude 
Dynamics  and  Control  Subsystem.  Also,  the  torque  components  are  with  respect  to  the 
spacecraft  body  fixed  frame. 

It  should  be  noted  that  the  total  torque  component  along  the  pitch  axis  is 
periodic  and  secular,  with  a  period  equal  to  the  orbit  period.  Hence  the  net  effect  of  this 
torque  is  a  build  up  of  Angular  Momentum.  On  the  other  hand,  the  components  along  the 
other  axes  have  cyclic  components,  resulting  in  no  net  change  in  the  spacecraft  angular 
momentum  over  a  full  orbit. 
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Figure  4-13:  The  three  components  of  torque  due  to  total  external  disturbance  torques  as 

a  function  of  time,  with  no  control. 
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Figure  4-14:  The  three  body  components  of  angular  velocity  due  to  total  external 
disturbance  torques  as  a  function  of  time,  with  no  control. 
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Figure  4-15:  The  three  components  of  torque  due  to  solar  pressure,  gravity 
gradients  and  aerodynamic  pressure  disturbances  as  a  function  of  time,  with  active  control. 

C.  INTERNAL  DISTURBANCE  TORQUES 

There  are  several  other  kinds  of  disturbance  torques  that  are  internal  to  the 
spacecraft.  First,  these  torques  may  be  due  to  moving  parts  such  as  rotating  wheels, 
circulating  fluids  or  tracking  devices.  Fortunately,  the  ADCS  designer  has  some  control 
over  these  torques.  For  example,  if  a  disturbance  torque  in  one  component  is  found  to  be 
much  larger  than  the  others,  the  component  can  be  changed  in  order  to  meet  higher 
specifications.  However,  there  is  a  trade-off  here.  A  change  can  be  made  to  meet  higher 
standards  only  at  the  expense  of  more  mass  and/or  cost.  Next,  an  internal  torque  may 
arise  due  to  thruster  mis-alignment  with  the  center  of  gravity.  This  will  only  occur  during 
thruster  firing  and  will  be  corrected  in  a  closed  loop  Attitude  Dynamics  and  Control 
Subsystem.  The  concern  here  is  that  the  disturbance  will  not  cause  the  spacecraft  to 
exceed  the  minimum  pointing  requirements.  [Ref.  6:  p.  352] 

Of  greater  importance  is  slosh  and  mechanical  machinery  torques.  They  can, 
however,  be  controlled  by  proper  selection  of  hardware.  Table  4-3  summarizes  the 
common  internal  disturbance  torques.   [Ref.  6:  p. 3 52] 
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Disturbances 

Effect  on  Vehicle 

Typical  Values 

Uncertainty  in  center  of 
gravity  (eg) 

•  Unbalanced  torques  during 
firing  of  coupled  thrusters 

•  Unwanted  torques  during 
translation  thrusting 

1  to  3  cm 

Thruster  misalignment 

Same  as  eg  uncertainty 

0.1°  to  0.5° 

Mismatch  of  thruster 
outputs 

Similar  to  eg  uncertainty 

±5% 

Rotating  machinery  (pumps, 
tape  recorders) 

Torques  that  perturb  both 
stability  &  accuracy 

Dependent  on  spacecraft 
design;  may  be 
compensated  by  counter- 
rotating  elements 

Liquid  sloshing 

Torques  due  to  fluid  motion 
and  variation  in  center  of  mass 
location 

Dependent  on  specific 
design;  may  be  controlled 
by  bladders  or  baffles 

Dynamics  of  flexible  bodies 

Oscillatory  resonance  at 
bending  frequencies,  limiting 
control  bandwidth 

Depends  on  spacecraft 
structure 

Thermal  shocks  on  flexible 
appendages 

Attitude  disturbances  when 
entering/leaving  eclipse 

Depends  on  spacecraft 
structure.  Worst  for 
gravity  gradient  systems 
with  long  inertia  booms 

Table  4-3 :  Principal  Internal 

disturbance  Torques  which  space 

craft  designers  have  a  great 

deal  of  control  over. 
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V.  ATTITUDE  DYNAMICS  AND  CONTROL  SUBSYSTEM 
HARDWARE  SELECTION  AND  DESCRIPTION 

A.  SIZING  THE  ADCS  REACTION  WHEELS 

1.  Torque  and  Angular  Momentum  Requirements  Due  to  Yaw 
Steering 

In  Chapter  III  the  yaw  steering  requirements  of  the  spacecraft  were  investigated 
and  it  was  concluded  the  Yaw  command  angle  could  be  stated  as: 


y/(P,S)  =  tan 


-i 


-sinQg) 

.  Xan(S) 


(5-1) 


where. 


P  is  the  angle  in  the  orbital  plane  between  the  point  with  minimum  distance  (sub- 
solar  point)  to  the  Sun  and  the  satellite.  And,  5  is  the  angle  between  the  orbital  plane  and 
the  direction  of  the  Sun  (see  Figure  5-1). 

By  looking  at  Figure  3-4  it  is  apparent  that  the  angular  velocity  and  angular 
acceleration  maximums  occur  for  the  minimum  value  of  8.  And,  since  our  minimum  8  is 
equal  to  5°  we  will  single  this  curve  out  and  find  the  maximum  angular  velocity  and 
maximum  angular  acceleration  the  spacecraft  must  perform  in  the  yaw  direction.  Since 
this  requirement  is  the  most  demanding;  if  we  can  show  it  will  be  met  then  any  other 
torque  and  angular  momentum  requirement  will  also  be  satisfied. 


Figure  5-1  is  a  graph  of  ^(/?,5°)  =  tan 


-sinpg) 
tan(5°)  _ 


for  8  =  5°.  It  can  readily  be 


shown  that  y/  (0,5° )  which  is  the  angular  velocity  is  given  by 


y/\P$)*  (-0.1 806) 


cos(2.76e-4*r) 


1  +  130.646  *  [sin(2.76e  -  4  *  t)] 


(5-2) 


and,  y/  (/?,5° )  which  is  the  angular  acceleration  is  given  by 
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^■■(A5-)*(4.98e-5)- 


sin(2.7fc?-4*0 


1  +  1 30.646  *  [sin(2.76e  -  4  *  /)]2 


+ 


(5-3) 


(0.013) 


[cos(2.76e  -  4  *  Q]2  *  sin(2.756g  -  4  *  Q 
1  +  (1 30.646)  *  [sin(2.76e  -  4  *  /)] 


Figures  5-2  and  5-3  are  graphs  of  the  angular  velocity  and  angular  acceleration 
respectively  for  8  =  5°. 


Figure  5-1 :  Graph  of  y/(/3,5°) ,  where  vy  is  in  degrees  and  time  is  in  seconds. 
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Figure  5-2:  Angular  Velocity  in  units  of  degrees/second. 
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Figure  5-3:  Angular  Acceleration  in  units  of  degrees/second 
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From  the  graphs  it  is  observed  that  the  maximum  angular  velocity  occurs  at  t=0, 
11400  and  22800  seconds  and  the  maximum  angular  acceleration  occurs  at  t  =  0+, 
1 1400  ~ ,  1 1400  +  ,  and  22800  "  seconds.  The  respective  values  are  as  follows: 

y/(/?,5°)max  *  0.18°  /  sec  *  3.1 42e-  3  rad  I  sec 

¥  (#5°  )max  *  0.00037°  /  sec2  a  6.4577e  -  6  rarf  /  sec2 

The  maximum  acceleration  of  the  satellite  has  to  be  realized  by  accelerating  the 
reaction  wheel  in  the  opposite  direction.   The  maximum  torque  T  max  needed  follows  from 

the  formula: 

^max  =  Izz*v"™*  =  041 13)  *(6.4577e  -6  rad  I  sec2)  =  0.091 1     Nm  (5-4) 

And,  the  maximum  angular  momentum  B  max  is  given  by: 

A™  =  izz*¥'™*  =  (141 13)*(3.142<?  -  3  rad  I  sec)  =  44.34     Nms  (5-5) 

The  total  acceleration  Aa>  of  the  reaction  wheel  due  to  the  rotation  can  simply  be 
calculated  with  the  following  relation: 

Izz*y/m**.  =  I<o*£i(0  =>  Act)  =   zz   - max  (5-6) 
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where,  we  find  Im  from  the  formula  for  the  moment  of  inertia  for  a  hollow  cylinder 

with  a  relative  small  wall  thickness,  /  =  mR2 .    So,  for  a  reaction  wheel  with  a  mass  of 
1 0.5  kg  and  a  radius  of  0. 1 965m: 

=>    Im  «  0.4054     kg  m2 


such  that. 


Aco  =  /zz*^max  *  9.69  rad  I  sec  a  93  rpm 
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In  conclusion,  the  reaction  wheel  to  be  selected  must  meet  these  requirements  for 
torque,  momentum  and  acceleration. 

2.  Torque  and  Angular  Momentum  Requirements  Due  to  External 
Disturbances 

The  maximum  torque  and  angular  momentum  requirements  due  to  external 
disturbances  were  calculated  and  found  to  be: 

Tlmax  =  Se-4     Nm  Himax  =  0.1718     Nms 

T2max=22e-3    Nm  H2max  =  9A9S    Nms 

Timax  =  6e-4    Nm  H3max  =  0.19033     Nms 

As  can  be  seen,  these  torque  and  angular  momentum  requirements  are  not  as 

demanding  as  those  resulting  from  yaw  steering.     However,  the  external  disturbances 

produce  a  secular  torque  which  builds  up  angular  momentum  with  time.  As  can  be  seen 
from  Figure  4-16;  a  secular  torque  results  in  the  ^direction  and  creates  a  buildup  of 

angular  momentum  as  given  by 

(  Tvavg )  x  (1  month)  =  362  Nms 

over  the  course  of  one  month.  This  build  up  must  be  dumped  by  either  the  Magnetic 
Acquisition  /  Despin  system  or  the  spacecraft  thrusters.  Not  accounting  for  these  secular 
disturbances  appropriately  could  result  in  a  premature  mission  due  to  unexpected  demands 
of  spacecraft  fuel. 

B.  SIZING  THE  MAGNETIC  TORQUERS 

Magnetic  torquers  use  the  Earth's  magnetic  field,  B ,  and  electrical  current  through 
the  torquer  to  create  a  magnetic  dipole  ( D)  that  results  in  torque  (  T )  on  the  vehicle  given 
by: 

D-T-  (5-7) 

Magnetic  torquers  used  for  momentum  dumping  must  equal  the  peak  disturbance 
plus  a  margin  to  compensate  for  the  lack  of  complete  directional  control. 

From  equation  5-8  we  can  estimate  the  worst  case  earth  magnetic  field  for  our 
orbit  given  by  using  the  following  values: 
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r  =  1.737  el  m 

0  =  0° 

juE  =  8.1  e25  g-cm3  =  8.1  e\5  Tm3 


|5(U^f(l  +  3sin^)1/2  (5-8) 

ii      r 


«     \  =  \55e-6    T 

max 


Next,  we  calculate  the  torque  rod's  magnetic  torquing  ability  (dipole)  to  counteract 
the  worst  case  external  disturbance; 


as, 


T     =  2.2e-3     Nm 

max 


D  =  —  =  1419.35     A-m2 


B 


which  is  achievable.  The  earth's  field  is  cyclic  at  twice  the  orbital  frequency;  thus, 
maximum  torque  is  available  only  twice  per  orbit.  A  torquer  of  2500  to  3000 
A-m2 capacity  will  provide  sufficient  margin  to  counteract  these  expected  maximum 
disturbance  torques  and/or  provide  a  momentum  dumping  capacity  for  the  reaction 
wheels. 

C.  OTHER  SENSOR  SELECTION  ISSUES 

The  driving  requirement  for  the  remaining  sensors  and  actuators  are  the  attitude 
knowledge  and  pointing  requirements  as  stated  in  Chapter  6  under  the  system  level 
requirements  section.  All  sensors  and  actuators  selected  meet  these  requirements  with  a 
considerable  margin  built  in. 

D.  HARDWARE  SELECTION 

1.  Spacecraft  Control  Processor 

The  Satellite  Control  Processor  unit  selected  is  produced  by  Honeywell  and  is 
scheduled  to  be  used  on  the  International  Space  Station.  Honeywell's 
Multiplexer/Demultiplexer  (MDM)  is  a  high  speed  interface  and  computer  that  controls 
sensors,  effectors  and  other  unique  devices.  [Ref.  29] 
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The  MDM  collects  signals  from  sensors,  effectors  and  subsystems;  then  processes 
the  information  and  transmits  signals  and  commands.  For  example,  the  MDM  directs  the 
solar  arrays  to  point  at  the  sun.  [Ref.  29] 

The  MDM  is  rugged  and  adaptable  to  a  wide  range  of  applications.  Five  standard 
I/O  circuit  cards  cover  most  applications,  but  this  can  be  modified  if  needed.  Also,  the 
MDM  size  can  be  modified  based  on  the  number  of  I/O  circuit  cards  needed.  [Ref.  29] 

To  enhance  reliability  and  supportability,  software  and  hardware  self-test  features 
can  detect  and  report  94%  of  all  possible  failure  modes.  In  addition  Honeywell  can  help 
develop  I/O  cards  and  user  application  programs  that  meet  individual  needs.  Figure  V-4  is 
an  overview  of  the  MDM  architecture.  Features  of  the  MDM:  [Ref.  29] 

1 .  The  basic  MDM  consists  of  a  chassis,  power  supply,  I/O  Control  Unit,  I/O 
backplane,  SX  backplane  and  I/O  circuit  cards  selected  to  meet  specific  requirements. 

2.  The  I/O  Control  Unit  includes  the  processor,  memories,  1553  bus  interface 
function,  I/O  Circuit  Card  Control  logic,  and  Analog-to-Digital  converter. 

3.  Communication  is  via  independent,  dual  redundant  MIL-STD-1553B  channels, 
enabling  the  MDM  to  operate  in  Bus  Controller,  Bus  Monitor  or  Remote  Terminal  mode. 
Firmware  supports  functions  of  each  1553  channel  independently  of  the  user  application 
software.  The  MDM  is  capable  of  serial  asynchronous  communication  via  RS449  with 
RS485  driver/receivers  and  provides  1 6  bits  of  parallel  I/O  functionality.  The  MDM  self- 
test  capabilities  include  Built-in  Test  (BIT)  software  routines  and  supporting  hardware 
Built-in  Test  Equipment  (BITE). 
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2.  Magnetic  Acquisition  and  De-spin  System 

The  Magnetic  Acquisition  and  Despin  (MAD)  system  is  designed  for  use  on 
spacecraft  for  the  reduction  of  spin  rates  and  momentum  dumping  and  is  built  and 
manufactured  by  ITHACO  Space  Systems.  The  spin  rates  can  be  encountered  as  a 
normal  part  of  a  mission  or  can  result  from  ground  control  errors  or  onboard  equipment 
failures.  Under  normal  conditions,  the  launch  vehicle  leaves  the  satellite  spinning  at  the 
end  of  orbit  insertion.  If  the  spacecraft  is  to  be  stabilized,  the  spin  rate  imparted  by  the 
launch  vehicle  must  be  corrected.  The  unexpected  activation  of  a  thruster  or  a  momentum 
wheel  failure  would  also  result  in  unacceptable  spin  rates.  Corrective  action  can  often  be 
taken  if  the  vehicle  rates  can  be  reduced  in  some  manner.  The  Magnetic 
Acquisition/Despin  system  performs  this  function.  [Ref.  30,  Ref.  31] 

The  MAD  system  consists  of  three  Torque  rods,  an  electronics  controller,  and  a 
three-axis  magnetometer.  One  Torque  rod  is  aligned  with  each  axis  of  the  spacecraft  and 
is  mounted  away  from  any  instruments  sensitive  to  magnetic  fields.  In  addition,  Torque 
rods  are  separated  from  one  another  in  order  to  avoid  cross  coupling  effects.  The  use  of 
Torque  rods  allows  complete  momentum  management  of  the  spacecraft  without 
consumables  in  almost  any  orbit,  thus  providing  a  clean,  contamination  free  environment. 
Even  in  high  orbits,  Torque  rods  can  make  a  significant  contribution  to  the  total 
momentum  management  problem.  By  eliminating  or  reducing  the  need  for  consumables, 
the  use  of  Torque  rods  can  reduce  mission  costs  significantly.  Next,  the  electronics 
controller  can  be  mounted  at  any  convenient  location.  The  electronics  should  be 
accessible,  however,  since  a  trimming  operation  at  the  spacecraft  level  is  required  to 
compensate  for  the  unavoidable  magnetometer  to  torque  rod  coupling.  The 
magnetometer  should  be  mounted  so  that  one  Torque  rod  is  aligned  with  each  of  the 
magnetometer  axes.  The  magnetometer  probe  should  be  mounted  as  far  away  from  the 
Torque  rods  as  possible  to  niinimize  magnetic  coupling.  The  maximum  rate  from  which 
the  MAD  system  can  achieve  capture  depends  upon  the  bandwidth  that  can  be  used.  The 
maximum  bandwidth  in  any  given  system  depends  on  the  amount  of  coupling  between  the 
Torque  rods  and  the  magnetometer.  Thus,  maximum  separation  is  an  advantage  for  the 
available  capture  capability  of  the  MAD  system.  [Ref.  30,  Ref.  31] 
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The  MAD  system  relies  upon  a  "minus  B-dot"  control  law.  Each  Torque  rod  is 
driven  to  a  value  determined  by  the  negative  rate  of  change  of  the  magnetic  field  along  the 
rod.  The  effect  is  that  the  Torque  rods  behave  like  magnetic  brakes  and  lock  onto  the 
Earth's  magnetic  grid.  From  this  terminal  state,  the  spacecraft's  attitude  control  system 
can  capture  control  or  excess  momentum  can  be  dumped.  The  system  specifications  are 
as  follows: 


SIZE 

WEIGHT 

POWER 

Control  Electronics 

6"  x  6"  x  2" 

0.9  kg 

2.0  W 

Magnetometer 

2.8"  x  3.8"  x  4.7" 

0.5  kg 

0.7  W 

Torque  Rods 

19.5  " 

0.9  kg 

N/A 

Table  5-1 :  Specifications  of  Magnetic  Acquisition  /  Despin  System 
3.  Earth  Sensor  /  Sun  Sensor  Assembly 

a.  Overview 

The  Earth  Sensor/Sun  Sensor  Assembly  (ES/SSA)  is  an  innovative 
combined  Earth  horizon  and  two  axis  sun  sensor  which  offers  unprecedented  performance 
in  terms  of  accuracy,  mass,  power  consumption,  reliability,  and  automation.  This  sensor 
autonomously  determines  three-axis  attitude  information  at  medium  earth  orbits, 
geosynchronous,  or  super  geosynchronous  altitudes,  as  well  as  during  transfer  orbit.  [Ref. 
32] 

b.  Theory  of  Operation 

The  Earth  Sensor  portion  of  the  ES/SSA  is  a  scanning  sensor  with  two  sets 
of  two  infrared  (IR)  Fields  of  View  which  scan  about  the  axis  of  rotation  (scan  axis)  at  a 
rate  of  2  revolutions  per  second.  Each  of  the  four  IR  scans  possesses  a  1.5°  instantaneous 
Field  of  View.  As  the  scan  progresses,  each  of  the  instantaneous  Fields  of  View  crosses 
the  relative  "cold"  of  space  to  the  "warmth"  of  Earth  and  back  into  space.  Roll  attitude  is 
derived  from  differences  in  the  sensed  Earth  chord  widths.  Pitch  attitude  is  derived  from 
the  phase  angle  of  the  sensed  Earth  crossings.  [Ref.  32] 

The  sensor  has  also  been  designed  to  operate  in  transfer  orbit  and/or 
constant  altitudes  ranging  from  11,000  km  through  77,000  km.     At  geosynchronous 
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altitudes  all  four  Earth  scans  cross  the  Earth's  disc  at  nominally  3  and  6  degrees  above  the 
equator  and  3  and  6  degrees  below  the  equator.  [Ref.  32] 

The  Sun  Sensor  portion  of  the  ES/SSA  incorporates  two  static,  mutually 
orthogonal  sensing  (axes)  arrays  to  detect  the  orientation  of  the  sun.  Each  axis  is 
comprised  of  a  slit  allowing  a  plane  of  light  to  pass  through  and  intersect  with  an 
orthogonally  mounted,  high  resolution,  linear  photodiode  array.  The  point  of  intersection, 
which  indicates  sun  angle,  is  obtained  by  reading  the  energy  stored  in  the  array.  [Ref.  32] 

Each  ES/SSA  is  comprised  of  a  scanning  Earth  sensor  and  a  two-channel 
sun  sensor  housed  in  a  single  aluminum  enclosure.  Mounting  is  achieved  by  an 
asymmetric  three  point  mount  which  allows  for  easy  alignment  and  shimming  if  required. 
[Ref.  32] 

The  Earth  sensor  operates  over  an  attitude  range  of  ±  24°  in  pitch  and  ± 
14°  in  roll,  while  the  Sun  Sensor  typically  operates  with  sun  locations  of  ±  60°  in  azimuth 
and  ±  30°  in  elevation,  although  fields  as  large  as  128°  x  128°  may  be  accommodated.  All 
processing  and  calibration  data  utilized  to  derive  pitch  and  roll  output  from  the  Earth 
Sensor,  as  well  as  sun  azimuth  and  elevation  from  the  Sun  Sensor  is  available  on-board  the 
sensor.  Output  occurs  over  a  digital  interface  governed  by  the  MIL-STD-1553  data  bus, 
EIA  RS-422A,  or  user  specified  interface.  Upon  asserting  a  data  request  the  ES/SSA 
outputs  all  Earth  Sensor  and  Sun  Sensor  information  including  diagnostics  and  status 
information.  Table  5-2  lists  the  Earth  Sensor  /  Sun  Sensor  Assembly  Specifications.  [Ref. 
32] 
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EARTH  SENSOR 

Roll/Pitch  Accuracy  (Includes  all  long 
term  &  diurnally  varying  errors) 

±  0.02°  (3a) 

Transfer  Orbit  Roll  /  Pitch  Accuracy 

±0.1°  (3a) 

Attitude  Range 

±14°  Roll,  ±24°  Pitch 

SUN  SENSOR 

Range 

Up  to  128°  x  128° 

Accuracy 

±0.01°  (3a) 

OVERALL  INSTRUMENT 

Manufacturer 

ITHACO  Space  Systems 

Power  Consumption 

<  5.5  Watts 

Mass 

<  5.5  lb  (2.5  kg) 

Dimensions 

6.4  x  3.9  x  8.0  in 

Output  Data  Rate 

4  Hz 

Data  Interface  Options 

MIL-STD-1553,  EIARS- 

422A,  or  per  User  Spec 

Qualification  Life  (Including  Particle 
Radiation) 

22.5  yrs 

Qualification  Temperature 

-34°  C  to +71° 

Vibration 

22.5  grms  (Random) 

Parts  Program 

S  Level  and  ML-STD-1 547 

Reliability 

<163xl0"9  failures  /  hr 

Mounting 

Asymmetric  3  Point  Flange 
Mount 

Table  5-2:  The  Earth  Sensor  /  Sun  Sensor  Assembly  Specifications 
4.  The  Reaction  Wheel  Assembly 

a.  Application 

The  Type  E  Reaction  Wheel  Assembly  provides  a  reliable  source  of 
reaction  torque  and  angular  momentum  storage  for  attitude  control  of  medium  to  large 
spacecraft.  Accelerating  or  decelerating  a  flywheel  with  an  integral  motor  provides  a 
means  of  controlled  momentum  exchange  with  the  spacecraft,  which  is  advantageous  for  a 
variety  of  attitude  control  schemes.  [Ref.  33] 
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Featuring  high,  bi-directional  torque  capability,  a  set  of  Type  E  Reaction 
Wheels  is  ideal  for  agile  three-axis  reaction  control  applications.  With  a  momentum 
storage  capacity  of  50  N-m-s  and  an  available  reaction  torque  of  up  to  0.3  N-m,  the  Type 
E  Reaction  Wheel  can  satisfy  a  broad  range  of  spacecraft  requirements.  [Ref.  33] 

b.  Design  Description 

1.  Structure 

The  Type  E  Reaction  Wheel  consists  of  a  precision  balanced  solid  flywheel 
suspended  on  ball  bearings  and  driven  by  a  brushless  DC  motor.  The  aluminum  housing  is 
vented  to  space  to  eliminate  viscous  drag  torque  from  the  flywheel  due  to  windage. 
Removable  domed  covers  provide  access  to  the  wheel  interior  for  inspection  with  the 
flywheel  installed.  [Ref.  33] 

2.  Suspension 

The  flywheel  is  suspended  on  high  capacity  angular  contact  ball  bearings. 
A  low  vapor  pressure  lubricant  synthesized  specifically  for  high  vacuum  spaceflight 
applications  is  used  to  lubricate  the  bearings,  which  incorporates  a  flight  proven  extreme 
pressure  additive  to  prevent  bearing  wear  while  operating  at  low  speeds  below  50  rpm 
when  an  EHD  film  is  insufficient  to  prevent  ball  to  race  contact.  Labyrinth  seals  and 
barrier  film  coatings  ensure  long  life  operation  by  preventing  loss  of  lubricant  from  the 
bearings  due  to  evaporation  and  surface  migration.  A  life  test  of  the  Reaction  Wheel  in  a 
bi-directional  reaction  wheel  application  frequently  crossing  through  zero  speed  has  been 
operating  successfully  since  April  1994.  In  addition,  ITHACO  is  offering  an  On-orbit 
grease  replenishment  design  for  all  their  Reaction  Wheels  in  order  to  reliably  increase 
minimum  ball  bearing  operating  lifetime  from  eight  to  fifteen  years.  The  newly  developed 
device  utilizes  a  Kapton  heater  and  a  simple  thermostat  circuit  to  drive  a  sealed  grease 
reservoir  to  a  specific  temperature.  The  high  thermal  expansion  coefficient  of  the  grease 
results  in  the  injection  of  a  controlled  amount  of  grease  from  the  reservoir  directly  into  the 
bearing.  The  system  is  unpressurized  and  contains  no  active  moving  parts,  resulting  in 
reliable,  fail-safe  performance.  [Ref.  33,  Ref.  34] 

3.  Motor 

A  discretely  commutated,  ironless  armature  brushless  DC  torque  motor 
drives  the  flywheel  in  the  Reaction  Wheel.  Maximum  inertia  to  weight  ratio  is  realized 
with  the  large  diameter  motor  components  and  the  unique  feature  of  the  ironless  armature 
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design  which  strategically  places  all  of  the  mass  of  the  iron  and  magnets  on  the  rotor.  The 
ironless  armature  design  also  eliminates  cogging  for  smooth  transitions  during  direction 
reversals.  Hall  sensors  integrated  directly  into  the  armature  sense  the  alternating  flux  field 
from  the  permanent  magnets  to  provide  commutation  information  for  the  motor  driver,  as 
well  as  a  high  resolution  tachometer  signal  for  speed  and  direction  of  rotation  information. 
[Ref.  33] 

4.  Motor  Driver 

The  motor  driver  responds  proportional  to  an  analog  torque  command  and 
provides  linear  motor  torque  to  the  flywheel  over  the  entire  speed  range.  A  digital 
interface  per  MIL-STD-1553  or  other  custom  interface  can  be  easily  incorporated  in  the 
existing  board  space.  Analog  and  digital  tachometers  are  provided  to  accommodate  a 
variety  of  user  needs.  Due  to  the  high  torque  and  current  in  the  Reaction  Wheel,  the 
motor  driver  is  packaged  in  a  separate  enclosed  electronics  box  for  flexible  thermal 
management.  The  motor  drivers  are  designed  to  minimize  electrical  interference,  and  have 
demonstrated  superb  electromagnetic  compatibility  by  meeting  the  requirements  of  MIL- 
STD-461  and  MIL-STD-462.  [Ref.  33] 

5.  Specifications 

The  performance  specifications  for  the  Reaction  Wheel  (E  -  Wheel)  are 
presented  here  in  tabular  and  graphical  form.  Higher  angular  momentum  capacities  and 
higher  torque  versions  can  be  made  available  as  non-standard  designs.  Within  the  present 
Reaction  Wheel  case  sizes,  higher  reaction  torque  can  be  accommodated  with  a  rescaled 
motor  driver.  [Ref.  33] 


78 


ITHACO  "E  WHEEL"  REACTION  WHEEL  ASSEMBLY 

Speed  Range  (rpm) 

0  to  ± 3850 

Angular  Momentum  Capability 

50  N-m-s  (36.8  ft-lb-sec) 

Reaction  Torque 

300  mN-m  (42  oz-in) 

Weight  of  Wheel  Assembly 

10.5  kg  (23.1  lb) 

Weight  of  Driver  Assembly 

3.63  kg  (8.0  lb) 

Dimensions  of  Wheel  Assembly 

39.3  dia  x  16.6  cm  (15.5  x  6.6  in) 

Dimensions  of  Driver  Assembly 

17.8  x  17.8  x  10.2  cm  (7.0  x  7.0  x  4.0  in) 

Tachometer  (Pulses  per  revolution) 

72 

Static  Imbalance 

3.6  gram-cm  (0.05  oz-in) 

Dynamic  Imbalance 

1 20  gram-cm 2  (oz-in 2 ) 

Estimated  Lifetime 

1 5  years 

Table  5-3:  Ithaco  Reaction  Wheel  Specifications 
5.  The  Satellite  Inertial  Reference  System 

a.  Heritage 

Honeywell  has  been  involved  in  many  successful  space  programs  including 
launch  vehicles,  satellites,  space  probes  and  manned  space  flight.  Recent  experiences  have 
included  the  TOS  ,  Centaur  launch  systems,  MSX  and  Space  Station  Freedom. 

The  Miniature  Inertial  Measurement  Unit  (MIMU)  product  is  an 
integration  of  components  from  several  Honeywell  programs  into  a  standard  product. 
These  other  programs  have  provided  the  hardware  and  software  knowledge  and 
experience  to  develop  a  family  of  systems  consisting  of  three  axis  rotation  sensors  and 
optional  three  axis  linear  acceleration  sensors  in  a  small,  lightweight  radiation  hardened 
product.  The  MIMU  uses  the  Honeywell  GG  1320  ring  laser  gyroscope  which  is  in  high 
volume  production  the  their  Minneapolis  ring  laser  gyro  production  facility. 

b.  Miniature  Inertial  Measurement  Unit 

The  MIMU  family  was  designed  to  be  a  low  cost,  low  risk  radiation 
hardened  product,  to  be  used  as  the  inertial  measurement  system  in  a  broad  range  of 
satellite  applications.  The  system  contains  three  GG  1320  ring  laser  gyros  mounted 
orthogonally  as  the  rotation  sensors  and  three  optional  Allied  Signal  QA  3000 
accelerometers  as  the  linear  acceleration  sensors.  A  block  diagram  of  the  MIMU  is  shown 
in  Figure  5-5.  [Ref.  35] 
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Figure  5-5:  The  Miniature  Inertial  Measurement  Unit 


80 


The  Honeywell  MIMU  Model  YG9666  B 

Size 

250  in3 

Weight  (w/o  magnetic  shield) 

7.61b 

Power 

33.7  Watts 

Rate  Capability  (deg/sec) 

±375 

Acceleration 

25  g 

Life  (yrs,  operation  dependent) 

15 

Gyro-GG1320 

Bias  Stability 

(la,  deg/hr;  short  term) 

0.008 

Angular  Random  Walk 

(deg/  4hr) 

0.005 

Scale  Factor  Error 
(la,  ppm;  short  term) 

5 

Readout  Noise 
(urad,  la) 

5 

Accelerometer 

QA  3000 

Bias  Stability 

(la,  mg;  short  term) 

0.04 

Scale  Factor  Error 
( 1  a,  ppm;  short  term) 

80 

Output  Data 

Angle  LSB  (uxad) 

1.0 

Velocity  LSB  (mm/sec) 

0.0753 

Temperature 

Operating  (°C) 

-30  to  +70 

Non-Operating  (°C) 

-45  to  +75 

Radiation  (krads  total  dose) 

100 

Table  5-4:  The  MIMU  Performance  Parameters 

The  MIMU  offers  excellent  long  term  stability,  very  high  reliability  and  a 
wide  dynamic  range.  Recent  reliability  studies  of  field  life  data  has  shown  the  life  of  the 
gyroscope  to  be  in  excess  of  300,000  hours.  This  gives  the  MIMU  a  high  probability  of 
success  for  missions  of  15  years  or  more.  [Ref.  35] 
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The  processor  electronics  assembly  consists  of  the  Honeywell  radiation 
hardened  GVSC,  RH-1750  processor,  a  16  bit  processor  utilizing  the  full  MIL-STD- 
1 750A  instruction  set.  The  software  performs  the  functional  control  as  shown  in  the  block 
diagram  of  Figure  5-6.  This  assembly  also  contains  a  standard  RS422  interface  operating 
at  a  200  Hz  data  rate.  MIMU  can  also  be  configured  to  use  an  optional  MIL-STD-1553B 
interface.  [Ref.  35] 

The  standard  RS422  interface  consists  of  ten  1 6-bit  data  words.  These 
data  words  are  shown  in  Table  V-5. 


Wordl 

Frame  Timer  (MSH) 

Word  2 

Frame  Timer  (LSH) 

Word  3 

X  Compensated  Gyro  Angle 

Word  4 

Y  Compensated  Gyro 

Word  5 

Z  Compensated  Gyro  Angle 

Word  6 

X  Compensated  Velocity 

Word  7 

Y  Compensated  Velocity 

Word  8 

Z  Compensated  Velocity 

Word  9 

Multiplexed  BIT  Status 

Word  10 

Checksum 

Table  5-5:  MIMU  Output  Data  Message 

The  data  rate  is  1  Mbit/sec.  The  multiplexed  BIT  Status  word  consists  of 
256  words  multiplexed  at  one  word  output  each  200  Hz  period.  The  BIT  status  word  is 
identified  by  the  LSH  of  the  frame  timer  word.  [Ref.  35] 

The  accelerometer  electronics  are  also  located  on  the  processor  assembly. 
The  QA3000  Q-Flex  accelerometer  is  a  torque-rebalance  sensor  with  fused  quartz  flexures 
with  a  permanent  magnetic  torquer,  a  capacitive  pickoff  system,  and  self-contained  servo 
electronics.  The  accelerometer  outputs  a  current  proportional  to  the  acceleration.  The 
accelerometer  electronics  on  the  processor  assembly  digitizes  this  current  (three  channels) 
which  is  then  read  by  the  processor.  Mathematical  modeling  of  bias,  scale  factor,  and  axis 
misalignment  is  performed  over  the  temperature  range  to  enhance  the  performance.  [Ref. 
35] 

The  gyro  electronics  assembly  contains  the  electronics  for  control  of  the 
three  Ring  Laser  Gyros.  The  RLG  readout  post  amplifier  provides  the  gain  to  amplify  the 
gyroscope's  photodetector  outputs  to  levels  which  can  be  digitized.   These  readout  pulses 
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are  then  accumulated  and  are  read  by  the  processor.  Pulses  are  accumulated  for  positive 
angular  rotation  and  for  negative  angular  rotation.  Net  rotation  is  then  calculated  by  the 
processor.  [Ref.  35] 

To  over  come  a  lock-in  mechanism  at  low  rates,  the  gyro  is  dithered,  which 
involves  subjecting  the  gyroscope  to  a  small  angular  oscillation  through  a  controlled 
mechanical  resonance.  The  M1MU  gyros  are  oscillated  at  three  separate  frequencies 
(nominally  525,  575,  and  625  Hz).  The  baseline  means  of  removing  the  dither  oscillation 
from  the  gyro  signal  is  through  a  procedure  called  "dither  stripping".  The  procedure 
involves  measuring  the  angle  between  the  gyro  block  and  the  gyro  mounting  base  and 
subtracting  this  angle  from  the  gyro  output  signal.  The  output  of  a  piezoelectric  sensor 
attached  to  the  dither  spring  is  passed  through  an  A/D  converter  to  the  digital  processor, 
scaled  and  corrected  for  phase,  and  subtracted  from  the  gyro  signal.  The  result  is  a  wide 
bandwidth  signal  with  the  content  at  the  dither  frequency  removed.  An  alternate  means  of 
removing  dither  is  filtering.  Filtering  is  used  when  low  readout  noise  is  required,  and  the 
corresponding  latency  penalty  can  be  tolerated.  The  dither  control  electronics  are  also 
located  on  the  gyro  electronics  assembly.  [Ref.  35] 

The  gyro  electronics  assembly  also  contains  the  path  length  control  (PLC) 
circuitry.  In  order  to  maintain  the  RLG  laser  at  its  peak  operating  point  (peak  power), 
PLC  is  utilized  by  using  a  PZT  which  converts  a  voltage  input  to  a  path  length  distance. 
The  lasing  power  is  at  maximum  when  the  laser  beam  path  length  is  equal  to  an  integral 
number  of  wavelengths.  [Ref.  35] 

The  gyro  electronics  assembly  also  contains  the  RLG  current  control 
circuitry.  The  current  in  each  leg  of  the  RLG  is  precisely  controlled  to  minimize 
temperature  and  gas  flow  induced  bias  changes.  [Ref.  35] 

The  Low  Voltage  Power  Supply  modules  convert  the  spacecraft  +28  Vdc 
to  the  isolated  secondary  voltages  required  for  operations  (+5  Vdc,  +15  Vdc,  and  -15 
Vdc).  [Ref.  35] 

The  High  Voltage  Power  Supply  (HVPS)  module  provides  +2200  Vdc  to 
the  gyroscope  cathodes  during  startup.  After  startup  the  HVPS  provides  -750  Vdc  to 
maintain  the  laser  currents.  The  PLC  circuitry  requires  +250  Vdc  provided  the  HVPS. 
[Ref.  35] 

6.  HD-1003  Star  Tracker 

The  HD-1003  Star  Tracker,  produced  by  Hughes  Danbury  Optical  Systems,  Inc.  is 
designed  to  meet  current  and  future  needs  for  a  lightweight,  reliable  and  accurate  star 
tracker.   A  simple  optical  system  ASIC  based  electronics  and  high  throughput  firmware 
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architecture  make  the  HD-1003  mission  adaptable,  even  while  in  operational  orbit.  [Ref. 
36] 

The  HD-1003  attitude  determination  system  provides  continuous  three-axis 
attitude  data  over  the  entire  orbit.  An  autonomous  all-stellar  attitude  determination 
system  is  available,  minimizing  spacecraft  computation  and  ground  support.  Since  the 
HD-1003  has  no  moving  parts,  high  reliability  and  long  life  are  achieved.  Further,  since 
the  HD-1003  system  can  provide  multiple  star  updates  every  one  tenth  of  a  second  for 
LOS  rates  up  to  1.0  deg/sec,  reliance  on  the  inertial  reference  unit  can  be  reduced  or 
eliminated.  [Ref.  36] 

Some  key  features  of  the  HD-1003  include: 

A.  Performance 

1 .  Six-star  tracking  capability 

2.  Six  arc-sec  rms  per  star 

3.  Two  arc-sec  rms  multi-star  average 

4.  10  Hz  update  rate 

5.  8  x  8  9  or  20  degree  circular  FOV 

6.  Parts  count  has  been  reduced  by  5  times. 

B.  Growth 

1 .  12  MIPS  throughput  with  50%  margin 

2.  Expandable  memory 

3.  Uplinkable  code 

4.  Small  or  large  space-craft  compatible 

C.  Survivability 

1 .  All  electronic  parts  are  hardened  to  100  Krad  or  higher 

2.  Built-in  radiation  shielding 

3.  Space  qualified  proton/debris  rejection  software 

4.  Expected  life  on  orbit  is  15  to  18  years 

5.  HDOS  has  completed  a  2  year  test  program  to  verity  the  HD-1003's 
performance  in  high  radiation  environments. 

6.  Test  results  and  design  have  been  independently  verified  by  Lockheed 
Martin,  Aerospace  and  the  Naval  Research  Lab. 
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D.  Producibility 

1 .  Modular  design 

2.  ASIC-based  electronics 

3.  No  moving  parts,  sun  safe 

4.  Safety  engineered 

The  HD-1003  Star  Tracker  benefits  from  over  Twenty  years  of  experience  in 
leading  attitude  determination  technology.  Hughes  Danbury  produced  the  world's  first 
CCD-based  star  sensor  to  be  qualified  in  space,  the  first  "multi-star"  tracker,  and  the  first 
CCD  upgrade  to  NASA's  Standard  Fixed  Head  Star  Tracker.  The  Hughes  Danbury  team 
continues  to  satisfy  evolving  demands  by  utilizing  the  latest  technology  and  innovative 
packaging.  [Ref.  36] 
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Figure  5-6:  CCD  Focal  Plane  and  Microprocessor  Controller  Sample  and  Centroid  up  to 

Six  Star  Images  Simultaneously. 
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Performance  Category 

HD-1003  Narrow 
FOV 

HD-1003  Wide  FOV 

Size  (with  light  shade)  (L  x  W  x  H) 

16x6.2x4.4  in. 

7x6.2x4.4  in 

Weight  (with  light  shade) 

8.0  1b 

7.1  lb 

Power  (average  at  28  Vdc) 

9  Watts 

9  Watts 

Communication  Interface 

MIL-STD-1553B 

MIL-STD-1553B 

Field  of  View  (deg) 

8°x8° 

20°  circular 

Sensitivity  (5000  K  star)  (  mv ) 

+6.0 

+4.7 

Overall  Accuracy,  Each  Axis 

*  Pitch  /  Yaw  (per  star)  (arc-sec, 

6 

30 

rms) 

2 

10 

*  Multi-star  average  (arc-sec,  rms) 

40 

50 

*  Roll  (5-stars)  (arc-sec,  rms) 

0.25 

0.25 

*  Magnitude  (±) 

Update  (frame)  Rate  (Hz) 

10 

4 

Stars  Simultaneously  Tracked 

6 

6 

Bright  Object  Rejection  Angle** 

*  Sun  (deg) 

35° 

40° 

*  Earth  (deg) 

25° 

30° 

*  Moon  (deg) 

25° 

25° 

Acquisition  Time  (6  stars)  (sec.  1-sigma) 

6 

6 

Mean  Time  Between  Failures  (Hours) 

1,000,000  hr 

1,000,000  hr 

Environments: 

Temperature 

*     Nominal  (  °C) 

-5  to  +45 

-5  to  +45 

*     Survival  (  °C) 

-15  to +55 

-15  to +55 

Vibration  (Random)  (g  rms) 

12.2 

12.2 

EMC  Level 

MIL-STD-461C 

MIL-STD-461C 

**  Light  shade  is  modular  to  accommodati 

i  various  mission  requirem 

ents 

Table  5-6:  HD-1003  Star  Tracker  Performance  Capability 
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7.  Global  Positioning  System 

A  Global  Positioning  System  receiver  is  required  to  provide  the  spacecraft  with  an 
autonomous  capability  to  accurately  determine  its  position.  The  Motorola  GPS 
Spaceborne  Receiver  system  was  selected  for  this  mission.  This  system  is  capable  of 
determining  the  spacecraft's  position  to  within  1 5  meters  for  Autonomous  Position  and  2 
to  5  meters  for  Differential  Position;  provided  that  the  Selective  Availability  is  not 
implemented.  This  system  also  provides  the  ability  to  determine  spacecraft  attitude  to 
accuracies  of  0. 1  °  in  pitch,  roll  and  yaw.  In  addition,  the  system  can  determine  pitch,  roll 
and  yaw  rates  to  accuracies  of  0.5°/second.  These  accuracies  are  not  sufficient  to  meet 
the  crosslink  pointing  knowledge  requirements  but  are  suitable  to  provide  a  backup  source 
of  attitude  information  for  basic  spacecraft  control  functions  in  the  event  of  failures  in  the 
primary.  Table  5-7  identifies  the  characteristics  of  the  Motorola  GPS  Spaceborne 
Receiver. 

Additional  key  features  of  the  Motorola  GPS  Spaceborne  Receiver  include: 

1 .  Space  qualified  digital  design. 

2.  Full  spaceborne  doppler  capability. 

3.  Autonomous  operation  with  doppler  capability  at  1  second  rate. 

4.  Non-destruct  doppler  measurements  at  1  second  rate. 

5.  Radiation  hardened  1  PPS  clock  output  synchronized  to  GPS  time. 
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Multiple  Configurations  Available 

*Ll,C/AandP-code 
♦L1/L2,  C/AandP-code 

*  Y-code-  Anti-Spoof  and  Selective  Availability 

*  Multiple  Antennas 

*  6  to  1 8  channels 

Physical  /  Environment 

Size 

6.0"  x  5.2"  x  2.7" 

Weight 

3.71b 

DC  Power 

7  Watts 

Vibration 

25  grms 

Shock 

3000  gat  10  kHz 

Temperature 

-20°  C  to  +75°  C 

Performance  Characteristics 

Receiver  Accuracies  (over  C/N  o  Range  of  32  to  42  dB/Hz): 

Pseudo-Range  (10  sec) 

0.2  m 

Pseudo-Range  Interchannel  Bias 

<0.1  m 

Carrier  Phase  (1  sec) 

0.01m 

Navigation  Solution  Accuracy  (GDOP  =  3,  Selective  Availability  not 
implemented): 

Autonomous  Position 

<15m(lc) 

Differential  Position 

2  to  5  m  (la) 

Time  Offset 

100  ns 

Time-to-First-Fix 

Mean  TTFF 

<  60  sec 

Maximum  TTFF  (P  f  >  0.9) 

<  90  sec 

User  S/C  Ephemeris 

<  10  km,  lOm/s 

User  Clock  Error 

<1  xlO-7 

Orbital  Dynamics  (typical  /  spacecraft  orbit): 

Velocity 

8,000  m/s 

Acceleration 

10m/s2 

Table  5-7:  Characteristics  of  the  Motorola  GPS  Spaceborne  Receiver 
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8.  Sun  Sensor 

a.  Overview 

Two  Sun  sensors  will  be  used;  one  on  each  solar  array  substrate  assembly. 
These  two  sun  sensors  will  be  the  primary  method  for  determining  the  yaw  angle  and  solar 
array  rotation  angle  required  to  maintain  solar  tracking.  The  Sun  Sensor  Assembly 
selected  is  a  two  axis  digital  sun  sensor  which  offers  unprecedented  performance  in  terms 
of  accuracy,  mass,  power  consumption,  reliability  and  automation.  The  sensor 
autonomously  determines  two-axis  attitude  information  with  respect  to  the  sun  line  from 
low  Earth  through  super  geosynchronous  orbits. 

b.  Theory  of  Operation 

The  Sun  Sensor  Assembly  incorporates  two  static,  mutually  orthogonal 
sensing  (axes)  arrays  to  detect  the  orientation  of  the  sun.  Each  axis  is  comprised  of  a  slit 
allowing  a  plane  of  light  to  pass  through  it  intersecting  with  an  orthogonally  mounted, 
high  resolution,  linear  photodiode  array.  The  point  of  intersection,  which  indicates  sun 
angle,  is  obtained  by  reading  the  energy  stored  in  the  array. 

The  Sun  Sensor  Assembly  typically  operates  with  sun  locations  of  ±  64°  in 
azimuth  and  ±  64°  in  elevation.  All  processing  and  calibration  data  utilized  to  derive  sun 
azimuth  and  elevation  from  the  Sun  Sensor  Assembly  occurs  on  board  the  sensor.  Output 
occurs  over  a  digital  interface  governed  by  the  MIL-STD-1553  data  bus,  EIA  RS-422A, 
or  user  specified  interface.  Upon  asserting  a  data  request  the  SSA  outputs  all  SSA 
information  including  diagnostics  and  status  information.  Units  have  been  provided  for 
the  MSTI,  APEX,  and  Sea  Star  Programs. 
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VI.  YAW  CONTROLLING 

On  orbit  pointing  control  utilizes  three  reaction  wheels  to  provide  roll,  pitch  and 
yaw  attitude  control.  An  earth  sensor  provides  two  axis  attitude  data  for  roll  and  pitch 
and  a  sun  sensor  provides  one  axis  attitude  data  for  yaw.  The  yaw  control  is  similar  to 
both  roll  and  pitch  except  that  it  must  utilize  additional  information  in  order  to  maintain 
solar  array  sun  tracking.  Also,  changes  is  yaw  do  not  adversely  effect  the  communications 
payload  as  long  as  pointing  error  limits  about  the  roll  and  pitch  axis  are  maintained.  In  this 
section  the  peculiarities  of  yaw  steering  will  be  discussed. 

A.  YAW  DISTURBANCES 

Due  to  secular  disturbances,  momentum  accumulates  in  the  wheel,  and  wheel 
speeds  will  build  up  with  time.  In  normal  operations  these  secular  torques  in  the  yaw 
direction  will  be  compensated  by  operation  of  the  magnetic  torque  rods  In  the  back-up 
mode  momentum  would  be  dumped  through  a  series  of  minimum  impulse  thruster  pulses. 

B.  SUN  AND  NADIR  POINTING  GEOMETRY 

Figure  6-1  illustrates  the  Sun  and  Nadir  pointing  geometry  for  the  MUS 
spacecraft.  P  represents  the  azimuth  and  8  the  elevation  of  the  sun  line  relative  to  the 
spacecraft  orbit  plane.  P  is  referenced  to  the  X  axis.  The  X  -  axis  coincides  with  the 
velocity  vector  in  the  local  geodetic  reference  frame  (i.e.  orbit  nominal  frame).  The 
commanded  yaw  and  array  angles  vary  as  a  function  of  P  and  8  due  to  the  effects  of  the 
space  vehicles  orbital  motion,  precession  of  the  orbit  plane,  and  movement  of  the  sun  in 
the  earth  centered  inertial  frame.  These  commands  maintain  the  normal  to  the  solar  panels 
directed  toward  the  sun  and  the  sun  in  the  spacecraft  roll/yaw  plane.  Except  for  very  low 
sun  angles,  when  8  <  5.0°,  the  -X  side  of  the  spacecraft  is  always  oriented  towards  the 
sun.  [Ref.  27:  p.  458] 
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Figure  6-1 :  Pointing  Geometry 

C.  SUN  AND  NADIR  YAW  POINTING  STRATEGY 

The  Attitude  Dynamics  and  Control  Subsystem  is  designed  to  track  the  nominal 
sun/nadir  yaw  trajectory  for  sun  angles  (8)  greater  than  5.0°  from  the  orbit  plane.  The 
closed  loop  controller  uses  sun  sensor  data  as  the  yaw  attitude  reference.  When  the 
spacecraft  is  in  the  sun/earth  co-linearity  region,  defined  when  8  is  less  than  15°  and  the 
projection  of  the  sun  line  in  the  orbit  plane  is  within  15°  of  the  nadir  vector,  yaw  sensing  is 
restricted  due  to  either  poor  geometry  or  an  eclipse  condition.  In  these  regions,  a  yaw 
momentum  tracking  tach  loop  is  used  to  track  an  ideal  yaw  momentum  profile.  Predicted 
Euler  coupling  torques  are  fed  forward  to  achieve  a  high  degree  of  accuracy,  minimizing 
yaw  pointing  error  at  exit  from  the  region.  [Ref.  27:  p.  460] 
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When  the  sun  is  within  5.0°,  nominal,  of  the  orbit  plane,  the  ADCS  enters  a 
constant  yaw  mode  which  maintains  the  pitch  (array  rotational)  axis  on  the  orbit  normal 
and  an  orbit  rate  rotation  about  the  pitch  axis.  The  constant  yaw  mode  minimizes 
Reaction  Wheel  torque  and  momentum  demand  while  maintaining  specified  yaw  pointing. 
Table  6-1  summarizes  the  sun  and  nadir  yaw  pointing  strategy.  [Ref.  27:  p.  460] 


Yaw 

Sun  Rotation 

Sun  to 

Yaw  Control 

Remarks 

Tracking 

Angle  p(deg) 

Orbit  Plane 

Reference 

Mode 

Angle  (8) 

All  P  angles 

|8|>15° 

ST/Yaw  SS 

Output  Null 

Position 

Closed  loop  yaw  control  using 

Star  Trackers  (ST)  /  SS  derived 

yaw  &  yaw  rate 

Ideal 

|P-90°|>15° 

15°>|5|> 

ST/Yaw  SS 

Closed  loop  yaw  control  using 

Yaw 

and 

5.0° 

output  null 

ST/SS  derived  yaw  &  yaw  rate 

Tracking 

|p-270°|>15° 

position 

Mode 

|p  -  90°|  <  15° 

15°>|S|> 

ST/Calculated 

Poor  yaw  signal  from  SS  when 

or 

5.0° 

yaw  cmd,  no 

|5|  <  1 5°  due  to  array  position  or 

IP  -270°|  <  15° 

SS  yaw 
sensing 

no  SS  output  due  to  eclipse 

|p-  90°|  >  15° 

|5|  <  5.0° 

ST/Yaw  SS 

Constant  pitch  rate  control; 

and 

output  biased 

Closed  loop  yaw  control  using 

|P-270°|>15° 

at  8 

ST/SS  derived  yaw  and  yaw  rate. 

Constant 

yaw 

|p-90°|<15° 

|5|  <  5.0° 

ST/Zero  yaw 

Constant  pitch  rate  control;  Feed 

tracking 

or 

angle;  no  SS 

forward  tach  loop  yaw  rate 

mode 

IP-  270°|  <  15° 

yaw  sensing 

tracking  control;  Poor  SS  yaw 
signal  when  |P|  <  15°  due  to 
array  position  or  no  SS  output 
due  to  eclipse;  Minimum  yaw 
control  torque  required. 

Table  6-1 :  Sun  and  Nadir  Yaw  Pointing  Strategy 


95 


D.  SUN  AND  NADIR  ARRAY  CONTROL  STRATEGY 

The  ADCS  is  designed  to  track  the  nominal  sun/nadir  array  trajectory  for  sun 
angles  (8)  greater  than  5.0°  from  the  orbit  plane.  The  closed  loop  controller  uses  course 
sun  sensor  data  as  the  array  loop  null  reference.  When  the  spacecraft  is  in  eclipse,  a 
calculated  nominal  array  trajectory  drives  the  array  loop  with  position  feedback  provided 
by  the  measured  array  angle.  When  the  sun  is  within  5.0°,  nominal,  of  the  orbit  plane,  the 
ADCS  enters  a  constant  yaw  mode.  The  array  control  strategy  is  unchanged,  but  the 
constant  yaw  geometry  results  in  an  orbital  rate  rotation  of  the  array.  The  transition  into 
and  out  of  the  constant  yaw  mode  is  limited  to  the  segment  of  the  orbit  with  sun  rotation 
angles  between  120°  and  240°;  i.e.  between  2pm  and  10pm  spacecraft  local  time, 
precluding  the  need  for  yaw  inversion  maneuvers.  Table  6-2  describes  the  sun/nadir  array 
control  strategy.  [Ref.  27:  p.  462] 


Solar  Array 

Eclipse 

Sun  to  Orbit 

Array  Control 

Remarks 

Tracking 

Condition 

Plane  angle  8 

Reference 

Mode 

No 

|8|  >  5.0° 

ST/Pitch  SS  output 

No  eclipse  will  occur 

Ideal  Sun- 

eclipse 

null  position 

when|S|>  15° 

Nadir 

Array 

The  error  between 

Tracking 

Eclipse 

15°>|8|>5.0° 

Calculated  ideal 
array  angle  command 

measured  array 

position  and 

command  reference 

is  used  for  control 

No 

|S|  <  5.0° 

ST/Pitch  SS  output 

Array  sun  pointing 

Constant 

eclipse 

null  position 

errorS 

Array 

Rate 

|8|  <  5.0° 

Calculated  array 

Measured  array 

Tracking 

Eclipse 

angle  command  for 

constant  yaw 

tracking 

position  used  for 

control;  sun  pointing 

error  equals  8 

Table  6-2:  Sun  and  Nadir  Array  Control  Strategy 
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E.     NORMAL,  ON-STATION  CONTROL  FUNCTIONAL  BLOCK 
DIAGRAM 

The  Attitude  Dynamics  and  Control  Subsystem  provides  fully  autonomous  three 
axis  attitude  determination  and  control  throughout  the  Space  Vehicle's  operational 
lifetime.  Three  of  four  redundant  reaction  wheels  control  body  attitude,  and  redundant 
magnetic  torquing  coils  control  system  momentum.  There  is  no  thruster  utilization  for  the 
Normal,  On-Station  mode  of  operation.  Thrusters  are  used  primarily  to  provide  A  V  and 
control  actuation  during  stationkeeping  maneuvers,  but  also  provide  backup  momentum 
unloading  in  the  event  of  a  failure  of  the  magnetic  system.  Also,  control  loop 
compensation  is  provided  in  the  attitude  control  software  for  Euler  coupling  torques, 
magnetic  and  thruster  desaturation  torques,  and  wheel  friction  torques  to  optimize 
pointing  and  reduce  jitter.  [Ref.  27:  p.  266] 

While  the  system  is  fully  autonomous,  the  Control  Segment  has  full  override 
access  to  the  ADCS  equipment  configuration  and  operating  mode.  The  internal  state  of 
all  software  flags  is  available  in  telemetry  and  any  control  mode  or  functions  may  be  either 
inhibited  or  forced  by  a  ground  override  of  the  appropriate  flag(s).  Any  control  algorithm 
may  also  be  reprogrammed  by  ground  command.  The  system  is  capable  of  operation  for 
periods  of  at  least  210  days  with  no  ground  intervention.  Figure  6-2  illustrates  the 
Normal,  On-Station  Functional  Block  Diagram.  [Ref.  27:  p.  466] 
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Figure  6-2:  Normal,  On-Station  Functional  Block  Diagram 

F.       ATTITUDE    DETERMINATION    SOFTWARE    FUNCTIONAL 
OVERVIEW 

The  attitude  determination  software  comprises  two  primary  functions;  namely, 
attitude  sensor  data  processing  and  payload  ephemeris  data  processing.  The  attitude 
sensor  processing  conditions  the  raw  data  input  from  the  various  attitude  sensors, 
processes  the  data  to  derive  three-axis  spacecraft  attitude  and  attitude  rate,  and  generates 
the  yaw  and  solar  array  drive  reference  commands  to  drive  the  respective  control  loops  in 
the  co-linearity  regions  and  for  low  sun  angles  (8<5.0°).  The  ephemeris  processing, 
which  uses  Navigation  Payload  data  as  the  primary  reference  and  ground  uplink  data  as 
the  backup  reference,  generates  sun  angle  data  for  sun/nadir  pointing  mode  operation, 
sun/moon  interference  predictions  for  earth  sensor  processing,  and  a  magnetic  field  model 
for  the  momentum  unloading  software.  Figure  6-3  is  a  Functional  Overview  of  the 
Attitude  Determination  Software.  [Ref.  27:  p.  468] 
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G.    SUN  ANGLE  ESTIMATION  ((3  AND  5) 

The  ADCS  processes  both  Navigation  payload  data  as  well  as  ground  uploaded 
ephemeris  to  provide  a  reliable  basis  for  sun/nadir  pointing.  The  spacecraft  ephemeris 
obtained  from  the  on-board  navigation  payload  and  the  sun  ephemeris  generated  on-board 
are  used  to  compute  the  sun  rotation  (azimuth)  angle  (P)  and  the  out-of  orbit  plane  sun 
angle  (8)  from  which  eclipses  times  and  yaw/array  control  phase  is  determined.  [Ref.  27: 
p.  470] 

The  backup  mode  determines  8  using  the  orbit  right  ascension  of  ascending  node, 
Q,  and  the  orbit  inclination,  i,  along  with  the  on-board  generated  sun  ephemeris.  P  is 
estimated  from  the  calculated  8  angle  and  the  array  position  when  the  array  position 
relative  to  the  yaw  axis  is  suitable  and  attitude  pointing  is  within  tolerance.  Yearly 
variations  of  Q  and  i  are  represented  by  an  on-board  model,  thus  each  upload  covers  a 
time  span  of  one  year.  To  satisfy  the  autonomy  requirements,  however,  uploads  are 
required  every  few  months.  The  backup  estimator  is  run  at  all  times,  but  the  estimated  P 
and  8  are  only  used  when  payload  data  is  unavailable  or  invalid.  Figure  6-4  is  the  flow 
diagram  for  the  estimation  of  Sun  angles  P  and  8.  [Ref.  27:  p.  470] 
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Figure  6-4:  Sun  Angle  Estimation 
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H.  YAW  DATA  PROCESSING 

The  ADCS  utilizes  a  two-axis  sun  sensor  for  yaw  attitude  control  and  solar  array 
drive  control.  The  Sun  sensor  is  comprised  of  two  sensing  units,  one  located  on  each 
solar  array,  which  provide  redundant  detector  sets  and  processing  electronics  for  each 
sensing  axis.  Data  from  the  two  detector  units  is  differenced  to  provide  a  measurement 
that  is  linear  on  each  axis  over  approximately  ±  40°  and  covers  a  full  range  of  ±  125°, 
minimum,  per  axis.  For  yaw  determination,  the  sensor  data  must  be  weighted  as  a 
function  of  the  solar  array  orientation.  [Ref.  27:  p.  472] 
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VII.  SIMULATION  OF  THE  ATTITUDE  DYNAMICS  AND 
CONTROL  SUBSYSTEM 

A.  PROBLEM  STATEMENT 

The  MUS  spacecraft  is  a  Medium  Earth  Orbit  satellite  which  has  the  following 
moments  of  inertia  of: 

Iroll  =  13,337.7     kg-m2 
/p/toA  =  14,083.6     kg-m2 
7^  =  14,113.7     kg-m2 

The  pointing  accuracy  for  the  spacecraft  is: 

a)  0.1°  for  Roll  and  Pitch 

b)  1.0°  for  Yaw 

The  disturbance  torques  on  the  spacecraft  are  primarily  due  to  solar  pressure,  and 
can  be  approximated  by: 

r,  =  rro//=(lxl0-3)*cos(a>o0  N-m 

Ty  =  Tpilch  =  (1.5  x  lO3)  *  cos(coj)         N  -  m 
T:  =  Tyaw=(\x\0-')*cos(o)ot)  N-m 

The  thruster  misalignment  torque  about  any  axis  can  be  approximated  by: 

7>  =  lxl(T4  N-m 

Finally,  the  design  of  the  attitude  control  system  for  the  spacecraft  must  meet  the 
given  specifications  when  the  environmental  torques  are  imposed  on  the  spacecraft. 
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B.    THREE  AXIS  STABILIZATION 

The  control  torques  along  the  axes  of  three-axis-stabilization  systems  are  provided 
by  various  combinations  of  momentum  wheels,  reaction  wheels,  and  thrusters.  Broadly, 
however,  there  are  two  types  of  three-axis-stabilization  systems:  a  momentum  biased 
system  with  a  momentum  wheel  along  the  pitch  axis,  and  a  zero-momentum  system  with  a 
reaction  wheel  along  each  axis.  [Ref.  25:  p.  130] 

In  momentum  wheel  systems,  the  angular  momentum  along  the  pitch  axis  provides 
gyroscopic  stiffness.  In  these  systems,  the  pitch  and  roll  axes  are  controlled  directly  and 
the  yaw  axis  is  controlled  indirectly  due  to  a  yaw  sensor.  In  reaction  wheel  systems,  all 
three  axes  are  controlled  independently,  thus  requiring  a  yaw  sensor  also.  The  momentum 
wheel  can  be  fixed  or  gimbaled  about  one  or  two  axes.  The  control  torque  along  the  pitch 
axis  is  provided  by  the  change  in  the  speed  of  the  momentum  wheel.  The  torque  along  the 
roll  axis  is  provided  by  thrusters,  a  reaction  wheel,  or  by  changing  the  gimbal  angles  in  the 
gimbaled  momentum  wheel.  [Ref.  25:  p.  130] 

1.  Small  Angle  Approximation 

When  the  satellite's  body  frame  orientation  remains  very  close  to  the  orbit  nominal 
frame,  a  small  angle  approximation  is  achievable  as  well  as  being  able  to  disregard  higher 
order  terms.  This  greatly  simplifies  the  problem.  A  linearized  analysis  of  the  three-axis- 
stabilization  system  is  given  in  this  section.  The  linearized  equations  of  motion  of  a  three- 
axis-stabilized  spacecraft  are  derived  first.  The  angular  velocities  of  the  spacecraft  can  be 
expressed  in  terms  of  the  orbital  rate,  co0  ,  and  the  attitude  error  angles,  i//*9,and  ^, 
which  are  known  as  yaw,  pitch,  and  roll  errors,  respectively,  as  follows: 
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(7-1) 


The  angular  momentum  of  the  system  can  be  written  in  the  form 


H  =  H  +  H„ 


(7-2) 


where, 


K=[k  k  hz]T 

Hh  =  [lxxo)x     Iyycoy     4»z] 


(7-3) 


are  the  angular  momentum  of  the  wheels  and  of  the  spacecraft,  respectively,  where 
I^^Iyy,  and  4.  are     the     principal     moments     of    inertia     of    the     spacecraft     and 

cox,cov,  and  at.  are  the  angular  velocity  components  of  axes  x,  y,  and  z  as  given  by  Eq. 

(7-1).       Introducing    Eq.    (7-2),    in    conjunction    with    Eqs.    (7-3)    and    (7-1)    into 

Mc  =  HCrei+  at  x  Hc  and  ignoring  nonlinear  terms,  we  obtain  the  moment  equation  in  the 
form 


M  = 


K  +  L(w  +  vj)  +  hy(if>  -  cow)  -  hx(e  -  co0)  -  [iyy  -  /„)fi)0(^  -  (oDw) 


(7-4) 
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2.  Large  Angle  Requirements 

Unfortunately,  the  yaw  steering  requirements  of  the  satellite  does  not  allow  us  to 
either  make  small  angle  approximations  or  to  disregard  higher  order  terms.  This  greatly 
complicates  the  equations  of  motion  which  are  presented  below.  As  previously  discussed, 
the  angular  velocities  of  the  spacecraft  can  be  expressed  in  terms  of  the  orbital  rate,  od0 , 
and  the  attitude  angles,  (j>,  0,  andy/,  which  are  known  as  roll,  pitch,  and  yaw  respectively, 
as  follows: 
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(7-5) 


The  angular  momentum  of  the  system  can  again  be  written  as  stated  in  Eq.  (7-2) 


H  =  K,  +  Hu 


(7-2) 


where, 


Hw=[hx     hy    h2]T 

Hb=\lxxO)x  IyyODy  I 'jD ,] 


(7-3) 
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are  the  angular  momentum  of  the  wheels  and  of  the  spacecraft,  respectively,  where 
1^,1^,,  and  I=     are     the     principal    moments     of    inertia    of    the     spacecraft     and 

cox,60y,  and  oo.  are  the  angular  velocity  components  of  axes  x,  y,  and  z  as  given  by  Eq. 

(7-5).       Introducing    Eq.    (7-2),    in    conjunction    with    Eqs.    (7-3)    and    (7-5)    into 

Mc  =  Hcrei+  oo  x  Hc  and  keeping  all  terms,  we  obtain  the  moment  equations  in  the  form: 


M  = 


hx  +  I^+Ocos^K  +  Ocosi^I^y/  +  y/cos(0)sm(0)h,  +  y/2  cos(#)sin(^)/_  - 
6>0cos(^)cos(y)/*r  -<y0cos(^)cos(^)/_^-ft>osin(^)sin(6>)sin(^)/?.  - 
oo0  sin(^)sin(<9)sin(^)4^  +  0sm(<p)hv  +  01  sin(<f>)I>y  -  y/cos(0)cos(0)hv  - 
y/ cos(0)cos(0)  I yy,0-ooo  sin(<ft)cos(y/)hy  -  oo0  sini^cosiy/)!^ 
+  oo  0  cos(</>)s\r\(0)sm(y/)hY  +  ooa  cos((p)sm(0)sm(y/)Ivy,0 

hy  +  Iyy9-6sm{(f))hx  -  0sm((f))lxx<f)->ry/cos(0)cos((l))hx  +  y/  cos(0)cos((f))I ^(j) 
+  6;0sin(^)cos(v/)/?x  +<y<)sin(^)cos(^)/xr^-^0cos(^)sin(6')sin(^/)/2;t  - 
co0  cos(^)  sin(#)  s\n(y/)l ^if)  -  (j>h.  -  iftl^y/  +  y/ sm(0)h.  + y/2  sin(0)I=  + 
oo  0  cos(#)  sin(  y/)h,  +  co0  cos(#)  sin(  (//)!„  y/ 


hz  +  IJj/  +  (f>hy  +  <t>IyyQ-  y/sin(0)hy  -  y/sm(0)ly},0-ooocos(0)sm(y/)hy  - 
ooocos(0)sm(y/)Iyy0-0cos(<p)hx  -  0cos(0)lxx(j>-y/c.os(0)s\n((())hx  - 
ys  cos(0)sm(<f>)I  ^if*  +  oo0cos(<f>)cos(y/)hx  +  coocos(0)cos(y/)Ixx0  + 
co0  sin(^)  sin(#)  sin(y/)hx  +  oo0  sin(^)  sin(#)  siniy/)!^ 


(7-6) 


The  external  moments  arise  from  three  major  sources:  the  gravitational  gradient, 
the  solar  radiation  pressure,  and  the  control  moments  from  actuators.  Denoting  the 
individual  torques  by  MG,  Ms,  and     Mc,  respectively,  we  can  write 


M=MG  +  MS  +  Mc. 


(7-7) 
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where, 


MG  =  3o)20< 


>(4-A,y 

o 


(7-8) 


Different  types  of  active  control  systems  have  different  key  parameters  and 
algorithms.  Frequently,  three-axis  control  can  be  decoupled  into  three  independent  axes. 
The  most  basic  design  parameter  in  each  axis  is  its  position  gain,  Kp .   This  is  the  amount 

of  control  torque  which  results  from  a  unit  attitude  error  and  can  be  expressed  in  Nm/deg 

or  Nm/rad.    The  position  gain  is  selected  by  the  designer  and  must  be  high  enough  to 

provide  the  required  attitude  control  accuracy  in  the  presence   of  disturbances,   or 

T  / 
Kp  >   Yq  ,  where  Kp  is  position  gain,  Td  is  peak  disturbance  torque,  and  0£  is  allowable 

/  £ 

attitude  error.  In  its  simplest  form,  a  spacecraft  attitude  control  system  can  be  represented 
in  the  s-domain  as  a  J/,  plant  and  may  be  controlled  by  a  proportional  plus  derivative 

(PD)  controller  where  Tc  =  Kp0£  + Kr0£.  The  position  gain,  Kp, controls  system 
bandwidth  and  the  rate  gain,  Kr ,  controls  damping.  [Ref.  6:  p.  364] 

The  value  of  the  position  gain  also  determines  the  attitude  control  system 

bandwidth  and  speed-of-response.  The  bandwidth  is  given  by  co„  =  \Kp  I A  2 ,  where  /  is 

the  spacecraft  moment  of  inertia.  The  bandwidth  defines  the  frequency  at  which  control 
authority  begins  to  diminish.  Attitude  control  and  disturbance  rejection  are  effective  from 
0  frequency  (d.c.)  up  to  the  bandwidth.  Speed  of  response  is  approximately  the  reciprocal 
of  bandwidth.  Note  that  position  gain  is  inversely  proportional  to  allowable  error  and 
bandwidth  is  proportional  to  the  square  root  of  position  gain.  Therefore,  high  accuracy 
implies  high  position  gain  and  high  bandwidth.  However,  high  bandwidth  may  cause 
bending  resonances  to  affect  control  system  performance.  [Ref.  6:  p.  365] 
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With  the  relations  given,  the  system  designer  can  estimate  required  position  gain 
from  his  estimates  of  disturbance  torque  and  accuracy  requirements.  He  can  use  this 
estimate  to  compute  control  system  bandwidth.  This  allows  him  to  specify  minimum 
bending  frequencies  as  discussed  below.  [Ref.  6:  p.  365] 

In  defining  algorithms  for  the  control  system,  we  must  also  consider  whether  the 
vehicle  will  have  flexible-body  effects  that  can  make  the  vehicle  unstable.  Spacecraft  with 
flexible  appendages  such  as  antennas,  booms,  and  solar  panels  may  produce  slight  warping 
at  their  natural  frequencies.  Control  torques  and  external-disturbance  torques  will  cause 
structural  vibrations,  in  some  cases  close  to  or  within  the  control  system's  bandwidth.  The 
lowest  natural  frequencies  of  flexible  components  should  be  at  least  an  order  of  magnitude 
greater  than  the  rigid  body  frequencies  before  we  can  neglect  flexibility. 

C.  THREE-AXIS  REACTION  WHEEL  SYSTEM 

A  three-axis  reaction  wheel  system  can  be  considered  as  a  combination  of  three 
independent  roll,  pitch,  and  yaw  control  systems,  see  Figure  7-1.  Each  axis  is  controlled 
by  varying  the  speed  of  the  reaction  wheel  in  response  to  the  attitude  error.  The  system 
requires  a  reaction  wheel  with  zero  nominal  angular  momentum  and  an  attitude  sensor  for 
each  axis.  [Ref.  25:  p.  149] 
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Figure  7-1 :  Block  diagram  of  a  three-axis  reaction  control  system 

1.  Small  Angle  Approximation 

In  this  section  the  control  of  the  spacecraft  is  explored  when  small  angle 
approximations  are  made  which  would  be  the  case  if  yaw  steering  was  not  being 
performed. 

Introducing  Eqs.  (7-4)  and  (7-8)  into  Eq.  (7-7)  and  neglecting  the  transverse 
inertia  of  the  wheel,  the  equations  of  motion  for  a  three-axis  reaction  wheel  system  are 


Kx  +  M„  =  I J  +  ]*a>l(l„  -  4)  -  a>0hyy  +  \-hy  -  *,(/„  -  /„  +  /„)]* 
+  'ehz-(o0hz  +  hx 


A 


M^  +  Msy  =  lj9  +  3o)ld(lxx  -  I~)  +  a>oK<t>  ~  K4>  +  a>oKv  +  K¥x  +  K  (7_9) 

Ma  +  MB  =  ij, + [©;(/„  -  /„)  -  &jtyy + [hy + ©„(/«  +  4  -  4)]^ 

-  h'0  +  (o,hx  +  K 


X  O      X 
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Although  Eqs.   (7-9)  are  coupled,   because    hx,hv,hz,and  coa    are   small,   the 

coupling  terms  are  small.  If  the  coupling  terms  are  neglected,  the  equations  of  motion 
about  the  pitch,  roll,  and  yaw  axes  become  independent,  and  hence  they  can  be  controlled 
independently.  Therefore,  the  equations  of  motion  become: 

Mcx  +  Ma  =  I»d>  +  K  =  Tx 

MQ/  +  M„  =  lJ  +  hy  =  Ty  (7-10) 

Then,  the  control  torques  are  applied  by  letting  the  rates  of  change  of  the  angular 
moment  of  the  reaction  wheels  have 


K  =  Kg(r0O  +  e)  (7-11) 


Making  the  appropriate  substitution  leads  to  the  following: 


Ty  =  lJ+Ky[ry0+0]  (7-12) 


T:  =  IJj/  +  K\r:ij/  +  y/] 
Next,  entering  the  Laplace  domain  yields: 
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«*)- 

[J"S2 

+  Kr 

vrs  +  Kr) 

oQ- 

TA 

») 

'  [  V2 

+  *, 

tps  +  Kp] 

„,M 

U. 

') 

(7-13) 


[i^+KyTyS+K^ 


The  next  step  is  to  determine  the  position  gains  and  the  time  constants.    For  the 
roll  axis  we  can  say: 


co. 


£r  = 


CO, 


(7-14) 


Then,  using  the  final  value  theorem,  <f>ss  =  Unisys)] ,  we  can  state 


.        (2xlQ-4) 

0«  = 

A, 


and  since,    <f>ss  <  0.001745     rad 
K  =  0.5. 


K  >  0.1 146.     Therefore,  we  can  let 


\K. 


Next,  since  mr  =    — -  =  0.0061227     rad  I  sec  =  9.7446e  -  4     hz 


and,  since  rr  = 


2C 


co. 


r  =2060 


In  a  similar  manner  the  remaining  values  were  determined  and  are  listed  in  Table  7- 1 . 
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ROLL 

PITCH 

YAW 

Kr  =  0.5 

Kp  =  0.3 

Ky  =  0.1 

o)r  =  0.0061  rad/sec 

o)p  =  0.0046  rad/sec 

<yv  =  0.00266  rad/sec 

£r=    1 

<T,=  1 

C>  =  i 

rr  -  2060  sec 

r^  =  433  sec 

rv.  =  751 

Table  7- 1 :  Control  Law  Values 

The  final  step  is  to  code  in  Matlab  and  simulink.  The  Laplace  equations  I  used  for 
this  simulation  are  as  follows: 


Tx(s) -  K\rrs  +  \]^s)  +  a>0[l„  -  I„  +  I=]sy,(s)  =  [i^s2  +  4a>20(l)y  -  I=)]<?>(s) 


Ty(s)  ~  KpVps  +  {}°(S)  =  [ V2  +  Wo(J~  -  L)fa) 


(7-15) 


T2(s)  -  Ky[zys  +  \]¥(s)  -  a}0(l„  +  I=-  I„)s4(s)  =  [l=s2  +  m2^  -  I^s) 


For  the  computer  simulation  I  used  Matlab  and  Simulink.  1  considered  three 
situations.  First,  a  step  input  due  to  thruster  misalignment.  Second,  the  sinusoidal  input 
due  to  solar  pressure  torques;  and  third,  an  impulse  input  of  approximately  one  order  of 
magnitude  larger  than  the  others;  see  simulink  diagrams  and  Matlab  code  in  Appendix  F. 

The  first  simulation  was  the  pitch  response  due  to  a  step  input.  As  can  be  seen 
from  Figure  7-2,  a  step  input  occurred  at  time,  t=100.  Then,  the  control  system  was  able 
to  counterbalance  the  input  with  in  1300  seconds  and  remain  with  in  a  pointing  tolerance 
of  0.1°.  But,  this  requires  a  continual  supply  of  angular  momentum  from  the  reaction 
wheel.  A  momentum  dump  will  need  to  occur  after  850,000  seconds  or  every  9.8  days  for 
the  momentum  wheels  selected.  Next,  with  a  sinusoidal  torque  input  the  system  remained 
within  pointing  accuracy  and  fluctuated  between  ±0.05°.  Finally,  with  an  impulse  input  of 
approximately  one  order  of  magnitude  larger  than  the  other  torques  the  system  stabled  out 
after  1600  seconds  and  also  remained  within  the  required  pointing  accuracy,  see  Figure  7- 
4. 
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The  second  simulation  was  of  the  roll  and  yaw  components,  which  are  coupled 
(see  Appendix  F).  Other  than  the  coupling  difference  the  simulation  process  was  the  same 
as  for  the  pitch.  Graphs  of  the  simulation  can  be  seen  in  Figures  7-5  through  7-8.  For  the 
sinusoidal  input  for  roll  the  angular  momentum  did  not  increase  with  time  but  instead  was 
sinusoidal.  This  is  significant  because  a  momentum  dump  will  not  be  required.  For  the 
sinusoidal  input  for  yaw  a  momentum  dump  must  occur  every  800,000  seconds  or  every 
9.25  days. 

The  final  simulation  was  that  for  the  contingency  mode  of  operation  which 
simulated  a  yaw  control  failure  (see  Appendix  F).  For  this  I  just  used  a  fourth  reaction 
wheel  canted  between  the  roll  and  yaw  directions.  Then,  I  ran  all  three  inputs  to  show 
that  the  system  was  able  to  control  the  spacecraft  and  remain  within  pointing  accuracies. 
This  information  is  displayed  in  Figures  7-9  through  7-12. 

Finally,  in  estimating  the  fuel  required  for  momentum  dumping  I  tended  to  be 
conservative.  For  the  estimation  I  am  assuming  that  both  pitch  and  yaw  angular 
momentum  must  be  dumped  every  10.5  days  and  9.25  days  respectively.  Therefore,  we 
will  need  approximately  80  pulses  in  a  year. 

Next,  from  the  following  relationships  we  can  determine  the  amount  of  propellant 
mass  used  in  each  burn: 


hPg 


Therefore,  the  total  propellant  mass  used  in  one  year  for  thruster  only  momentum 
dumping  is  estimated  to  be  1.5  kg. 
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Figure  7-2:  Pitch  Response  Due  to  a  Step  Input  Disturbance  Torque. 
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Figure  7-3:  Pitch  Response  Due  to  a  Sinusoidal  Disturbance  Torque. 
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Figure  7-4:  Pitch  Response  Due  to  an  Impulse  Disturbance  Torque. 
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Figure  7-5:  Roll  Response  Due  to  a  Sinusoidal  Disturbance  Torque. 
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Figure  7-6:  Yaw  Response  Due  to  a  Sinusoidal  Disturbance  Torque. 


x  10 


Roll  Torque  vs  Time 


o- 


0  1000        2000        3000        4000         5000        6000         7000        8000        9000        10000 

Time  (sec) 


x  10 


Roll  angle  phi 


degree 

s 

r 

c 
o 
a 
in   *. 

• 

|1 

CD 

c 
«I 
=;  0 

OC 

E 

c 
• 
E 
o 

fc-2 

a 

?-3 

«*  0  1000        2000        3000        4000        5000        6000         7000        8000        9000        10000 

Time  (sec) 


x  10 


1000        2000        3000        4000        5000        6000         7000        8000        9000       10000 

Time  (sec) 
Angular  Momentum  vs  Time 


-1 


N-m-s 

Figure  7-7:  Roll  Response  Due  to  an  Impulse  Disturbance  Torque. 
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Figure  7-8:  Yaw  Response  Due  to  an  Impulse  Disturbance  Torque. 
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Figure  7-9:  Roll  Response  Due  to  a  Step  Input  Disturbance  Torque  -  Back  Up  System. 


x  70"* 
g 

Yaw  Torque  vs  Tim 

e 

4 

9 

n 

c 
o 

2  0.06 


0.04 


0.02 


> 
E 

c 

o 

5 


*-0.5 


1000        2000         3000 


4000         5000         6000 

Time  (sec) 

Yaw  angle  psi 


7000         8000        9000        10000 


" 

degret 

>s 

1000        2000        3000 


4000         5000         6000 

Time  (sec) 

Angular  Momentum  vs  Time 


7000         8000        9000        10000 


N-m-s 

0  1000        2000        3000        4000         5000        6000         7000         8000        9000        10000 

Time  (sec) 


Figure  7-10:  Yaw  Response  Due  to  a  Step  Input  Disturbance  Torque  -  Back  Up  System. 
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Figure  7-11:  Roll  Response  Due  to  an  Impulse  Disturbance  Torque  -  Back  Up  System. 
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Figure  7-12:  Yaw  Response  Due  to  an  Impulse  Disturbance  Torque  -  Back  Up  System. 
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2.  Yaw  Steering  Requirements 

In  this  section  the  attitude  dynamics  and  control  of  the  spacecraft  is  explored  when 
small  angle  approximations  are  not  made  and  all  higher  order  terms  are  included  in  the 
equations  of  motion.  This  is  a  complexity  that  is  a  result  of  the  requirements  for  yaw 
steering. 

Introducing  Eqs.  (7-6)  and  (7-8)  into  Eq.  (7-7),  the  equations  of  motion  for  a  three 
axis  reaction  wheel  system  are 

Mcx  +  Ma  =  I^ij)  +  6cos((f>)I=ii/  +  if/1  cos(#)sin(^)/_  -  cou  cos(^)cos(y/)/_y/  - 
co0  sin(#)sin(0)sin(^)/_^  +  02  sin(^)/VT  -  ^/cos(^)cos(^)/KV,6'  - 
^osin(0)cos(^)/w0  +  <yocos(0)sin(6Osin(^)/j9- 
3co20  sin(^)cos(^)[cos(^)]"/_.  +  3co20  sin(^)cos(^)[cos(#)]  I,y  +  hx 
+  6cos{(p)hz  +  ii/cos(0)sm((f>)h:  -  coocos(0)cos(y/)h:  - 
<y0sin(^)sin(#)sin(y/)/7_  +  0sm(<f>)hv  -  i}/ cos(0)  cos(<fi)hv  - 
0)o  sin(^)  cos(iy)hy  +  co0  cos(^)  sin(#)  sm(y/)hv 

M^  +  M^  =  Iyy0-0sin(0)Ixx0+  i}/cos(0)cos((/))Ixj  +  coosm(^)cos(i//)Ixj- 

-  co0  cos(^)  sin(^)  sin( ^r)/JCX^  -  <f£B\jf  +  y/2  sin(#)/„  + 
cd0  cos(#)  sin(^)/=^  -  3co20  sin(#)cos(#)cos(^)/,,  + 

3a>20  sin(^)cos(^)cos(^)/xl  +  hy  -  0sm($)hx  +  y/ cos(0)  cos(^>)hx  + 

G)osm(0)cos(y/)hx  -ct)ocos(<fi)sm(0)sm(y/)hx  - ifkz  +  y/sin(0)h:  + 
a)  0  cos(0)  sin(y/)hz 

(7-17) 

Ma  +  MS2  =  I  J,  +  <f>lj  -  y,  sm(0)lj  -  co0  cos(0)  sm{¥)lj  -  0cos(<f>)lj  - 

y/  cos(#)  sin(  0)7^0  +  co0  cos(^)  cos(^)  1 ' J$  +  cov  sin(^)  sin(6>)  s\n{if/)  I  ^ip 

-  3co20  sin(6>)cos(^)sin(^)/xx  +  3co20  sin(0)cos(0)sin(0)/„,  +  hz  +  if)hy  - 
ij/sm(0)hv  -  co0  cos(0)sin(^)/7v  -  0cos{<f>)hx  -  ij/ cos(0)  sin(<p)hx  + 
co0  cos((f>)cos(y/)hx  +  0)o  sSn{(f>)sm{0)sm(y/)hx 

These  equations  of  motion  are  somewhat  overwhelming.  In  the  process  of 
developing  a  control  law  the  first  step  was  to  transform  these  equations  of  motion  into  a 
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state  space  formulation  and  develop  a  Matlab/Simulink  model  for  the  motion  of  the 
spacecraft  with  no  control.  The  state  space  formulation  used  was  the  following: 


*1  =  <l> 

JC,  =  JC,  =  (j) 

x,  =  e 

x4  =  Xj  =  0 

x5  =  y/ 

x6  =  x5  =  i// 

xx=x2 

~-l~ 

x2  =  xl  =  <}>  = 

J». 

(Ocos(0)I„y/  +  ^cos^sinC^)/^  -  co()cos((/>)cos(y/)Iz2ij/  - 

coosin(</>)sm(0)sm(y/)l:::i//  +  &  sin^)/^,  -  y/cos(0)cos(</>)Ivy,0- 
co0  sin(^)  cos(  \f/}lj&  +  0)o  cos(^)  sin( 0)  s\n{y/)Ixy0  - 
3co20  sin(^)cos(0)[cos(<9)f  Is  +  3co20  sin(^)cos(^)[cos(<9)]2 1„ 
+  hx  +  0cos(<f>)h,  +  y/cos(0)sin(^>)hz  -  6>0cos(^)cos(y/)/i,  - 
co0  sin(^)sin(0)sin(^)/z,  +  9sm(<f>)hy  -  y/cos(0)cos(^)hy  - 
co osin(0)cos(y/)hv  + co ocos(0)sin(O)sm(y/)hv  -  Mxd) 


h  =  X4 

K4  =  Xj  =  0  — 

-1 

Jyy. 

(7-18) 


(-6sm((f>)I ^(j)  +  tj/cos(0)  cos^)/^  +  a>0  sin(^)  cos^)/^  - 

co0  cos(^)sin(^)sin(^-)/xr^  -  <t>IJi/  +  co0  cos(0)sm(y/)Iz.i//  -  <j)hz 

-  3o)20  sin(#)cos(#)cos(^)/,_  +  3o)20  sm(0)cos(9)cos(<p)Ixx  +  hv 

-  6s\n((j))hx  +  {f/cos{6)cos>{<f>)hx  +  co0  sm{^)cos(if/)hx  +  \j/2  sin(#)/,_  - 
cooo,os{<f))sm(0)s\n{\i/)hx  +  y/sm(0)h,  +co„cos(0)sm(y/)h,  -  M2d) 
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x5  —  x6 


X6=X5   =   W  = 


-1 


(flyfi  -  ¥  s^WyyO  -  to 0  cos(0)  sin(  yr)/^  -  6cos(<f>)lJ 


y/  cos(6>)  sin^)/^  +  co0  cos(»  cos(y/)I ^  +  co0  sin(^)  sin(0)  sin(^)4^ 
-  W0  sin(6>)cos(6>)sin(^)/xx  +  3^  sin(6>)cos(6>)sin(^)/v:v  +  hz 
+  <f>hy-¥  sm(0)hv  -  co0  cos(6>)  sin(»/?y  -  9cos(</>)hx  -  ^  cos(0)  sin(^)^ 
+  0)o  cos(<f>)cos(i//)hx  +  co0  sin(0)sin(0)sin(^)/?x  -  M,d) 


Then,  after  developing  the  simulink  model,  I  analyzed  the  dynamics  of  the 
spacecraft  as  a  function  of  time  when  a  disturbance  torque  was  applied.  This  can  be  seen 
in  Figures  7-13  through  7-15,  also,  see  Appendix  G.  Note  that  without  control  the 
spacecraft  will  become  unstable  very  quickly. 
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Figure  7-13:  Disturbance  Torques  Applied  to  the  Spacecraft  -  Large  Angle  Scenario. 
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Figure  7-14:  The  three  body  components  of  angular  velocity  due  to  a  disturbance  torque 

with  no  control  -  Large  Angle  Scenario. 
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Figure  7-15:   Spacecraft  Angular  Response  Due  to  Disturbance  Torque  -  Large 

Angle  Scenario. 
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D.  ADCS  CONTROL  AND  PERFORMANCE  ANALYSIS  -  YAW 

STEERING 

In  order  to  keep  the  antenna  beam  pointing  accuracy  and  the  solar  array  panel 
pointing  accuracy  between  acceptable  limits  (0.1°  and  1°  respectively),  a  control  system 
must  determine  the  required  movements  of  the  reaction  wheel.  This  controller  (control 
loop)  could  be  developed  in  many  different  ways.  In  our  case  the  angle  and  the  angular 
speed  measured  by  the  sensors  are  used.  In  this  way  the  system  corrects  the  deviation 
based  on  the  actual  deviation  and  the  rotation  speed.  The  control  loop  is  given  in  Figure 
7-16. 
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Figure  7-16:  The  Control  Loop 

The  block  scheme  of  the  system  model  is  shown  in  Figure  7-17.  The  satellite 
dynamics  are  described  by  the  equations  7-16.  The  BAPTA  controller  controls  the  actions 
of  the  BAPTA,  which  again  influences  the  motions  about  the  y-axis.  There  are  three  RW 
controllers  each  controlling  a  reaction  wheel.  The  OBC  has  two  sensor  inputs,  the  orbital 
angle  P  and  the  angle  between  the  orbital  plane  and  the  Sun-Earth-vector  8. 
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Figure  7-17:  System  Model 

For  the  computer  simulation  Matlab  and  Simulink  were  used  (see  Appendix  G). 
And,  three  situations  were  considered.  First,  a  step  input  due  to  thruster  misalignment. 
Second,  the  sinusoidal  input  due  to  solar  pressure  torques;  and  third,  an  impulse  input  of 
approximately  one  order  of  magnitude  larger  than  the  others  (see  simulink  diagrams  and 
Matlab  code  in  Appendix  G).  Several  combinations  of  the  three  scenarios  were  also 
tested.  When  modeling  the  simulations,  disturbance  torque  values  and  scenarios  were 
used  that  rarely  would  occur  on  the  satellite.  This  was  done  in  an  effort  to  show  the 
robustness  of  the  satellite  control.  In  general,  we  are  trying  to  look  at  the  worst  case 
scenarios.  Also,  the  control  values  selected  are  indicated  in  Table  7-2. 
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ROLL 

PITCH 

YAW 

£r  =8.5 

Kp  =9.0 

Ky=\ 

rr=2060 

rp=2000 

rv=751 

Table  7-2:  Control  Law  Values 

The  first  simulation  imparted  impulse  torques  along  the  roll,  pitch  and  yaw  axes. 
The  magnitudes  selected  were  greater  than  those  expected  to  actually  occur  on  orbit  and  it 
be  would  extremely  unlikely  to  have  3  simultaneous  impulse  torques.  As  can  be  seen  in 
Figure  7-18,  the  impulse  torques  occurred  at  time,  t=100  seconds.  Figure  7-20  illustrates 
the  angular  responses  for  roll,  pitch,  and  yaw.  A  thorough  study  of  the  responses  reveal 
the  greatest  error  occurs  from  the  effects  of  yaw  steering  and  not  the  impulse  torques. 
And,  the  required  pointing  accuracies  have  been  maintained  throughout  the  orbit.  Next, 
with  a  step  input  along  all  three  axes,  the  angular  responses  were  evaluated  and  found  to 
be  within  operational  limitations.  It  should  be  noted  that  the  greatest  error  occurs  during 
this  scenario.  This  is  because  the  step  disturbance  is  occurring  at  the  same  time  the  yaw 
inversion  is  taking  place.  So  in  effect,  we  have  an  additive  disturbance  among  the  two. 
Then,  a  sinusoidal  disturbance  torque  was  applied  along  all  three  axis  simultaneously. 
This  can  be  seen  in  Figures  7-21  through  7-23.  And,  once  again  the  pointing  accuracies 
fell  within  limits.  Finally,  a  simulation  of  sinusoidal,  impulse,  and  step  disturbance  torques 
were  applied  along  the  roll,  pitch,  and  yaw  axes  respectively.  And,  in  this  scenario  all 
pointing  accuracy  requirements  were  maintained. 

Throughout  the  simulations  the  maximum  error  occurred  when  the  yaw  inversion 
took  place.  It  is  at  this  point  that  the  yaw  angular  velocity  is  greatest.  This  is  due  to  the 
fact  that  the  yaw  angle  and  angular  rate  are  coupled  into  the  roll  and  pitch.  So,  when  yaw 
rate  of  change  is  greatest  we  see  our  greatest  error  along  all  three  attitude  axes.  Also,  the 
yaw  angular  velocity  and  angular  acceleration  is  not  continuous  at  their  maximum  and 
minimum  values.  This  in  turn  creates  a  rapid  change  condition  which  then  allows  for  more 
error  to  be  introduced. 

The  final  simulation  was  a  sinusoidal  input  along  the  pitch  axis  and  was  performed 
to  illustrate  the  value  of  disturbance  torque  that  would  cause  the  system  to  exceed 
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pointing  accuracy  requirements.  As  can  be  seen  in  Figures  7-24  through  7-26,  a  sinusoidal 
input  of  0.01  N-m  caused  the  pitch  axis  pointing  accuracy  to  exceed  0.1°.  Recall  that  the 
maximum  disturbance  torque  anticipated  is  0.001  N-m.  Therefore  a  built  in  safety  margin 
exists. 
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Figure  7-18:  Impulse  Disturbance  Torque  for  Control  of  Yaw  Steering  Scenario  #1 
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Figure  7-19:  Spacecraft  Angular  Velocities  During  the  Control  of  Yaw  Steering  Scenario 

#1  (Time  in  seconds). 
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Figure  7-20:   Spacecraft  Angular  Response  Due  to  Impulse  Disturbance  Torque  -  Yaw 

Steering  Scenario  #1  (Time  in  seconds). 


131 


x  10 


-1 


-1 


0.5 


Disturbance  Torques  -  Tx,  Ty,  Tz 


1  1.5 

Time  (sec) 


0                              0.5                              1                               1.5                              2                              2.5 

Time  (sec)                                                                  .  „■» 
^                                                                                  '                                                                 x  10 
x  10 

^^~ 

0                              0.5                              1                               1.5                              2                              2.5 

Time  (sec)                                                                      ..4 
*                                                                                  '                                                                 x  10 
x  10 

"""        ~~~ 

^ 

2.5 

x  10" 


Figure  7-21:  Sinusoidal  Disturbance  Torque  for  Yaw  Steering  Scenario  #2. 
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Figure  7-22:  Spacecraft  Angular  Velocities  During  the  Control  of  Yaw  Steering  Scenario 

#2. 
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Figure  7-23:   Spacecraft  Angular  Response  Due  to  Sinusoidal  Disturbance  Torque  -  Yaw 

Steering  Scenario  #2. 


133 


Disturbance  Torques  -  Tx,  Ty,  Tz 


5 


E 

£    C 


-1 


0.01 


-0.01 


-1 


0.5 


1  1.5 

Time  (sec) 


0                              0.5                               1                               1.5                              2                              2.5 

Time  (sec)                                                              x  1Q4 

0                              0.5                               1                               1.5                              2                              2.5 

Time  (sec)                                                              x  1Q4 

x  10 


2.5 

4 


Figure  7-24:   Sinusoidal  Disturbance  Torque  for  Yaw  Steering  Scenario  #3. 
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Figure  7-25:  Spacecraft  Angular  Velocities  During  the  Control  of  Yaw  Steering  Scenario 

#3. 
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Figure  7-26:  Spacecraft  Angular  Response  Due  to  Sinusoidal  Disturbance  Torque  -  Yaw 

Steering  Scenario  #3 . 
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VIII.  CONCLUSION 

A.  SUMMARY  OF  DESIGN  PROCESS 

In  this  thesis  the  preliminary  design  of  an  Attitude  Dynamics  and  Control 
Subsystem  for  a  Medium  Earth  Orbit  communications  satellite  was  performed.  The 
process  focused  on  the  six  critical  steps  in  designing  an  ADCS  which  are  stated  as  follows: 

1 .  Design  control  modes. 

2.  Define  or  derive  system  level  requirements. 

3.  Select  type  of  spacecraft  control  by  attitude  control  mode. 

4.  Quantify  the  disturbance  environment. 

5.  Select  and  size  ADCS  hardware. 

6.  Define  determination  and  control  algorithms. 

Also  presented  was  the  important  aspect  of  how  the  ADCS  is  inter-related  with  the 
other  spacecraft  subsystems.  And,  how  many  of  the  requirements  and  design  drivers  of 
the  ADCS  originate  in  the  other  spacecraft  subsystems. 

The  highlights  and  conclusions  of  the  ADCS  design  process  are  summarized  as 
follows.  First,  the  control  modes  of  the  spacecraft  were  determined  and  defined.  The 
control  modes  determined  to  be  most  appropriate  for  our  satellite  were  the  following: 

1 .  Orbit  insertion  mode 

2.  Acquisition  mode 

3.  Sun  acquisition/sun  hold  mode 

4.  Earth  acquisition  mode 

5.  Normal  mode 

6.  Slew  mode 

7.  Contingency  mode 

8.  Orbit  transition  mode 
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The  second  design  milestone  was  the  determination  of  system  level  requirements. 
Ten  ADCS  requirements  were  presented  which  are  summarized  below  as  follows: 

1.  The  ADCS  must  be  fully  operational  without  the  aid  of  ground  base  control. 

2.  The  operational  lifetime  of  the  ADCS  must  meet  or  exceed  ten  years. 

3.  The  ADCS  must  maintain  a  UHF  antenna  beam  pointing  accuracy  of  less  than 
or  equal  to  0.1°. 

4.  The  ADCS  must  have  the  capability  to  provide  continuous  yaw  steering  and 
solar  array  tracking  such  that  the  solar  array  pointing  accuracy  is  better  than  or  equal  to 
5°. 

5.  The  spacecraft  must  be  fully  operational  during  eclipse. 

6.  The  spacecraft  must  be  stabilized  and  under  control  during  all  phases  of  the 
mission. 

7.  The  ADCS  system  should  have  system  redundancy  such  that  there  are  no  single 
point  failures. 

8.  The  ADCS  should  have  a  ground  override  capability. 

9.  The  spacecraft  must  have  a  cross-link  antenna  pointing  accuracy  of  0.07°  which 
requires  a  spacecraft  attitude  knowledge  requirement  of  0.02°  in  roll,  pitch,  and  yaw. 

10.  The  ADCS  must  have  a  positional  knowledge  requirement  of  ±  350  meters. 

The  next  design  milestone  was  deteirnining  the  type  of  spacecraft  control  by 
attitude  control  mode  and  presenting  a  functional  description  of  the  ADCS.  During  this 
stage  we  determined  what  kind  of  sensors  and  actuators  were  to  be  used  and  then  showed 
how  they  were  to  be  used  in  each  mode  of  operation.  Also  introduced  in  this  stage  of  the 
design  was  the  motion  of  the  satellite  and  solar  arrays  due  to  the  yaw  steering 
requirement. 

Then,  the  design  focused  on  quantifying  the  disturbance  environment.  In  this 
process  the  derivation  of  the  Moments  due  to  solar  pressure,  gravity  gradient, 
aerodynamic  and  geomagnetic  were  presented.  And,  an  estimated  disturbance  value  was 
attained.    Next,  the  selection  of  the  ADCS  hardware  was  performed.    Each  sensor  and 


138 


actuator  was  described  and  the  reasons  for  selection  of  each  component  were  given.   The 
sensor  and  actuator  selection  was  driven  solely  by  the  requirements  initially  stated. 

The  final  step  in  the  design  process  was  to  define  determination  and  control 
algorithms.  Chapter  VI  examined  yaw  controlling  and  presented  explanations  dealing  with 
control  strategy  and  control  software.  Then,  in  chapter  VII  several  simulations  of  the 
ADCS  were  presented  and  the  development  of  the  control  law  was  explained. 

B.  SUMMARY  OF  SPACECRAFT  ADCS  DESIGN 

The  ADCS  chosen  for  the  MUS  satellite  was  a  Zero  Momentum  three  axis  system 
which  can  provide  any  combination  of  slew  maneuvers  in  the  roll,  pitch,  and  yaw 
directions.  The  ADCS  can  maintain  Nadir  pointing  accuracies  to  within  0.1°.  At  the  same 
time  it  can  continuously  rotate  in  the  yaw  direction  while  the  solar  arrays  are 
simultaneously  rotating  in  order  to  keep  a  perpendicular  orientation  towards  the  sun.  This 
will  greatly  decrease  the  complexity  of  the  solar  array  control  mechanisms  and  reduce 
spacecraft  mass. 

As  previously  mentioned  the  spacecraft  ADCS  has  eight  control  modes.  Each 
control  mode  has  peculiar  operating  characteristics  which  enable  the  spacecraft  to 
maintain  operations  for  ten  or  more  years.  The  operational  requirements  of  the  control 
modes  were  satisfied  by  the  utilization  of  various  ADCS  components  as  presented  in 
Figure  8-1. 
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Figure  8-1 :  Block  Diagram  of  the  Attitude  Dynamics  and  Control  Subsystem 

The  ADCS  components  selected  were  those  that  satisfied  the  subsystem 
requirements  as  stated  in  the  previous  section.  The  components  selected  for  the  MUS 
spacecraft  are  listed  in  Table  8-1 . 
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Unit 

Quantity 

Mass  (ea) 

Avg  Power 

ADCS  Control 
Processor 

2 

5.7  kg 

25  W 

Control  Electronics 

1 

0.9  kg 

2.0  W 

Magnetometer 

1 

0.5  kg 

0.7  W 

Torque  Rods 

3 

0.9  kg 

— 

Reaction  Wheel 

1 

14.3  kg 

28  W 

Reaction  Wheel 

3 

10.5  kg 

30  W 

Star  Trackers 

3 

3.63  kg 

9  W 

Inertial  Reference 
Unit 

2 

3.45  kg 

34  W 

Sun  Sensor 

2 

0.64  kg 

1  W 

GPS 

2 

1.68  kg 

7W 

GPS  Antenna 

4 

0.2  kg 

— 

Earth/Sun  Sensor 

2 

2.5  kg 

5.5  W 

Table  8-1 :  SpacecraJ 

t  ADCS  Components 

After  the  ADCS  component  selection  had  taken  place,  the  satellite  bus  design  was 
able  to  proceed  by  incorporating  the  equipment.  The  final  ADCS  configuration  for  the 
MUS  design  is  shown  in  Figure  8-2. 
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Figure  8-2:  ADCS  Component  Layout 

Finally,  a  simulation  of  the  spacecraft  control  was  performed.  When  expected 
disturbance  torques  and  yaw  steering  requirements  were  imposed  on  the  model;  pointing 
accuracy  requirements  were  maintained.  The  limits  of  the  subsystem  were  tested  and  it 
was  concluded  that  it  would  require  disturbance  torques  approximately  one  order  of 
magnitude  larger  than  what  was  expected  to  cause  spacecraft  operating  limits  to  be 
exceeded. 
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APPENDIX  A 


ACRONYMS 


ADCS 

AFSCN 

BAPTA 

CG 

CMG 

CPU 

COTS 

dB 

DoD 

DT 

DSCS 

DuSD  -  Space 

EHF 

EPS 

ES 

ES/SSA 

GEO 

GG 

GHz 

GPS 

HF 

IRU 

Kbps 

LEO 

MEO 

MHz 

MILSATCOM 

MOI 

MSS 

MUS 

Nm 

Nms 

PMOI 

Pri 

RF 

RW 

SA 

SCOC 

Sec 

SHF 

SOCC 


Attitude  Dynamics  and  Control  System 

Air  Force  Satellite  Control  Network 

Bearing  and  Power  Transfer  Assembly 

Center  of  Gravity 

Control  Moment  Gyro 

Central  Processing  Unit 

Commercial  Off  The  Shelf 

Decibels 

Department  of  Defense 

Disturbance  Torque 

Defense  Satellite  Communications  System 

Deputy  Undersecretary  of  Defense  for  Space 

Extremely  High  Frequency 

Electrical  Power  Subsystem 

Earth  Sensor 

Earth  Sensor  /  Sun  Sensor  Assembly 

Geostationary  Earth  Orbit 

Gravity  Gradient 

Gigahertz 

Global  Positioning  System 

High  Frequency 

Inertial  Reference  Unit 

Kilobits  per  second 

Low  Earth  Orbit 

Medium  Earth  Orbit 

Megahertz 

Military  Satellite  Communication 

Moments  of  Inertia 

Mobile  Satellite  Services 

Mobile  Users  Subscribers 

Newton  Meter 

Newton  Meter  Second 

Principal  Moments  of  Inertia 

Primary 

Radio  Frequency 

Reaction  Wheel 

Solar  Arrays 

Satellite  Control  Operations  Center 

Secondary 

Super  High  Frequency 

Satellite  Operational  Control  Center 
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SS  Sun  Sensor/Star  Sensor 

TT&C  Telemetry,  Tracking  and  Control 

UFO  UFH  Follow  On 

UHF  Ultra  High  Frequency 

VHF  Very  High  Frequency 


144 


APPENDIX  B 
GLOSSARY 

Angular  Momentum.  The  vector  product  of  the  position  vector  and  linear  velocity  of  a 
particle  in  motion  relative  to  an  axis. 

Cyclic  Torque.  Torque  that  repeats  itself  in  cycles  such  that  a  build  up  of  angular 
momentum  does  not  occur. 

Desaturation.  The  process  of  removing  built  up  angular  momentum  from  the  reaction 
wheels. 

Disturbance  Torque.  A  torque  applied  to  the  spacecraft  such  that  its  attitude  is 
disturbed. 

Geostationary  Orbit.  A  geosynchronous  orbit  with  an  inclination  of  0°  or  near  0°, 
positioning  the  satellite  directly  above  the  equator.  [Ref  16:  p.  2] 

Geosynchronous  Orbit.  An  orbit  with  a  period  equal  to  the  earth's  rotational  period  of 
1,436  minutes  and  an  altitude  above  the  earth  of-  19300  nm.  Most  geosynchronous 
satellites  move  in  near  circular,  near  equatorial  orbits.  [Ref  16:  p.  2] 

Global  Coverage.  Coverage  of  all  latitudes  and  longitudes,  90°  N  to  90°  S. 

ICO.  A  medium  earth  orbit  constellation  consisting  of  10  active  satellites  and  2  on  orbit 
spares,  which  will  operate  at  10,400  km,  providing  worldwide  coverage.  The  system  will 
consist  of  2  orbital  planes  with  an  inclination  of  45°.  Proposed  by  ICO  Global 
Communications,  IOC  for  the  system  is  planned  for  the  4th  quarter  of  1999,  with  FOC 
scheduled  for  the  3rd  quarter  of  2000.  ICO,  which  will  use  TDMA  for  modulation,  will 
support  4,500  simultaneous  handheld  users  per  satellite,  and  can  also  support  aeronautical, 
fixed,  mobile  and  pager  users.  Satellite  service  life  is  projected  at  10  years.  The  system 
will  operate  in  S-band  and  C-band,  supporting  data  rates  of  2.4  or  4.8  kbps.  Space 
segment  cost  is  estimated  to  be  2.6  Billion  (including  7  earth  stations  and  3  gateway 
switches),  estimated  handset  cost  is  $1,000,  and  projected  operating  cost  is  $1.00- 
$2.00/minute  in  addition  to  a  $40.00  per  month  access  charge.  [Ref.  16:  p.  3] 

Low  Data  Rate.  For  the  purposes  of  the  Mobile  User  Study,  64  Kbps  and  below. 
[Ref.  16:  p.  4] 

Medium  Earth  Orbit.  Implies  an  orbital  altitude  between  Low  Earth  Orbit  (LEO  ~ 
lOOnm  to  540nm)  and  Geosynchronous  Orbit  (GEO  ~  18225nm).  Orbit  of  choice  for 
many  proposed  Mobile  Satellite  Systems.  [Ref.  16:  p.  5] 
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Mobile  Satellite  Service  (MSS).  The  military  use  of  satellite  based  commercial 
communications  service  providing  digital  voice,  data,  paging  or  fax  services  to  users 
equipped  with  small  (typically  handheld)  terminals.  For  the  purpose  of  the  MUS,  this 
definition  refers  to  Military  use  of  PCS. 
[Ref.  16:  p.  5] 

Mobile  User.  A  highly  mobile  and  or  disadvantaged  user  that  requires  a  small  terminal 
and  or  antenna  (includes  hand  held  terminals  and  terminals  for  platforms  such  as  trucks, 
submarines,  small  ships,  boats,  and  aircraft).  [Ref.  16:  p.  5] 

Nadir.  An  earth  pointing  orientation  or  in  the  direction  of  the  earth. 

Odvssev.  A  medium  earth  orbit  constellation  consisting  of  12  active  satellites  with  spares 
launched  as  required,  which  will  operate  at  10,355  km  (circular  orbit),  providing 
worldwide  coverage  at  all  latitudes.  The  system  will  consist  of  3  orbital  planes  with  an 
inclination  of  50°.  Proposed  by  TRW,  the  system  is  planned  for  IOC  in  the  4th  quarter  of 
2001,  with  FOC  planned  for  the  4th  quarter  of  2002.  ODYSSEY  will  use  CDMA  for 
modulation,  and  will  support  9,500  simultaneous  handheld  users  per  satellite,  as  well  as 
fixed  and  transportable  users.  Space  segment  cost  is  projected  to  be  3.2  Billion,  estimated 
handset  cost  is  $500-5700,  and  projected  operating  cost  is  less  than  $1.00/minute.  [Ref. 
16:  p.  6] 

Personal  Communications  System  (PCS).  Satellite  based  commercial  communications 
service  providing  digital  voice,  data,  paging  or  fax  services  to  users  equipped  with  small 
(typically  handheld)  terminals.  [Ref.  16:  p.  6] 

Polar  Covera2e.  Coverage  that  extends  to  the  extreme  northern  and  southern  polar 
regions  (90°  N  to  90°  S).  [Ref.  16:  p.  6] 

Secular  Torque.  A  torque  that  is  continuous  thus  causing  a  build  up  of  Angular 
Momentum  to  occur. 

Slew.  A  reorientation  of  the  spacecraft  about  the  roll,  pitch  or  yaw  axes. 

A  -  Area. 

P  -  The  angle  in  the  orbital  plane  between  the  point  with  minimum  distance  (sub  solar 
point)  to  the  sun  and  the  satellite. 

P  -  The  Geomagnetic  Field. 
Cd  -  Coefficient  of  drag. 
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dr  -  Vectors  from  the  center  of  mass  to  each  external  piece's  center  of  pressure. 

Ms  -  Moment  due  to  solar  radiation  pressure. 

ps  -  Photons  that  are  specularly  reflected  from  a  surface. 
pd  -  Photons  that  are  diffusely  reflected  from  a  surface. 

pa  -  Photons  absorbed  by  a  surface. 

p  -  Atmospheric  density. 

8  -  The  angle  between  the  orbital  plane  and  the  direction  of  the  sun. 

8,  -  Critical  value  of  8  where  an  alternative  control  scheme  must  take  place. 

Fa  -  Radiation  force  due  to  absorbed  photons. 

Fd  -  Force  due  to  solar  radiation  pressure. 

Fs  -  Radiation  force  due  to  specularly  reflected  photons. 

i_  -  Satellite  orbit  inclination. 

/     -  Vector  distance  from  center  of  mass  to  center  of  pressure. 

77  -  Solar  Array  angle  of  rotation. 
uv  -  Unit  vector  in  velocity  direction. 

s_  -  Solar  radiation  unit  vector. 

Q-  Right  Ascension  of  the  ascending  node. 

T(a)  -  Spacecraft  torque  due  to  areodynamic  effects. 

T(s)  -  Spacecraft  torque  due  to  gravity  effects. 

T(m)  -  Spacecraft  torque  due  to  the  magnetic  field  effects. 
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Tx  -  Torque  about  the  spacecraft's  roll  axis. 

T  -  Torque  about  the  spacecraft's  pitch  axis. 

Tr  -  Torque  about  the  spacecraft's  yaw  axis. 

v  -  Earth  centered  true  anomaly. 

Xs  -  The  satellite's  roll  axis. 

Ys  -  The  satellite's  pitch  axis. 

Zs  -  axis  -  Yaw  axis  or  axis  parallel  to  spacecraft  UHF  antenna 

y/  -  Yaw  rotation  angle. 
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APPENDIX  C 
MATLAB  CODE  FOR  YAW  STEERING 

%  Danny  Busch  21  August  1997 

%AA4871 

%  This  m  file  plots  the  Yaw  Rotation  Angle  and  Solar 

%  Array  Rotation  angle  as  a  function  of  time  and  as 

%  a  function  of  the  Sun  Rotation  Angle  in  order  to 

%  maintain  solar  array  tracking. 

a=17378, 
u=398600.5; 
t=0: 1:22798; 
x=(u/(aA3))A(l/2); 
B=(x)*t; 

figure(l) 

d=l*(pi/180); 

Psi=(atan((-sin(B)/tan(d))))*(180/pi); 

subplot(2,l,l); 

plot((B*(180/pi)),Psi); 

hold  on 

d=5*(pi/180); 

Psi=(atan((-sin(B)/tan(d))))*(180/pi); 

plot((B*(180/pi)),Psi), 

d=15*(pi/180); 

PsKatan((-sin(B)/tan(d))))*(  1 80/pi); 
plot((B*(180/pi)),Psi); 

d=30*(pi/180); 

Psi=(atan((-sin(B)/tan(d))))*(180/pi); 

plot((B*(180/pi)),Psi); 

d=45*(pi/180); 

Psi=(atan((-sin(B)/tan(d))))*(  1 80/pi), 
plot((B*(180/pi)),Psi), 

d=60*(pi/180); 

Psi=(atan((-sin(B)/tan(d))))*(  1 80/pi), 
plot((B*(180/pi)),Psi); 

d=75*(pi/180); 
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Psi=(atan((-sin(B)/tan(d))))*(180/pi), 
plot((B*(180/pi)),Psi); 

d=90*(pi/180), 

Psi=(atan((-sin(B)/tan(d))))*(  1 80/pi); 
plot((B*(180/pi)),Psi),grid; 
xlabel('Sun  Rotation  Angle  (degrees)'), 
ylabel('Yaw  Command  Angle  (degrees)'); 
titleCNominal  Yaw  Trajectories'); 
axis([0,360,-100,100]); 


gtext('d=l'); 

gtext('d=5'), 

gtext('d=15'); 

gtext('d=30'); 

gtext('d=45'); 

gtext('d=60'); 

gtext('d=75'); 

gtext('d=90'); 

hold  oflf 


subplot(2,l,2); 

d=l*(pi/180); 

Psi=(atan((-sin(B)/tan(d)))), 

n=atan(-(sin(Psi)./tan(B))); 

plot((B*(180/pi)),(n*(180/pi))); 

ylabel(' Solar  Array  Rotation  Angle  in  degrees'); 

xlabel('Sun  Rotation  Angle  in  Degrees'); 

title('Solar  Array  Trajectories'); 

hold  on 

d=5*(pi/180); 

Psi=(atan((-sin(B)/tan(d)))), 
n=atan(-(sin(Psi)./tan(B))); 
plot((B*(180/pi)),(n*(180/pi))); 

d=15*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))), 
plot((B*(180/pi)),(n*(l  80/pi))); 

d=30*(pi/180); 

Psi=(atan((-sin(B)/tan(d)))); 

n=atan(-(sin(Psi)./tan(B))); 
plot((B*(180/pi)),(n*(180/pi))); 
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d=45*(pi/180), 
Psi=(atan((-sin(B)/tan(d)))), 
n=atan(-(sin(Psi)./tan(B))); 
plot((B*(180/pi)),(n*(180/pi))); 

d=60*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))); 
plot((B*(180/pi)),(n*(180/pi))); 

d=75*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))), 
n=atan(-(sin(Psi)./tan(B))); 
plot((B*(180/pi)),(n*(180/pi))); 

d=90*(pi/180); 

Psi=(atan((-sin(B)/tan(d)))); 

n=atan(-(sin(Psi)./tan(B))); 

plot((B*(180/pi)),(n*(180/pi))),grid; 

axis([0,360,-100,100]); 

gtext('d=l'); 
gtext('d=5'); 

gtext('d=15'); 
gtext(*d=30'), 
gtext('d=45'); 
gtext('d=60'); 
gtextCd=75'); 
gtext(*d=90'); 

figure(2) 

t=t*(l/60);  %  Convert  time  to  minutes 

d=l*(pi/180); 

Psi=(atan((-sin(B)/tan(d))))*(  1 80/pi); 
subplot(2,l,l); 
plot(t,Psi); 

hold  on 
d=5*(pi/180); 

Psi=(atan((-sin(B)/tan(d))))*(  1 80/pi); 
plot(t,Psi); 

d=15*(pi/180); 
Psi=(atan((-sin(B)/tan(d))))*(  1 80/pi); 
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plot(t,Psi); 

d=30*(pi/180), 

Psi=(atan((-sin(B)/tan(d))))*(  1 80/pi); 
plot(t,Psi); 

d=45*(pi/180); 

Psi=(atan((-sin(B)/tan(d))))*(  1 80/pi), 
plot(t,Psi); 

d=60*(pi/180); 

Psi=(atan((-sin(B)/tan(d))))*(  1 80/pi), 
plot(t,Psi); 

d=75*(pi/180); 

Psi=(atan((-sin(B)/tan(d))))*(  1 80/pi), 
plot(t,Psi); 
axis([0,380,-100,100]); 

d=90*(pi/180); 

Psi=(atan((-sin(B)/tan(d))))*(180/pi), 

plot((B*(180/pi)),Psi),grid 

xlabel('Time  in  minutes'); 

ylabel('Yaw  Command  Angle  (degrees)'); 

titleCNominal  Yaw  Trajectories'); 

%gtext('d=l'); 

%gtext('d=5'); 

%gtext('d=15'); 

%gtext('d=30'); 

%gtext('d=45*); 

%gtext('d=60'); 

%gtext('d=75'); 

%gtext('d=90'); 

hold  off 

subplot(2,l,2); 

d=l*(pi/180); 

Psi=(atan((-sin(B)/tan(d)))); 

n=atan(-(sin(Psi)./tan(B))); 

plot(t,(n*(l  80/pi))); 

ylabel('Solar  Array  Rotation  Angle  in  degrees'), 

xlabel('Time  in  minutes'); 

title('Solar  Array  Trajectories'); 

hold  on 
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d=5*(pi/180); 

Psi=(atan((-sin(B)/tan(d)))), 
n=atan(-(sin(Psi)./tan(B))); 
plot(t,(n*(180/pi))),grid; 

d=15*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)  ./tan(B))); 
plot(t,(n*(180/pi))); 

d=30*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))), 
n=atan(-(sin(Psi)./tan(B))); 
plot(t,(n*(180/pi))); 

d=45*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))); 
plot(t,(n*(180/pi))); 

d=60*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))); 
plot(t,(n*(180/pi))); 

d=75*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))), 
n=atan(-(sin(Psi)./tan(B))); 
plot(t,(n*(180/pi))); 

d=90*(pi/180); 

Psi=(atan((-sin(B)/tan(d)))), 

n=atan(-(sin(Psi)./tan(B))); 

plot(t,(n*(180/pi))), 

axis([0,380,-100,100]), 

%gtext('d=l'); 

%gtext('d=5'); 

%gtext('d=15'); 

%gtext(*d=30'), 

%gtext('d=45'); 

%gtext(,d=60'); 

%gtext('d=75'); 

%gtext('d=90'); 
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figure(3) 

subplot(2,l,l); 
d=5*(pi/180); 

Psi=(atan((-sin(B)/tan(d))))*(  1 80/pi); 
plot(t,Psi); 

fori=l:length(Psi), 

Psi(i)=90; 
end; 

hold  on 

plot(t,Psi),  grid; 

xlabel(Time  in  minutes'); 

ylabel('Yaw  Command  Angle  (degrees)'); 

titleCNominal  Yaw  Trajectories'), 

%gtext('d=5'); 

%gtext('Yaw  Angle  =  Constant  90  degrees  for  d  <  5'); 

axis([0,380,-140,140]) 

y=[120,-120,120,-120]; 
q=[t(l),t(10695),t(12045),t(22799)]; 

subplot(2,l,2); 
plot(q,y) 

hold  on 

d=5*(pi/180); 

Psi=(atan((-sin(B)/tan(d)))); 

n=atan(-(sin(Psi)./tan(B))); 

plot(t,(n*(180/pi))),  grid; 

axis([0,380,-150,150]); 

xlabel('Time  in  minutes'); 

ylabel('Solar  Array  Rotation  Angle  in  degrees'); 

title('Solar  Array  Trajectory  for  d  <  5  degrees'); 

%gtext('d=5'); 
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APPENDIX  D 
MATLAB  CODE  TO  COMPUTE  DIRECTION  COSINE  MATRIX 

%  Danny  K.  Busch  4  October  1997 

%  Thesis 
% 

%  This  is  a  supplemental  .m  file  to  compute  the  initial  dcm 
%  from  the  orbit/nominal  frame  to  the  principal  axis  frame. 
%  I  am  making  the  assumption  that  the  orbit/nominal  frame 
%  is  lined  up  body  axes  when  this  starts.  So  this  is  the  initial 
%  dcm  from  both  the  body  axes  and  orbit  frame  to  the  principal 
%  axes. 

%  From  the  spreadsheet  we  have  the  Inertia  Matrix. 

Iold=[13337  -28.97  0,-28.97  14083  17.72;0  17.72  14113]; 

[evec,princ]=eig(Iold) 

%  Where  princ  is  the  matrix  containing  the  principal  moments 
%  of  inertias,  and  Evec  is  the  transpose  of  the  dcm  to  get  there. 

dcm=evec' 


%  To  prove  this,  multiply  dcm*Iold*dcm'  to  get  from  the  old 
%  (non-principal)  system  to  the  new  (principal)  system 

Inew=dcm*Iold*dcm' 
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%  Danny  K.  Busch  4  October  1997 

%  Thesis 

%  Plot  the  body  frame  vs  the  Principal  Inertia  Axis 

dcm=[0.9992,0.0388,0.0009;0.0354,-0.9026,0.4290,-0.0158,0.4287,0.9033] 
x=[0,0,0;dcm(l,l),dcm(l,2),dcm(l,3)]; 
y=[0,0,0;dmc(2, 1  ),dcm(2,2),dcm(2,3)]; 
z=[0,0,0;dcm(3,l),dcm(3,2),dcm(3,3)]; 

u=[0,0,0;l,0,0]; 
v=[0,0,0;0,l,0]; 

w=[0,0,O;0,0,l]; 

figure(l); 

plot3(x(:,l),x(:,2),x(:,3),V,y(:,l),y(:,2),y(:,3);r,,z(:,l),z(:,2),z(:,3);r,) 

hold  on 

plot3(u(:,l),u(:,2),u(:,3),V,v(:,l),v(:,2),v(:,3),V,w(:,l),w(:,2),wC,3),V) 

hold  off 

grid; 

xlabel(,X,),ylabelCY,),zlabel(,Z'); 

title('Body  Frame  vs  Principal  Inertial  Frame'), 
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APPENDIX  E 
SIMULATION  OF  THE  TOTAL  EXTERNAL  DISTURBANCES 

%  Danny  K.  Busch  27  July  1 997 

%  Thesis 

%  ADCS  Design 

% 

%  Gravity  Gradient  +  Solar  Pressure  +  Aerodynamic  Pressure  Torques 

% 

%  Simulates  the  attitude  motion  of  a  s/c  vs  time. 

%  That  is  the  orientation  of  the  b  frame  wrt  to  the 

%  orbit/nominal  frame.  Without  attitude  control. 

% 

function  xdot=total(x), 

%  Define  constants 

11=45.01961; 
12=35.53275; 
13=59.65711; 

Ky=(I3-I2)/Il; 
Kr=(I3-Il)/I2; 
Kp=(I2-Il)/I3; 

%  Semi-major  axis  of  sun/earth  system 
asun=1.49598e8; 

%  Sun  constant 
us=1.327ell, 
to=0; 
Mos=pi; 

%  Arguement  of  perigee 
w=0.4363323; 

%  Longitude  of  the  ascending  node 
omega=(pi/2); 

%  Inclination 
inc=(pi/4); 

%  Eccentricity  of  Sun/Earth  system 
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es=.  016708; 

%  Eccentricity  of  Earth/Satellite  system 
ee=0; 

ue=3.986012e5; 

%  Radius  of  the  Earth 
Re=6378; 

%  Generaic  cause  normalize  later,  vector  from  sun  to  earth 
r=10; 

%  R=dist  fin  center  of  earth  to  center  of  mass  of  s/c 

%  User  friendly  conversions 

wl=x(l); 

w2=x(2); 

w3=x(3); 

all=x(4); 

al2=x(5); 

al3=x(6); 

a21=x(7); 

a22=x(8); 

a23=x(9); 

a31=x(10); 

a32=x(ll); 

a33=x(12); 

rl=x(13); 

r2=x(14); 

r3=x(15); 

rdotl=x(16); 

rdot2=x(17); 

rdot3=x(18); 

time=x(19), 

R-sqrt(r  1  A2+r2A2+r3  A2); 
v=sqrt(rdot  1  A2+rdot2A2+rdot3  A2); 
nudot=v/R: 


%[eRs]h=r  vector  fm  sun  to  earth  expressed  in  heliocentric 

%frame 

%[eRs]I=r  vector  frn  sun  to  earth  expressed  in  earth 
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%centered  inertial 

%[eRs]o=  r  vector  frn  sun  to  earth  expressed  in  the 

%orbit  frame 

%[eRs]b=  r  vector  fin  sun  to  earth  expressed  in  the 

%body  frame 

%  Find  nu  and  d  as  a  function  of  time 

Moearth=0; 

toearth=0; 

%a  =  semimajor  axis  of  satellite  earth  system 

a=21000; 

ne=(ue/aA3)A(l/2); 

%  ne  mean  orbital  rate 

%M-Mo=n(t-to)  general  equation  mean  anomoly 

Me=ne*time; 

%nu  earth 

nu=Me  +  2*ee*sin(Me)  +  1.25*((ee)A2)*sin(2*Me); 

ns=sqrt(us/asunA3  ); 

%Ms-Mos=ns(time-to);  mean  anomoly  for  earth  sun  system; 

Ms=ns*time+pi; 

%orbital  parameter  nu  for  sun  earth  system 

d=Ms+2*es*sin(us)+1.25*(esA2)*sin(2*Ms); 

%Now  we  have  nu  and  d    or  nu  for  sun  and  earth 

%  R  vector  from  sun  to  earth  in  earth  centered  inertial 

eRsI=[-r*cos(d);-r*cos(.4084)*sin(d);-r*sin(.4084)*sin(d)], 

%  multiply  by  -1  so  we  R  points  from  earth  to  sun  according  to  Agrwal 

%  Define  dcm  from  earth  centered  inertial  to  rotating  orbital  frame 

J  =  [cos(nu)  sin(nu)  0;-sin(nu)  cos(nu)  0;0  0  1]; 

K  =  [cos(w)  sin(w)  0,-sin(w)  cos(w)  0;0  0  1]; 

L  =  [1  0  0;0  cos(inc)  sin(inc);0  -sin(inc)  cos(inc)]; 

Q  =  [cos(omega)  sin(omega)  0;-sin(omega)  cos(omega)  0,0  0  1]; 

oCi  =  J*K*L*Q, 

%  Now  transform  to  the  orbit  frame 
eRso  =  oCi*eRsI; 
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%define  dcm  from  orbit  to  body  frame 

bCo  =  [all  a21  a31;al2  a22  a32;al3  a23  a33]; 

%  Now  tranform  R  vector  to  orbit  frame 

eRsb  =  bCo*eRso; 

%  Now  compute  the  unit  vector  from  earth  to  sun 
Shatb  =  eRsb/(norm(eRsb)); 
s=  Shatb; 

%  Now  we  have  S  hat  expressed  in  the  b  frame  as  a 
%  function  of  d  and  nu,  which  is  a  function  of  time 

%  Define  n  hat  vectors  in  the  body  frame 


nl  =  | 

;0;1;0]; 

n2  = 

;0.866;  0.5;  0]; 

n3  =  | 

"0.866;  -0.5;  0]; 

n4  =  j 

;o,-i,o]; 

n5  =  | 

-0.866;  -0.5;  0], 

n6  =  | 

-0.866;  0.5;  0]; 

n7  = 

;0;1;0]; 

n8  =  | 

;0;-l;0]; 

n9  =  | 

;-l;0;0]; 

nlO  = 

[0;  0;  1]; 

nll  = 

[0;0;-l]; 

nl2  = 

[0;  0;  -1]; 

nl3  = 

[0,0;  1]; 

%  Now  define  the  dr  vectors  from  center  of  mass  to  center 
%  of  pressure 

drl  =  [0.11;  0.5;  0]; 
dr2  =  [0.543;  .25,  0]; 
dr3  =  [0.543; -.25;  0]; 
dr4  =  [0.11; -0.5;  0]; 
dr5  =  [-0.543;  -0.25;  0]; 
dr6  =  [-0.543;  0.25;  0]; 
dr7  =  [-0.756;  0.5;  0]; 
dr8  =  [-0.756;  -0.5;  0]; 
dr9  =  [-0.39;  0;  0]; 
drlO  =  [-.756;  0;  0.5]; 
drl  1  =  [-756;  0; -0.5]; 
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drl2  =  [0.11;0;-0.5]; 
drl3  =  [0.11;  0;  0.5]; 

Fe=1358; 
C  =  2.9979E8; 
P  =  Fe/C; 
A  =  zeros(l,13); 

%  F  =  PA(nhat  dot  shat)(2*nhat) 
%  Calculate  the  (nhat  dot  shat) 

hat  =  zeros(l,13); 
hat(l)  =  sum(nl.*s), 
ifhat(l)<0; 

hat(l)=0, 
else 

hat(l)=hat(l), 
end 

hat(2)  =  sum(n2.*s); 
ifhat(2)<0; 

hat(2)=0; 
else 

hat(2)=hat(2); 
end 

hat(3)  =  sum(n3.*s); 
ifhat(3)<0, 

hat(3)=0; 
else 

hat(3)=hat(3); 
end 

hat(4)  =  sum(n4.*s), 
ifhat(4)<0; 

hat(4)=0; 
else 

hat(4)=hat(4); 
end 

hat(5)  =  sum(n5.*s); 
ifhat(5)<0; 

hat(5)=0; 
else 

hat(5)=hat(5); 
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end 

hat(6)  =  sum(n6.*s); 
ifhat(6)<0; 

hat(6)=0; 
else 

hat(6)=hat(6); 
end 

hat(7)  =  sum(n7.*s); 
ifhat(7)<0, 

hat(7)=0; 
else 

hat(7)=hat(7); 
end 

hat(8)  =  sum(n8.*s); 
ifhat(8)<0; 

hat(8)=0; 
else 

hat(8)=hat(8); 
end 

hat(9)  =  sum(n9.*s); 
ifhat(9)<0; 

hat(9)=0; 
else 

hat(9)=hat(9); 
end 

hat(10)  =  sum(nlO.*s); 
ifhat(10)<0; 

hat(10)=0; 
else 

hat(10)=hat(10); 
end 

hat(ll)  =  sum(nll.*s); 
ifhat(ll)<0; 

hat(ll)=0; 
else 

hat(ll)=hat(ll); 
end 

hat(12)  =  sum(nl2.*s); 
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ifhat(12)<0; 

hat(12)=0; 
else 

hat(12)=hat(12); 
end 

hat(13)  =  sum(nl3.*s), 
ifhat(13)<0; 

hat(13)=0; 
else 

hat(13)=hat(13); 
end 

%  hat  =  nhat(i)  dot  shat(i) 
%  hat  =  [-o — 00--0--0] 

%  Define  the  area  matrix 

A  =  [.5  .5.5.5.5  .5.331  .331  1.77.331  .331  .724.724]; 

ps  =  .4, 
pd  =  .4; 

%  From  Agrwal's  Book  p.  134 

%  Ms  =  PA(n  dot  s)*r  X  [(l-ps)*s  +  2(ps  +  l/3pd)n] 

Ms  =  zeros(3,13); 

%  Assign  intermediate  variables 


G  =  zeros(3,13); 
G(:,l)  =  P*A(l)*hat(l)*drl; 
G(:,2)  =  P*A(2)*hat(2)*dr2; 
G(:,3)  =  P*A(3)*hat(3)*dr3; 
G(:,4)  =  P*A(4)*hat(4)*dr4; 
G(:,5)  =  P*A(5)*hat(5)*dr5; 
G(:,6)  =  P*A(6)*hat(6)*dr6; 
G(:,7)  =  P*A(7)*hat(7)*dr7; 
G(:,8)  =  P*A(8)*hat(8)*dr8, 
G(:,9)  =  P*A(9)*hat(9)*dr9; 
G(:,10)  =  P*A(10)*hat(10)*drlO; 
G(:,l  1)  =  P*A(1  l)*hat(l  l)*drl  1; 
G(:,12)  =  P*A(12)*hat(12)*drl2; 
G(:,13)  =  P*A(13)*hat(13)*drl3; 


T  =  zeros(3,13); 
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T(:,l)  =  ((l-ps)*s)  +  ((2*(ps  +  0.333*pd))*nl); 
T(:,2)  =  ((l-ps)*s)  +  ((2*(ps  +  0.333*pd))*n2), 
T(:,3)  =  ((l-ps)*s)  +  ((2*(ps  +  0.333  *pd))*n3); 
T(:,4)  =  ((l-ps)*s)  +  ((2*(ps  +  0.333  *pd))*n4), 
T(:,5)  =  ((l-ps)*s)  +  ((2*(ps  +  0.333*pd))*n5); 
T(:,6)  =  ((l-ps)*s)  +  ((2*(ps  +  0.333  *pd))*n6), 
T(:,7)  =  ((l-ps)*s)  +  ((2*(ps  +  0.333*pd))*n7); 
T(:,8)  =  ((l-ps)*s)  +  ((2*(ps  +  0.333*pd))*n8), 
T(:,9)  =  ((l-ps)*s)  +  ((2*(ps  +  0.333*pd))*n9); 
T(:,10)  =  ((l-ps)*s)  +  ((2*(ps  +  0.333*pd))*nl0), 
T(:,l  1)  =  ((l-ps)*s)  +  ((2*(ps  +  0.333*pd))*nl  1); 
T(:,12)  =  ((l-ps)*s)  +  ((2*(ps  +  0.333*pd))*nl2); 
T(:,13)  =  ((l-ps)*s)  +  ((2*(ps  +  0.333*pd))*nl3); 

Z  =  cross(G,T), 

%  We  now  have  a  3  by  13  matrix  containing  torques  in  bl  direction  for 
%  thirteen  areas  in  the  first  row,  etc ... 


%  Solar  pressure  moments  in  b  direction 


Mspl  =  sum(Z(l, 
Msp2  =  sum(Z(2, 
Msp3  =  sum(Z(3, 


)); 
)); 
)); 


%  (l/2)*p*Cd*A*v*v*(uvhat  X  Icp) 
%  Define  p  as  a  function  of  R 

ifR<6578; 

p=4e-10; 
elseifR<6678, 

p=5e-ll; 
elseif  R  <  6778; 

p=1.5e-ll; 
elseif  R  <  6878; 

p=5e-12; 
elseif  R  <  6978 

p=2e-12 
else 

p=8e-13, 
end 

%  Define  Cd  as  a  constant  given  on  P. 90  of  handout 

Cd=3.2; 


164 


%  Define  Area  constants  in  the  bl,  b2  and  b3  directions.  This  is  the  total 
%  cross  sectional  area  in  the  bl  hat  direction,... b2  hat  direction  and  b3 
%  hat  direction. 

Abl=1.77; 
Ab2=884; 
Ab3=1.61; 

Abody=[Abl  Ab2  Ab3]'; 

%  Now  build  the  velocity  vector  for  the  I  frame 

RDOTi=[rdotl  rdot2  rdot3]'; 

%  Now  convert  velocity  from  I  frame  to  B  frame 

RDOTb=bCo*oCi*RDOTi; 

%  Now  find  the  unit  vector  of  RDOTb,  uvhat. 

uvhatb=RDOTb./norm(RDOTb); 

%  Now  compute  area  perpendicular  to  uvhat. 

AT=abs(Abody  *uvhatb); 

%  Compute  total  area 

ATOTAL=sum(AT); 

%  Now  compute  (uvhatb  X  Icp) 

Icp=zeros(3,13); 

Icp=[drl  dr2  dr3  dr4  dr5  dr6  dr7  dr8  dr9  drlO  drl  1  drl2  drl3]; 

Maero=(l/2)*p*Cd*ATOTAL*v*v*(cross(uvhatb,Icp)); 

Maerototal=[sum(Maero(l,:))  sum(Maero(2,:))  sum(Maero(3, :))]', 


Mspl=sum(Z(l, 
Msp2=sum(Z(2, 
Msp3=sum(Z(3, 


)); 
)); 
)); 


Mlgg=((3*ue)/RA3)*(I3-I2)*al2*al3; 
M2gg=((3*ue)/RA3)*(I1-I3)*al3*all; 
M3gg=((3*ue)/RA3)*(I2-Il)*all*al2; 
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xdot(l)=  ((Mspl+Maerototal(l)+Mlgg)/Il)  -  Ky*w2*w3; 

xdot(2)=  ((Msp2+Maerototal(2)+M2gg)/I2)  -  Kr*w3*wl; 

xdot(3)=  ((Msp3+Maerototal(3)+M3gg)/I3)  -  Kp*wl  *w2; 

xdot(4)=al2*(w3-a33*nudot)-al3*(w2-a32*nudot); 

xdot(5)=al3*(wl-a31*nudot)-all*(w3-a33*nudot), 

xdot(6)=al  1  *(w2-a32*nudot)-al2*(wl-a3 1  *nudot); 

xdot(7)=a22*(w3-a33*nudot)-a23*(w2-a32*nudot); 

xdot(8)=a23  *(wl-a3 1  *nudot)-a2 1  *(w3-a33  *nudot); 

xdot(9)=a2 1  *(w2-a32*nudot)-a22*(wl  -a3 1  *nudot); 

xdot(10)=a32*(w3-a33*nudot)-a33*(w2-a32*nudot), 

xdot(l  I)=a33*(wl-a3 1  *nudot)-a3 1  *(w3-a33*nudot); 

xdot(12)=a3 1  *(w2-a32*nudot)-a32*(wl-a3 1  *nudot), 

xdot(13)=x(16), 

xdot(14)=x(17), 

xdot(15)=x(18); 

xdot(16)=((-ue)/RA3)*rl; 

xdot(17)=((-ue)/RA3)*r2; 

xdot(18)=((-ue)/RA3)*r3; 

xdot(  1 9)=Msp  1  +Maerototal(  1  )+M  1  gg; 

xdot(20)=Msp2+Maerototal(2)+M2gg, 

xdot(2 1  )=Msp3+Maerototal(3)+M3gg; 
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APPENDIX  F 


SMALL  ANGLE  SIMULATION 


Pitch  Response  Simulation 


H3* K     time    ] 

Clock  To  Workspace2 


TQpitch     To  Workspace3 


ED 

Steplmput 


Sine  Wave 


impulseb 


rad  to  degrees     To  Workspace 


K    »p   3 

To  Workspace? 


Impulse 
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%  Danny  K.  Busch  4  October  1997 

%  Thesis 

%  This  is  the  initial  conditions  used  with  the  simulink. 

wo=2.75578e-4; 

Ixx=13337, 

Iyy=  14083, 

Izz=14113; 

Kp=3; 

Tp=433; 

Ky=l; 

Ty=751; 

Kr=5; 

Tr=2060: 
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%  Danny  Busch  4  October  1997 

%  Thesis 

%  This  file  plots  torque,  angle  (theta),  and  reaction 

%  wheel  momentum  vs  time. 

figure(l) 

subplot(3,l,l); 
plot(time,  TQpitch),  grid, 
title('Pitch  Torque  vs  Time'); 
xlabel('Time  (sec)'), 
ylabel('Torque  (N-m)'); 
orient  tall 

subplot(3,l,2); 

plot(time,  theta),  grid, 

title('Pitch  Angle  vs  Time'); 

xlabel('Time  (sec)'); 

ylabel('Pitch  Angular  Response  (deg)'); 

orient  tall 

subplot(3,l,3) 

plot(time,hp),  grid; 

title('Pitch  Wheel  Angular  Momentum  vs  Time'); 

xlabel('Time  (sec)'); 

ylabel(' Angular  Momentum  (N-m-s)'), 

orient  tall 
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Roll  Yaw  Response  Simulation 


TOvll 


ED 

Steplmput 


To  Workspace3 


Clock  To  Workspace2 


Sing  Wave 
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impulseb 


Impulse 
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Transfer  Fen 
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Sum 


PID 


PID  Controller 
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Gam 


Integrator 
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►y/y* ►£    "y 
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%  Danny  K.  Busch  4  October  1 997 

%  Thesis 

%  Plot  m  file  used  for  thesis  to  aid  in  plotting 

%  the  results  from  Roll  and  Yaw 

%  Plots  torque,  angle  and  momentum  vs  time 

figure(l) 

subplot(3,l,l); 

plot(time,TQroll),  grid; 

title('Roll  Torque  vs  Time'), 

xlabel('Time  (sec)1); 

ylabel(Torque  (N-m)'); 

orient  tall 

subplot(3,l,2), 

plot(time,phi),  grid; 

title('Roll  angle  phi'), 

xlabel('Time  (sec)'), 

ylabel('Roll  Angular  Response  (deg)'); 

orient  tall 

subplot(3,l,3) 

plot(time,hr),  grid; 

title('Angular  Momentum  vs  Time'), 

xlabel('Time  (sec)'); 

ylabel(' Angular  Momentum  (N-m-s)'), 

orient  tall 

figure(2) 

subplot(3,l,l); 

plot(time,TQyaw),  grid; 

title('Yaw  Torque  vs  Time'); 

xlabel('Time  (sec)'); 

ylabel('Torque  (N-m)'); 

orient  tall 

subplot(3,l,2); 

plot(time,psi),  grid; 

title('Yaw  angle  psi'); 

xlabel('Time  (sec)'); 

ylabel('Yaw  Angular  Response  (deg)'); 

orient  tall 

subplot(3,l,3) 

plot(time,hy),  grid; 

title(' Angular  Momentum  vs  Time'); 

xlabel('Time  (sec)'); 

ylabel(' Angular  Momentum  (N-m-s)'); 

orient  tall 


173 


This  Page  Intentionally  Left  Blank 


174 


APPENDIX  G 
LARGE  ANGLE  SIMULATION 

%  Danny  K.  Busch  3  October  1997 

%  Thesis 

% 

%  Using  S  Function  to  model  the  dynamics  of  the  Spacecraft 

%  Define  S-Function 

function[sys,  xo]=large(t,x,u,flag) 

%  VDPM  EOM  as  an  M-file 

Ixx=13337; 
Iyy= 14083; 
Izz=14113, 
wo=2.75578e-4; 

ifabs(flag)=l, 

%  Return  rates 

sys(l)=x(2); 

sys(2)=(-l/Ixx)*(x(4)*cos(x(l))*Izz*x(6)+x(6)*x(6)*cos(x(3))*sin(x(l))*Izz- 
wo*cos(x(l))*cos(x(5))*Izz*x(6)- 

wo*sin(x(l))*sin(x(3))*sin(x(5))*Izz*x(6)+x(4)*x(4)*sin(x(l))*Iyy- 
x(6)*cos(x(3))*cos(x(l))*Iyy*x(4)- 

wo*sin(x(l))*cos(x(5))*Iyy*x(4)+wo*cos(x(l))*sin(x(3))*sin(x(5))*Iyy*x(4)+u(4)+x(4) 
*cos(x(l))*u(3)+x(6)*cos(x(3))*sin(x(l))*u(3)-wo*cos(x(l))*cos(x(5))*u(3)- 
wo*sin(x(l))*sin(x(3))*sin(x(5))*u(3)+x(4)*sin(x(l))*u(2)- 
x(6)*cos(x(3))*cos(x(l))*u(2)- 

wo*sin(x(l))*cos(x(5))*u(2)+wo*cos(x(l))*sin(x(3))*sin(x(5))*u(2)- 
3  *wo*wo*sin(x(  1  ))*cos(x(  1  ))*cos(x(3))*cos(x(3))*Izz+3  *wo*wo*sin(x(  1  ))*cos(x(  1  ))*c 
os(x(3))*cos(x(3))*Iyy); 

sys(3)=x(4); 

sys(4)=(-l/Iyy)*(- 
x(4)*sin(x(l))*Ixx*x(2)+x(6)*cos(x(3))*cos(x(l))*Ixx*x(2)+wo*sin(x(l))*cos(x(5))*Ixx 
*x(2)-wo*cos(x(l))*sin(x(3))*sin(x(5))*Ixx*x(2)- 

x(2)*Izz*x(6)+x(6)*x(6)*sin(x(3))*Izz+wo*cos(x(3))*sin(x(5))*Izz*x(6)- 
x(4)*sm(x(l))*u(l)+x(6)*cos(x(3))*cos(x(l))*u(l)+wo*sin(x(l))*cos(x(5))*u(l)- 
wo*cos(x(l))*sin(x(3))*sin(x(5))*u(l)- 
x(2)*u(3)+x(6)*sin(x(3))*u(3)+wo*cos(x(3))*sin(x(5))*u(3)+u(5)- 


175 


3*wo*wo*sin(x(3))*cos(x(3))*cos(x(l))*Izz+3*wo*wo*sin(x(3))*cos(x(3))*cos(x(l))*I 
xx); 

sys(5)=x(6); 

sys(6)=(-l/Izz)*(x(2)*Iyy*x(4)-x(6)*sin(x(3))*Iyy*x(4)- 
wo*cos(x(3))*sin(x(5))*Iyy*x(4)-x(4)*cos(x(l))*Ixx*x(2)- 

x(6)*cos(x(3))*sin(x(l))*Ixx*x(2)+wo*cos(x(l))*cos(x(5))*Ixx*x(2)+wo*sin(x(l))*sin( 
x(3))*sin(x(5))*Ixx*x(2)+x(2)*u(2)-x(6)*sin(x(3))*u(2)-wo*cos(x(3))*sin(x(5))*u(2)- 
x(4)*cos(x(l))*u(l)- 

x(6)*cos(x(3))*sin(x(l))*u(l)+wo*cos(x(l))*cos(x(5))*u(l)+wo*sin(x(l))*sin(x(3))*sin( 
x(5))*u(l)+u(6)- 
3*wo*wo*sin(x(3))*cos(x(3))*sin(x(l))*Ixx+3*wo*wo*sin(x(3))*cos(x(3))*sin(x(l))*Iy 

y); 

elseifflag=0 

%  Return  initial  conditions 

sys=[6;0;6;6;0;0]; 

%  [continuous  states,  discrete  states,  outputs,  inputs,  discontinuous  roots, 
%  direct  feed  through] 

xo=[0;0;0;0;0;0]; 

elseifflag==3 

%  If  flag  =  3,  return  system  outputs 

sys(l)=x(l); 
sys(2)=x(2); 
sys(3)=x(3); 
sys(4)=x(4); 
sys(5)=x(5); 
sys(6)=x(6); 


else 
end; 


sys=[]; 
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%  Danny  K.  Busch  4  October  1997 

%  Thesis 

%  This  is  the  initial  conditions  used  with  the  simulink. 

wo=2.75578e-4; 
Ixx=13337; 
Iyy=14083; 
Izz=14113; 

Kr=l; 
Kp=2.0; 

Ky=2; 


%Kr=0.5; 
Tr=2060; 

%Kp=3; 
Tp=433; 

%Ky=.l; 
Ty=751; 
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%  Danny  Busch  1 3  October  1 997 

%  Thesis 

% 

%  m  file  used  to  simulate  the  motion  of  the  satellite 

function  act=psicmd(t) 

act=(-1.000507)*atan(11.43005*sin((0.0002755937)*(t))); 

end; 
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